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The spectroscopical dissociation scheme of CO gives for the heat of sublimation of carbon 
into §S carbon atoms L,=169,735 kcal./mole. This means an 84.58-kcal./mole average bond 
energy in diamond. The C—H bond energy of the methylene radical is 97.02 kcal./mole, the 
strength of the paraffinic C—C bond 82.59 kcal./mole, the olefinic double bond C=C is 
145.12 kcal./mole. The rupture of successive bonds in methane requires 100.76, 97.36, 97.02, 


97.02 kcal. /mole. 


INTRODUCTION 


HE determination of bond energies of 

organic compounds presupposes in most 
cases the knowledge of the heat of sublimation 
of carbon. This value can be given with great 
accuracy by the aid of the spectroscopical dis- 
sociation scheme of the CO molecule, as deter- 
mined by Schmid and Geré,"? and enables one to 
determine important quantities of organic chem- 
istry. An attempt to explain thermochemical 
data by means of the corresponding heat of 
sublimation value of carbon was made by 
Schmid* directly after the determination of the 
dissociation energy of CO, while Schmid and Geré 
were able to interpret the results of electronic 
collision experiments in organic compounds.' 


* This paper is written by J. G. Valatin on the basis of 
the notes left by Geré, and of the memory of some dis- 
cussions. Most of these notes are dated from 1944, and 
they had had to form the bases of a larger work. Some of 
the numerical values have been recalculated according to 
the more recently available thermochemical data, and 
these are on the responsibility of the compilator of the 
article. Results concerning some carbonyl compounds and 
a critical discussion of the more recent literature are 
contained in the following two articles. 

1R. Schmid and L. Gerd, Zeits. f. physik. Chemie (B) 
36, 105 (1937). 

2 R. Schmid and L. Gerd, Physik. Zeits. 39, 460 (1938). 

3R. Schmid, Math. és Természettudomanyi Ertesité 
54, 769 (1936). 

*R. Schmid and L. Gerd, Math. és Természettudomanyi 
Ertésité 56, 865 (1937); 57, 637 (1938). 


In view of the different opinions with respect 
to the value of the heat of sublimation of carbon, 
however, the values of bond energies of organic 
compounds given in the literature are very 
divergent. It seemed, therefore, worth while to 
review the subject again, and this paper contains 
part of this work. In this way, very satisfactory 
values result with respect to the C—C and C—H 
bond energies of hydrocarbons, which are also in 
agreement with the data relating to other organic 
compounds. In regard to aliphatic hydrocarbons 
an additivity of bond energies is found within a 
few percent, while the more recent thermo- 
chemical measurements make it possible to follow 
also the finer deviations from bond additivity. 
The large resonance energies supposed by some 
authors do not occur. 

In order to obtain data characterizing a single 
molecule, the AHf® values of heat of formation 
at 0°K are used. The theoretical values char- 
acterizing bond strengths refer to a vibrationless 
state and have to be calculated by subtracting 
half the sum of all vibrational frequencies of the 
molecular ground state. It is to be noted that 
calculations with these theoretical bond strengths 
relating to vibrationless states show an improve- 
ment of bond additivity, especially in compounds 
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containing also other elements besides carbon 
and hydrogen. 
HEAT OF SUBLIMATION OF CARBON 
In identifying the energy height of the ob- 
served predissociation at 89620+47 cm above 
the X'S (v=0, J=0) ground state in the CO 
spectrum with the energy height of the C(5S) 


+O(*P») atomic term combination," the height 
of the atomic term combination C(S)+O(*P2) 
is at 89400 cm. The dissociation energy of O» 
resulting from band spectra is 40990 cm-. Using 
a conversion factor 1 cm~'=2.8581 cal./mole, 
there results by means of these data and the heat 
of formation® of CO at 0°K: 


CO(gas)—>C (5S) +O(@P2) — 255.514 kcal./mole, 
O(?P:2)—302(gas)+ 58.577, 
C(c, graphite) + 402(gas)>CO(gas)+ 27.202, 


C(c, graphite) —>C (gas) — 169.735 kcal. /mole. 


The heat of sublimation of carbon into °S 
atoms is accordingly L.=169.735 kcal./mole at 
0°K. This is the value which plays the important 
part mentioned in the calculation of bond ener- 
gies of organic compounds, these bond energies 
being characterized in general by energy values 
counted from the °S energy level of carbon, the 
lowest energy state of the sp* configuration of 
tetravalent carbon being to be taken into 
account as final dissociation product. 

In CO carbon is divalent and the levels of the 
CO(X'Z) ground state converge according to the 
dissociation scheme quoted to the C(!D)+O('D) 
atomic term combination, the bond energy of 
the CO molecule being 233.39 kcal./mole. The 
dissociation energy of CO counted from the 
lowest C(@P)+O(P) atomic term combination 
is 158.91 kcal./mole. 

Using for the heat of transition of graphite 
into diamond*® AHf®=0.577 kcal./mole, for the 
heat of sublimation of diamond into °S carbon 
atoms at 0°K there results 


C(c, diamond)—>C(8S) (gas) — 169.158 kcal./mole. 


The average bond strength of the carbon-carbon 
linkage in diamond is 


[C—C Jaiamona = 169.158/2 = 84.58 kcal./mole. 


BOND ADDITIVITY AND DEPARTURES FROM 
BOND ADDITIVITY IN THE PARAFFIN 
HYDROCARBONS 


The heats of formation of the normal paraffins® 
show a regular increase of 3.673 kcal./mole per 


5D. D. Wagman, J. E. Kilpatrick, W. J. Taylor, K. S, 
Pitzer, and F. D. Rossini, Bur. Stand. J. Research 34, 
143 (1945). 

6E. J. Prosen, K. S. Pitzer, and F. D. Rossini, Bur. 
Stand. J. Research 34, 403 (1945). 


methylene radical at 0°K beginning from hexane. 
The variation of the heats of formation of the 
first members of the series is slightly different, 
the differences probably to be attributed first of 
all to the effect of the two methyl radicals at the 
ends of the chain which interact also with the 
second next methylene radicals. Apart from this, 
however, there is an apparent bond additivity. 

The knowledge of the value of the heat of 
sublimation of carbon makes it possible to 
determine the bond strengths. Slight deviations 
are found from the average bond strength in 
isomeric paraffin hydrocarbons, which yield also 
a measure for the departures from bond addi- 
tivity. This can be illustrated by an investigation 
of the interaction of methyl radicals on the 
same carbon atom in isomeric paraffins.’ 

Forming the differences between the heats of 
formation of 3-methylpentane and n-hexane, and 
those of 3-methylhexane and n-heptane, there 
results in average a difference of 0.34 kcal./mole. 
Approximately, this can be interpreted as the 
difference by which the sum of the hydrogen 
bond energies in a methyl radical exceeds three 
times the [C—H ]cue hydrogen bond energy in 
the methylene radical. With this supposition, 
the total bond energy of an m-paraffin corre- 
sponding to a dissociation into C(*S) and H(?S) 
atoms is 


(n—1)[C— Ch. paraffin + (2n+2)[(C —H Jone 
+0.68 kcal./mole. 


Adding the value 0.34 kcal./mole to the differ- 


7 The original calculations were performed without data 
for the isomeric octanes. Analogous calculations performed 
with respect to the octanes with slight differences give the 
same results. 
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ence of heat of formation of 2-methylhexane and 
n-heptane gives 1.46 kcal./mole as the excess 
energy of two CH; radicals at the end of a 
paraffin chain; adding 0.34 kcal./mole to the 
difference of heat of formation of 2,2-dimethyl- 
pentane and n-heptane gives 3.69 kcal./mole as 
the excess energy of three CH; radicals at the 
end of the chain. The selected figures are rather 
characteristic, and the average excess energy of 
a methyl radical is in this way roughly 0.3, 0.7, 
1.2 kcal./mole in the case of 1, 2, or 3 methyl 
radicals at the end of a paraffin chain. This can 
be considered also as a measure of the departure 
from bond additivity in the paraffins. 


BOND ENERGIES OF THE PARAFFIN AND 
MONOOLEFIN HYDROCARBONS 


In order to calculate the [C—H ]cne and the 
average [C—C] bond energies in the paraffin 


hydrocarbons, we can use the —30.98-kcal. /mole 


heat of formation value of m-hexane, the higher 
n-paraffins yielding the same values. With 
L.=169.735 kcal./mole and the spectroscopical 
dissociation energy of He, D(H2) =36115 
= 103.22 kcal./mole there results for n-hexane, 


n-hexane—6C (°S) + 14H (?S) 
—1771.93 kcal./mole. 


From the mentioned increase per methylene 
radical in the paraffin series 


—H —C] = 169.735+103.22 
+3.673 =276.63 kcal./mole. (a) 


Subtracting five times this equation from the 
corresponding equation for n-hexane 


5[C—C]+14[C —H 
= 1771.93 kcal./mole 


gives 4[C — H = 388.10 kcal./mole, that is, 
[C—H ]cu,=97.02 kcal./mole. 


The mean paraffinic C—C bond energy is from 
Eq. (a) 
[C—C]=82.59 kcal./mole. 


As to the double bond of the monoolefins, the 


—2.69 kcal./mole heat of formation value® for 


8 J. E. Kilpatrick, E. J. Prosen, K. S. Pitzer, and F, D. 
Rossini, Bur. Stand. J. Research 36, 559 (1946). 


TABLE I. Bond energies in methane. 


Heat of 

Energy formation 

required into C(5S) 

to rupture and H(?S) 
the bond atoms 

kcal./mole kcal./mole 
First C—H bond 97.02 CH 97.02 
Second C—H bond 97.02 CH: 194.04 
Third C—H bond 97.36 CH; 291.40 
Fourth C—H bond 100.76 CH, 392.16 


1-hexene gives for the difference in bond energies 
of n-hexane and 1-hexene 


2(C—H 
+ 30.98 — 2.69 = 131.51 kcal./mole, 


the same difference resulting also for the higher 
members of the two series. This relation com- 
bined with Eq. (a) gives for the double bond 


[C=C]=145.12 kcal./mole. 


METHANE 


With the given values of Z., D(H), and the 
15.987-kcal./mole heat of formation of methane 
there results 


CH.>C (5S) +4H(°S) — 392.16 kcal./mole, 


involving a mean bond energy of 392.16/4 = 98.04 
kcal./mole in methane. According to the pre- 
ceding the mean hydrogenic bond energy of the 
methylene radical is 97.02 kcal./mole, the excess 
energy in the methyl radical being 0.34 kcal./ 
mole. This means for the methyl radical 


CH;—C('S) +3H(?S) — 291.40 kcal./mole, 


and the difference of the bond energies of 
methane and the methyl radical, the energy 
required to rupture the fourth hydrogenic bond, 
is 


CH,—CH;+H — 100.76 kcal. /mole.® 


Table I shows the energies required to rupture 
the different bonds in methane, supposing that 
the bond energy of the CH radical is identical 
with [C—H Jcny. 

®The striking agreement with the results of G. B. 
Kistiakowsky and E. R. Van Artsdalen, J. Chem. Phys. 
12, 469 (1944) is to be noted. They obtained for this value 


100.8+1 kcal./mole at 0°K by direct reaction kinetical 
measurements. 
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Dissociation Schemes of Some Carbonyl Compounds 
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University for Technical and Economic Sciences, Budapest, Hungary 


(Received February 9, 1948) 


Dissociation schemes and bond energies are given for formaldehyde, acetaldehyde, acetone, 
and ketene. The distinction between dissociation term combinations which belong or do not 
belong to the lowest electronic state of the molecule is stressed, only energies corresponding 
to the former being in connection with bond energies. The survey of the different energetic 
possibilities may contribute to a better understanding of photo-chemical processes. 


N the preceding article of Geré' the bond 
energies of some hydrocarbons were given by 
the aid of the knowledge of the heat of sublima- 
tion value of carbon. These bond energies, 
together with some spectroscopical and thermo- 
chemical data, make it possible to determine the 
dissociation schemes of some carbonyl com- 
pounds, and enable a further examination of 
observed predissociations, dissociation continua, 
and photo-chemical effects. The investigation of 
these dissociation schemes allows also a survey 
of the energetic possibilities which have to be 
considered in connection with the corresponding 
photo-chemical processes and may contribute to 
the understanding of these processes. . 
The skeleton of the given dissociation schemes 
is the work of Geré. A part of the original 
interpretations have had, however, to be changed 
owing to the since available data, and the 
schemes have been further completed. The re- 


TABLE I. 

Process Q kcal./mole 
96.60 
C(graphite)—C('S) 169.73 
58.58 
H.—2H(2S) 103.22 
CO(X!Z)—C (5S) +O0(P) 255.51 
CO(a*II)—>C (5S) +O(P) 116.97 
CO(a#z)—+C (5S) +O(P) 97.23 
=CH—C(5S) ++H(?S) 97.048 
= (5S) +2H(?S) 194.08" 
—CH;—C(5S) 291.398 
CH,—C (8S) +4H(?S) 392.16 
(5S) +6H (2S) 665.51» 
(5S) +4H(?S) 531.39¢ 


¥ fr The slightly modified values are given in the following paper of 
atin. 
b With the AH/° value of E. J. Prosen, K. S. Pitzer, and F. D. Rossini, 
Bur. Stand. J. Research 34, 403 (1945). 
° With the AH/® value of J. E. Kilpatrick, E. J. Prosen, K. S.Pitzer, 
and F. D. Rossini, Bur. Stand. J. Research 36, 559 (1946). 


1L, Gerd, J. Chem, Phys. 16, 1011 (1948). 


sponsibility for the greater part of the expressed 
views devolves, accordingly, on Valatin. 


FORMALDEHYDE 


In order to obtain data characterizing a single 
molecule which can be directly compared with 
the results of spectroscopical and photo-chemical 
observations, heat of formation values at 0°K 
are to be used throughout, and every numerical 
value is given at this temperature. For the heat 
of formation of gaseous formaldehyde we use the 
value adopted by Long and Norrish? as a “‘best 
value” from data of Delépine and Badoche, and 
of Wartenberg and Lerner-Steinberg, 


H,CO(gas) = C(graphite) +H2+ 302 
— 25.3 kcal./mole, 


which yields by means of the data given in 
Table I. 


H.CO(gas) 2H (??S) + C(®S) +O(@P) 
— 356.8 kcal./mole. 


The greater part of the data of Table I is given 
or quoted in the preceding paper of Geré. 

Using the data of Table I and this 356.8- 
kcal./mole value for H2CO, the dissociation 
scheme of Table II is obtained, in which the 
energy heights of the different dissociation 
products are counted from the ground level of 
the HeCO molecule. As is seen from the dis- 
sociation process (6), the bond energy required 
to remove the oxygen atom is 162.7 kcal./mole. 
Supposing from the similarity of the vibrational 
frequencies that the ground state of the for- 
maldehyde molecule is connected with the a*II- 
state of CO, the dissociation process (7) is 


2L. H. Long and R. G. W. Norrish, Proc. Roy. Soc. 
A187, 337 (1946). 
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TABLE II. Dissociation products and dissociation energies TABLE III. Dissociation products and dissociation energies 


of formaldehyde. 


of acetaldehyde. 


Energy height 


Energy height 


Dissociation products kcal./mole Dissociation products kcal./mole 

1 H2+CO(X!2) -1.9 1 CH,+CO(X!3) —1.5 
2 H(?S)+HCO ~75-85 2 —CH;+HCO = 73-83 
3 2H(?S)+CO(X!3) 101.3 3 H(?S)+CH;CO 
4 136.6 4 87.7 
5 H.+CO(a*z) 156.4 5 =CH2+H2+CO(x!z) 93.4 
6 =CH2+0(P) 162.7 6 =CH2+H2CO 95.3 
7 239.8 7 H(?S)+CH.2CHO ~98 
8 —CH;+H(?S)+CO(X!z) 99.3 

9 CH,+CO(a'Il) 137.0 

10 CHy+CO(az) 156.8 
characteristic of the bond strengths of the 11 =CHCH;+0(P) ~175.1 
12 237.8 


hydrogen atoms, and for the mean bond energy 
of a hydrogen atom the result is 119.9 kcal. /mole. 
This figure contains, however, also the inter- 
action energy of the hydrogen atoms with the 
oxygen atom. The role of the CO(a*II) state in 
formaldehyde is meant in the sense that the 
vibrational levels of the formaldehyde ground 
state converge, directly or by a successive con- 
vergence of the levels of the H+CHO term 
combinations, to the energy height of process (7), 
the corresponding term combination giving also 
the asymptotic energy value of the potential 
function of the formaldehyde ground state.* 
Accordingly, it is evident that one cannot 
estimate bond energies from the energy height 
of the dissociation products (3). In view of the 
opinions found in the literature, it is to be 
emphasized that the energy heights of pre- 
dissociation continua which belong to different 
electronic states of the molecule, are, in general, 
not characteristic of the properties of the mo- 
lecular ground state, and cannot be brought into 
connection with the concept of bond energies. 
On the other hand, the energy heights of pre- 
dissociations play an important part in chemical 
reactions, and dissociation continua can also get 
a role by direct transition in the primary absorp- 
tion process of a photo-chemical decomposition. 
The energy height of the dissociation products 
(2) which play an important part in the pho- 
tolysis of formaldehyde has been estimated from 
the energy height (3), by means of the 26- 
kcal./mole and 16-kcal./mole activation energies 
of the reaction HCO—H(?S)+CO(X!) given 
by Gorin? and by Style and Summers.‘ As our 
* See also J. G. Valatin, Comptes rendus 226, 702 (1948). 
3E. Gorin, J. Chem. Phys. 7, 256 (1939). 


*D. W. G. Style and D. Summers, Trans, Faraday Soc. 
42, 388 (1946). | 


knowledge with respect to the HCO radical is 
very little, it cannot be decided whether the low 
~75-85-kcal./mole energy value of the dissocia- 
tion process (2) as compared with the formerly 
calculated 119.9-kcal./mole mean hydrogen bond 
energy means merely a strong reorganization 
energy of the HCO radical—with the opposite 
sign as in methane where the energy required to 
remove the first hydrogen atom is 100.8 kcal./ 
mole in spite of a mean bond energy of 98 
kcal./mole—or whether the ground state of for- 
maldehyde corresponds to an excited energy state 
of HCO, in which case process (2) is not charac- 
teristic to the ground state of H2CO, in the same 
way as process (3) is not. In any case, it is possible 
that HCO has an excited energy state which 
should be included in Table II, and which might 
play a role in the photolysis of formaldehyde. 
From the point of view of the formaldehyde 
photolysis observed between 3650—2600A, the 
processes (1), (2), and (3) of Table II have to be 
taken into account. The possibility of process 
(3) at shorter wave-lengths is to be emphasized, 
since this process is generally left out of con- 
sideration. The predissociation at 2750A ob- 
served by Price® corresponds also to these dis- 
sociation products, though this predissociation 
limit is often attributed to the removal of a 
single hydrogen atom.® The predissociation ob-. 
served at 1745A, corresponding to the “ex- 
tremely diffuse band’ of Price, could be at- 
tributed, according to Table II, either to the 
term combination H:+CO(a*z) or to 


5 W. C. Price, J. Chem. Phys. 3, 256 (1935). 
a : om” for instance, M, Burton, J. Chem, Phys. 7, 1072 
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+O(*P), the corresponding term combinations 
occurring also in the dissociation schemes of 
acetaldehyde and of acetone, in the spectra of 
which there are, according to Price, also pre- 
dissociations in the same region. In this connec- 
tion we can note also the work of Norrish and 
Noyes,’ who did not find any evidence for the 
formation of oxygen atoms in the far ultraviolet 
photolysis of formaldehyde. 

It can also be noted that as is seen from 
Table I, in the dissociation scheme of CO we 
have the atomic term combination C(?P)+O(?P) 
about in the same energy height as the a *- 
energy level of CO; in this way, at the places 
where the CO(a*Z) state occurs between the 
dissociation products in the given dissociation 
schemes, we have always also a dissociation 
limit corresponding to this atomic term combi- 
nation, this being not represented in the corre- 
sponding tables. 


ACETALDEHYDE 


Reducing the 44.0-kcal./mole heat of forma- 
tion value of acetaldehyde given by Bichowsky 
and Rossini® to 0°K, the data of Table I result 
for acetaldehyde 


CH;CHO-2C +4H(?S) +O(P) 
— 646.2 kcal./mole. 


TABLE IV. Dissociation products and dissociation energies 


of acetone. 
Energy height 
Dissociation products kcal./mole 
1 C.Hs+CO(X!2) 5.5 
3 
4 — 
6 =CH.+CH,+CO(X'2) 84.8 
7 =CH.+CH;CHO 86.3 
8 2—CH;+CO(X!2) 88.2 
9 H(?S) + —CH:COCH; ~98 
10 144.0 
11 C.H.+CO(a*z) 163.8 
12 =C(CH3)2+O0(P) ~179 
13 226.8 


7R. G. W. Norrish and W. A. Noyes, Proc. Roy. Soc. 
A163, 221 (1937). 

8F. R. Bichowsky and F. D. Rossini, The Thermo- 
chemistry of the Chemical Substances (Reinhold Publishing 
Corporation, New York, 1936), 


TABLE V. Dissociation products and dissociation energies 
of ketene. 
Dissociation products ye 
1 =CH2+CO(X!2) 65.6 
2 H(?S)+CHCO ~100 
3 =C=CH.+0(*P) ~184 
4 =CH.+CO(a'll) 204.1 


This value, together with the data of Table I, 


- enables us to calculate the dissociation scheme 


of acetaldehyde given in Table III. 

With respect to the bond energy of the O atom 
in acetaldehyde, process (11) gives ~175.1 
kcal./mole, where we have made the approxi- 
mation of identifying the energy required to 
remove two hydrogen atoms from a methyl 
radical in ethane with the energy required to 
the same process in a free methyl radical. The 
mean bond energy of the methyl radical and the 
hydrogen atom is from process (12) 118.9 kcal./ 
mole, slightly less than the mean bond energy of 
hydrogen in formaldehyde. The energy required 
to remove a hydrogen atom from the methyl 
radical in process (7) was estimated according 
to the corresponding value in Table I. 

The great variety of the dissociation term 
combinations (1)—(8) in the region of photolysis 
is remarkable, all of which can give rise in 
principle to predissociation effects, and have to 
be considered in investigating the variation of 
the absorption and of the acetaldehyde pho- 
tolysis in function of the wave-length of the 
absorbed radiation. Though many possibilities 
of free radical formation can be seen, their role 
in the chain reactions of the photolysis requires 
still further investigation. It is not probable 
that only processes (1) and (2) have to be 
considered in the acetaldehyde photolysis as it 
is often supposed in the literature. Also the 
eventual possibility of free methylene radical 
formation should be pointed out. 

The place of the maximum of the continuous 
absorption observed by Leighton and Blacet* 
corresponds to the —CH;+H(?S)+CO(X'2) 
term combination of process (8). The energy 
height of the dissociation products (2) has been 
estimated in the same way as in formaldehyde, 


®P. A. Leighton and F. E, Blacet, J. Am, Chem, Soc. 
55, 1766 (1933). 
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while the energy of process (3) should be slightly 
higher. 


ACETONE 


Using the 427-kcal./mole heat of combustion 
of liquid acetone, and the values of Table I, 
we obtain at 0°K 


(CH3)2CO-3C('S) +6H(?S) +O(P) 
—926.5 kcal./mole. 


The corresponding dissociation scheme of acetone 
is given in Table IV. 

Process (12) gives for the energy required to 
remove the oxygen atom ~179 kcal./mole, in 
the same way as sketched in the case of acetalde- 
hyde, while the mean bond energy of the CH; 
radicals is according to (13) 113.4 kcal./mole. 

In the region of the observed photolysis we 
have again a great variety of term combinations, 
part of the energetic possibilities being excludable 
by chemical arguments. It is worth while to 
mention a work of Herr and Noyes,!° who point 
out the eventual importance of the process (8) 
of the primary formation of two methyl radicals. 
For processes (3) and (4), the energy heights of 
which should be estimated from process (8), 
we have not given estimated values. 


KETENE 
With a 13.9-kcal./mole heat of formation for 
ketene, there results 


+2H(??S) +O(@P) 

—515.2 kcal./mole. 
Dissociation energies and dissociation products 
of ketene are given in Table V. 


1D. S. Herr and W. A. Noyes, J. Am. Chem. Soc. 62, 
2052 (1940). 


While for the rupture of the methylenic 
hydrogen bond we can estimate ~100 kcal./ 
mole, the bond energy of the oxygen atom is 
about 184 kcal./mole, taking 200 kcal./mole for 
the energy required to remove the two hydrogen 
atoms from ethylene, the approximate character 
of which calculation should, however, be stressed. 
The rupture of the double carbon bond requires, 
according to process (4), 204.1 kcal./mole, which 
energy also contains again the interaction corre- 
sponding to the oxygen atom. 

At the energy height (1) predissociation should 
be observed, and this we can identify with the 
predissociation observed by Norrish, Crone, and 
Saltmarsh" at 3850A ~74.2 kcal./mole. That 
the transition probability of the predissociation 
is smaller in the direct proximity of the term 
combination can be concluded also from the 
observation of Ross and Kistiakowsky,” who 
reported a photo-chemical decomposition of 
ketene at 3660A ~78.1 kcal./mole with a de- 
creasing quantum yield as compared with the 
effect at shorter wave-lengths. According to 
Burton, Davis, Gordon, and Taylor," the ketene 
photolysis at A>3000A gives evidence of the 
formation of free CH, radicals which show high 
reactivity with C.H,. Here, it is again to be 
emphasized that the energy height of a photolysis 
need not be characteristic to the bond energies, 
since it may be due also to predissociation effects 
or to a direct absorption transition into a 
continuum. 


uR. G. W. Norrish, H. G. Crone, and O. Saltmarsh, 
J. Chem. Soc. (1933), 1533. 

2 W. F. Ross and G. B. Kistiakowsky, J. Am. Chem. 
Soc. 56, 1112 (1934). 

13M. Burton, T. W. Davis, A. Gordon, and H. A, 
Taylor, J. Am. Chem. Soc. 63, 1956 (1941). 
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A discussion of different arguments on the dissociation energy of CO and on the role of the 
5S state of carbon is given. The consistency of the 169.735 kcal./mole value of the heat of 
sublimation of carbon into *S carbon atoms given by Geré in spite of some arguments of 
Long and Norrish is shown. The energy required to rupture ethane into two methyl radicals 
results to be {C—C}ethane=82.78 kcal./mole, in excellent agreement with the value which 
can be obtained by a cycle proposed by Van Artsdalen without the use of the heat of subli- 
mation value of carbon. The difference of the mean C—C bond energies in diamond and 
in the higher members of the paraffin series is 2.04 kcal./mole. The energy required to rupture 
the double bond in 1-heptene, {C=C}1-neptene=145.43 kcal./mole, in 2-methyl-1-butene 
{C=C}2-methyl-1-butene = 148.70 kcal./mole. In the acetylene hydrocarbons {C=C}1-pentyne 
= 198.67 kcal. /mole, {C =C Jacetylene = 194.28 kcal. /mole. 


S shown in the preceding paper of Gerd,! 
the dissociation scheme of Schmid and 
Gerd of the CO molecule gives an accurate value 
for the heat of sublimation of carbon, giving 
also the cue to many problems on bond energies 
of organic compounds. Since, however, up to the 
present this value has not been fully accepted, 
and the literature contains many different views, 
further discussion seems necessary. We also try 
here to reply to the arguments of Long and 
Norrish,? and to make some conclusions with 
respect to the bond energies in connection with 
the more recent results of thermochemistry. 


DISSOCIATION ENERGY OF CO 


.The dissociation scheme of the CO molecule’ * § 
given by Schmid and Gerd in 1936 is based on a 
fairly great number of spectroscopical observa- 
tions. There are atomic term combinations at 
the energy heights 11.55 ev, 11.11 ev, 9.61 ev 
and 8.87 ev where predissociations were ob- 
served. At the heights of the atomic term combi- 
nations at 11.11 ev and 10.13 ev term conver- 
gences of the and J' levels were observed, 
at 10.13 ev also a weakening effect at a definite 
rotational level of the A'II state being present.® 


* At present Institut Henri Poincaré, Paris, France. 

1L. Gerd, J. Chem. Phys. 16, 1011 (1948). 

?L. H. Long and R. G. W. Norrish, Proc. Roy. Soc. 
A187, 337 (1946). 

3 R. Schmid and L. Gerd, Zeits. f. Physik 99, 281 (1936). 

4R. Schmid and L. Gerd, Zeits. f. physik. Chemie B36, 
105 (1937). 

5 R. Schmid and L. Gerd, Phys. Zeits. 39, 460 (1938). 

® R. Schmid and L. Gerd, Zeits. f. Physik 106, 205 (1937). 


Besides, the atomic term combinations at 11.55 
ev, 8.87 ev, 8.19 ev and 6.90 ev explain the 
breaking off of emission bands.* Some of these 
effects have no equal independent validity as 
evidence and may be discussed ; the dissociation 
scheme gives account, however, of such a great 
variety of observed effects that it seems hardly 
possible to consider the close agreement of these 
data with the dissociation scheme of Schmid and 
Geré as merely accidental, not even from a 
purely spectroscopical point of view. 

According to this dissociation scheme the 
CO X' ground state dissociates into the C('D) 
+O(1D) atomic term combination with a bond 
energy of 10.13 ev. The general acceptance of 
the dissociation scheme was strongly hindered 
by the view according to which the ground state 
of a molecule should dissociate into the atomic 
ground states. Now it is already known from a 
great variety of molecules that this is not the 
case, and the corresponding view cannot be 
supported either by the results of molecular 
spectroscopy or by theoretical arguments.’ 

Long and Norrish adduce® against the dissocia- 
tion scheme of Schmid and Gerd that the 6.90 ev 
dissociation energy of CO finds no place in the 
empirical relationship of Barrow for the variation 
of bond strength with internuclear distance in 
homologous diatomic molecules. They do not 
take into consideration, however, that in this 


7J. G. Valatin, Proc. Phys. Soc. 58, 695 (1946). 
a tn H. Long and R. G. W. Norrish, Nature 158, 237 
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connection it is not the 6.90 ev energy distance 
of the C(?@P)+O(?P) atomic term combination 
which is to be taken into account, but the 10.13 
ev bond energy of the CO ground state corre- 
sponding to C(!1D)+O(‘D) which fits rather well 
to the more recent data of Barrow.® In view of 
the way of determination of the data for the 
homologous molecules, however, this agreement 
can only be regarded as support for the calcula- 
tions with respect to the bond energies of the 
homologous molecules, and not as a check of 
the bond energy of CO. 

In this connection the close agreement of the 
bond energy of CO with the bond energy of the 
X'Z, ground state of the isoelectronic N»2 mole- 
cule may be remarkable, which latter is according 
to the determination of Schmid and Ger6'® 9.76 
ev corresponding to a dissociation into the 
N(?D)+N(D) atomic term combination. 

Further support of this dissociation scheme of 
CO is given by the interpretation of the photo- 
chemical decomposition of CO observed by 
Faltings, Groth and Harteck." An explanation 
of electronic collision experiments was given 
according to this scheme by Schmid.” Hagstrum 
and Tate adduce™ against this interpretation 
according to their more recent measurements 
that the observed appearance potential of C* 
ions at 22.8+0.2 ev cannot be due to the 
C+(4P)+O(°P) term combination which is above 
this value. Their requirement according to which 
the appearance potential, if in error, should be 
higher than the actual dissociation limit, can be 
satisfied according to a note of Gerd by the 
C+?P)+O(4S) term combination at 22.35 ev. 
The departure is smaller than in O, the spectro- 
scopical dissociation scheme of which is generally 
accepted and where the appearance potential of 
the O* ions at 19.2+0.2 ev is to be interpreted 
by a dissociation limit at 18.7 ev, and is attrib- 
uted by Hagstrum and Tate to a potential curve 
with a maximum of 0.5 ev. The determination 


®I wish to express my sincere thanks here also to Dr. 
Barrow for the kind placing at my disposal of his more 
recent data. See also Proc. Phys. Soc. 58, 538, 717 (1946). 
(1945) F. Schmid and L. Gerd, Csillagaszati Lapok 6, 102 

1 R, F. Schmid and L. Geré, Proc. Phys. Soc. 58, 701 
(1946), 59, 502 (1947). 

2 R. Schmid, Zeits. f. Physik 99, 274 (1936). 
ase iy D. Hagstrum and J. T. Tate, Phys. Rev. 59, 534 


of the limit of error of dissociation energy values 
derived from collision experiments seems to need 
further investigations. 

In connection with the recent note of Hag- 
strum, according to which the dissociation 
energy of CO should be 9.6 ev, the predissociation 
at 9.61 ev corresponding to the C(@P)+O(*P) 
atomic term combination, it can be noted also 
from the purely spectroscopical side that this 
scheme cannot account for the observed conver- 
gence and predissociation effects at 10.13 ev and 
the predissociation at 8.87 ev. The hypotesis 
that the observed predissociation limit at 11.11 
ev, which is also a convergence limit of the a*Z 
levels, should lie 0.25 ev higher than the next 
atomic term combination, owing to a potential 
curve with a maximum, seems also to be very 
improbable on the ground of the observed effects 
at 11.11 ev. 


5S ENERGY LEVEL OF CARBON 


The lowest °S energy level of the tetravalent 
sp® configuration of carbon plays an important 
part in connection with the bond energies of 
organic compounds, this being the energy level 
from which bond energies have to be counted 
probably rather generally. According to some 
theoretical considerations it is rather the energy 
differences counted from the energy level of a 
hypothetical valence state of carbon that are 
characteristic to the bond strengths of organic 
compounds.'®!6 These energy differences are, 
however, by no means observable energy differ- 
ences and differ essentially from the well- 
definable bond energies of the chemical bonds. 
The vibrational levels of the ground state of an 
organic molecule converge to definite atomic and 
radical term combinations, defining by successive 
term convergences of the dissociation products 
atomic terms which have to be taken into 
account in connection with the bond energies of 
the molecular ground state. A given electronic 
state of a molecule can also be characterized by 
a potential function which depends only on the 
coordinates of the atomic nuclei and results by 

4H. D. Hagstrum, Phys. Rev. 71, 376 (1947). See also 
H. D. Hagstrum, Phys. Rev. 72, 947 (1947), in which 
article the observed effects in the CO spectrum are re- 
viewed, and J. G. Valatin, Phys. Rev. 74, 350 (1948). 

15 J. H. Van Vleck, J. Chem. Phys. 2, 20, 297 (1934). 


16 J. Duchesne, P. Goldfinger, and B. Rosen, Nature 
159, 130 (1947). 
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an averaging of the interactions between all the 
particles with respect to the fast electronic 
motions. The asymptotic values of the potential 
function taken at large internuclear distances 
also correspond to definite atomic and radical 
states, and the bond energies are characterized 
by the energy difference of these asymptotic 
energy values and of the energy at the place of 
the minimum of the potential function or of the 
energy of the lowest discrete energy level of the 
molecule. 

The energy of excitation of the °S state of 
carbon has escaped direct observation by atomic 
spectra for a long time, and its exact value, 
33735.2 cm, is only recently known by the 
measurements of Shenstone.!? The fact that the 
~33800 value of the °S excitation energy 
of carbon, deduced from the CO dissociation 
scheme of Schmid and Gerd, represented until 
its direct determination the best one of the 
proposed values,?*!* gives further support to 
the CO dissociation scheme of Schmid and Gerd. 
It should be also stressed, that though the 
approximate energy height of the C(®S)+O(*P) 
atomic term combination was indicated for 
Schmid and Gerd by the older estimated values 
of Edlén and of Bacher and Goudsmit,” after 
the determination of the CO-dissociation scheme 
the slightly different values!® proposed by Edlén 
and by Goudsmit were no more reconcilable 
with the observed predissociation effects in CO 
according to the scheme of Schmid and Gerd. 
We have here, in this way, the rather rare case 
of molecular spectra having yielded an accurate 
value for an atomic term distance which had not 
been until then directly observed in atomic 
spectra. 

As to the remark of Shenstone with respect to 
the lifetime of the 5S state of carbon, the question 
of the metastability of the ®S state does in no 
way affect the role, explained above, of the °S 
state in organic compounds, according to which 
the dissociation limits corresponding to the term 


17 A. G. Shenstone, Phys. Rev. 72, 411 (1947). 

18... H. Long and R. G. W. Norrish, Nature 157, 486 
(1946); J. G. Valatin, Nature 158, 237 (1946); B. Edlén, 
(lot 159, 129 (1947); and S. Goudsmit, Nature 159, 742 


947). 

19 R. Schmid and L. Gerd, Zeits. f. Physik 105, 36 (1937) ; 
B. Edlén, Zeits. f. Physik 84, 746 (1933); R. F. Bacher 
and S. Goudsmit, Phys. Rev. 46, 948 (1934); and see also 
J. G. Valatin, Comptes rendus 226, 702 (1948). 
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convergences and to the asymptotic values of 
the potential functions should correspond to the 
5S energy state of carbon. 

The value of the excitation energy of the °S 
state of carbon seemed to be the most decisive 
point in the discussion on the heat of sublimation 
of carbon with Long and Norrish. After the 
observations of Shenstone, the scheme of Long 
and Norrish based on Herzberg’s CO dissociation 
energy can no longer be held. In the following, 


- however, we have still to consider the arguments 


of Long and Norrish which are of thermochemical 
nature, and to show that they challenge in no 
way the scheme of Schmid and Gerd. 


HEAT OF SUBLIMATION OF CARBON; 
CONSISTENCY OF THE USED DATA WITH 
SOME THERMOCHEMICAL ARGUMENTS 


According to the paper of Geré,' the dissocia- 
tion scheme of CO gives together with recent 
thermochemical data L,=L2=169.735 kcal./ 
mole at 0°K for the heat of sublimation of carbon 
into ®S atoms. The importance of the sublimation 
heat value corresponding to ®°S carbon atoms 
in thermochemical calculations has also been 
stressed by Long and Norrish? who have denoted 
this value with Z2 and pointed out the untena- 
bility of some views which do not take into 
consideration the tetravalency of carbon and the 
role of the 5S state. 

The direct measurements of the heat of subli- 
mation of carbon show highly divergent results 
according to the methods used. With the primary 
process of the carbon sublimation into ®S carbon 
atoms, the results of the dynamic experiments 
which determine the rate of loss in weight of 
heated carbon filaments are in good agreement 
with the value deduced from the CO spectrum. 
The lower value obtained by the equilibrium 
measurements in the carbon arc has been inter- 
preted by Schmid and Geré, who constructed 
also vapor pressure curves which were in good 
agreement with these measurements,?? by an 
equilibrium between ®S carbon atoms and C2 
molecules. In these calculations Schmid and 
Gerd used a dissociation energy of 185 kcal./mole 
for the C. molecule into *S carbon atoms which 
they deduced from an extrapolation of the 

2” R. Schmid and L. Gerd, Mitteilungen der berg- und 


hiittenmannischen Abteilung der Universitat Sopron 9, 
173 (1937). 
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thermochemical C—C, C=C and C=C bond 
energies to a bond with four valencies, a value 
which could also be somewhat larger in the light 
of more recent thermochemical data. It should 
be emphasized that the available spectroscopical 


. evidence on the C; spectrum does not at present 


allow any conclusion as to the dissociation 
scheme of C2, and not even the relative position 
of the singlet and triplet term systems is 
known. The above state of Co, if it exists, may 
be a singlet state, and it may be that all the 
states of the free C, radical correspond to excited 
atomic states as the states of the free CH 
radical.”! 

The great divergence of the results of the 
different carbon heat of sublimation determina- 
tions shows in any case the presence of secondary 
processes, the details of which require further 
investigations. However, with respect to the 
interpretation of recent equilibrium experiments 
as yielding L;~170 kcal./mole corresponding to 
a sublimation into *P atoms, it should be men- 
tioned that this value, leading to the CO dis- 
sociation scheme of Gaydon and Penney,” 
seems to be excluded by other observations. One 
can bring to bear against D(CO) =11.11 ev all 
the predissociation effects observed below this 
value and the convergence of the J'= terms® to 
about 10.1 ev in the CO spectrum, and the 
photochemical decomposition process of CO." 
Besides the spectroscopical and photochemical 
evidence, in the electronic collision experiments 
in CO® there is an observed appearance potential 
of Ot ions at 23.1 ev, whereas the height of the 
lowest term combination C(?P)+Ot(4S) in the 
scheme of Gaydon and Penney is at 24.73 ev. 
Thermochemical arguments are also against this 
scheme since to substitute simply Z; for Ze in 
the present calculations would lead to a faulty 
theoretical picture. For Lz the CO dissociation 
scheme of Gaydon and Penney gives ~267 
kcal. /mole. 

Long and Norrish argue in favor of a higher 
L, heat of sublimation value of carbon also from 


(sat). Gerd and R. Schmid, Zeits. f. Physik 118, 210 
22 A. G. Gaydon and W. G. Penney, Proc. Roy. Soc. 

A183, 374 (1945); A. G. Gaydon, Dissociation Energies 

and Spectra of Diatomic Molecules (Chapman and Hall, 

Ltd. ndon, 1947). 

aed Valatin, }. Chem. Phys. 14, 568 (1946), 15, 336 


the results of Andersen, Kistiakowsky and Van 
Artsdalen*+*> with respect to the energy required 
to remove a hydrogen atom from methane. 
These authors arrive at 100.8 kcal./mole for this 
energy at 0°K and the striking agreement of 
this value with the calculations of Geré! show 
that these measurements can in no way be 
opposed to the Le value given by Geré. The 
calculations of Geré show that the energy re- 
quired to remove the first hydrogen of methane 
is slightly greater than the mean bond energy of 
the hydrogens which results to 98.04 kcal. /mole. 
Long and Norrish suppose the energy required 
to remove the first hydrogen from methane to be 
less than the mean bond energy, and conclude 
from this a larger value of Le. In regard to their 
supposition they refer also to theoretical con- 
siderations of Van Vleck.'® Calculations based 
on the Heitler-London theory, however, cannot 
be expected to give results which would be able 
to decide in questions regarding differences of a 
few percent, or in the improved form given by 
Voge.?® Apart from this, the conclusions of Van 
Vleck do not refer to the former bond energies, 
but to the “gross energy per bond” which is 
counted from an energy level corresponding to 
the hypothetical valence state of the carbon 
atom. 

It is also to be emphasized that the bond 
energies of the CH» and CH radicals occurring in 
organic compounds is also to be counted from a 
term combination corresponding to ®°S carbon 
atoms, as these are given for instance in methane 
in the paper of Geré. The bond energy of a CH 
radical would be accordingly 97.04 kcal./mole?? 
corresponding to a ‘2 state of CH, and this state 
of the CH radical cannot be identified with the 
CH(X°II) doublet state known from band spec- 
tra. In theoretical calculations these two states 
are often identified, and it is supposed that the 
dissociation energy of CH(X°I1) is 80 kcal./mole 
as given by Herzberg. According to the spectro- 
scopical observations of Geré and Schmid,” 
however, neither of the doublet states of CH are 
formed from carbon atoms corresponding to 


**H. C. Andersen and G. B. Kistiakowsky, J. Chem. 
Phys. 11, 6 (1943). 

%G. B. Kistiakowsky and E. R. Van Artsdalen, J. 
Chem. Phys. 12, 469 (1944). 

26H. H. V. Voge, J. Chem. Phys. 4, 581 (1936). 

27 The slight deviation from the value in Geré's paper is 
explained later in the text. 
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energy terms of the lowest s?p? configuration. In 
this way, it is also doubtful whether CH, 
possesses energy levels lower than the levels of 
the free methylene radical. 

Long and Norrish adduce ‘three lines of 
thermochemical evidence” in favor of their 
energy of excitation value of the °S state and of 
their value for Ze. All the three arguments start, 
however, from suppositions which do not seem 
justifiable. 

(a) The conclusions concerning the molecular 
structure of CH» are based on Herzberg’s men- 
tioned CH dissociation energy value, and refer 
to an investigation of Burton, Davis, Gordon 
and Taylor?* on the photolysis of ketene. The 
ketene photolysis can be explained satisfactorily 
by the production of methylene radicals which 
explanation seems to be preferred also by these 
authors, and does not involve the attempted 
consequences. 

(b) The stability of the CHO radical fits also 
into the scheme used throughout this paper.” 
Besides, owing to the interaction of the H atoms 
with the O atom, the smaller force constant of 
formaldehyde as compared with methane does 
not involve a smaller ‘dissociation energy” of 
the hydrogen atoms or the missing of an ‘‘energy 
of reorganization’”’ of the CHO radical. The 
possibility of predissociation effects is also to be 
emphasized.” 

(c) Long and Norrish thought to be able to 
conclude to a lower value of the C(5S) excitation 
energy from the fact that according to the 
observed force constants, bond lengths and ion- 
ization potentials the bond strength of COz is 
smaller than that of CO. But according to the 
scheme of Schmid and Geré we have also CO.— 
C(®S) +20(?P) — 380.86 kcal./mole, involving a 
mean bond strength of 190.43 kcal./mole, com- 
pared with the 233.39 kcal./mole bond energy 
of CO(X!Z) which corresponds to the dissociation 
into C(1D)+O('D) atoms. The calculated 162.7 
kcal./mole bond energy of the oxygen atom in 
formaldehyde” is again a very reasonable value, 
and does not make necessary any such conclu- 
sions. 


28M. Burton, T. W. Davis, A. Gordon, and H. A. 
Taylor, J. Am. Chem. Soc. 63, 1956 (1941). 
(1948) 6 and J. G. Valatin, J. Chem. Phys. 16, 1014 


BOND ENERGIES 


The most important quantity which can be 
given in order to characterize the strength of a 
chemical bond is the energy required to dissociate 
the molecule into definite excited or unexcited 


atoms and radicals which correspond to the 


asymptotic energy values of the potential func- 
tion of the given molecular state. This bond 
energy is, for instance, for the dissociation of 
methane into a H(2S) atom and a free CH; 
radical 100.8 kcal./mole. This is the energy 
required to rupture the corresponding bond, and 
the name bond energy is more suitable for the 
reason that dissociation energy means in spectro- 
scopy the energy height of the corresponding 
lowest term combination without respect to the 
possible role of excited atoms or radicals in the 
molecular ground state. On the other hand, the 
mean bond energy of the four hydrogen atoms 
in methane is only about 98 kcal./mole. Owing 
to the interaction of several atoms, this quantity 
cannot always be defined with respect to a 
definite bond, and for instance, the mean bond 
energy given for the C—O bond in CO, includes 
also an interaction energy of the O atoms and 
differs considerably either from the average 
interaction energy between an oxygen atom and 
the carbon atom, or from the energy required to 
remove an oxygen atom. 

Besides, from a theoretical point of view the 
energy values characterizing directly the po- 
tential function are of greater interest, leaving 
out of consideration the part of the dynamic 
effects which are less characteristic of the inter- 
action of the constituent atoms. From this point 
of view it would also be interesting to know 
what part of the slight differences in the bond 
energies of isomeric molecules is due to bonding 
forces, and what to the differences in the vibra- 
tional zero-point energies. 

In order to obtain values characterizing a 
single molecule, besides spectroscopical data we 
use in the following heat of formation values 
calculated to 0°K. 


PARAFFIN HYDROCARBONS 


The evaluation of the results of recent thermo- 
chemical measurements allows by the knowledge 
of the heat of sublimation value of carbon to 
make some further conclusions with respect to 
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the bond energies of hydrocarbons. In order to 
obtain heat of formation differences character- 
izing a definite bond, we propose to subtract 
from a given heat of formation value the corre- 
sponding values for the hydrocarbon which 
results by replacing the given radicals with 
hydrogen atoms. 

For instance, the heat of formation differences 
of 3-methyl-heptane and n-heptane, 3-methyl- 
hexane and n-hexane, and 2,4-dimethylhexane 
and 2-methylhexane give according to the data 
of Prosen, Pitzer and Rossini®® equally within a 
hundredth 3.98 kcal./mole, and this figure is 
characteristic of the bonds of a methyl radical 
on the third carbon atom of a paraffin chain. 
With the same hypothesis as in Ger6’s paper 
that the strength of this methyl bond equals the 
mean C—C bond resulting from the 3.673 
kcal./mole heat of formation difference of suc- 
cessive v-paraffins, the 3.98—3.67=0.31 kcal./ 
mole difference is the excess energy of the 
methyl radical as compared with the bonds in 
methylene, 


ACH;={C—H}cu;—[C—H 
=0.31 kcal./mole, 


where {C—H}cu; means the energy required to 
remove the first hydrogen of the methyl radical. 
The hypothetical character of this assumption 
is, however, to be emphasized, leaving a slight 
uncertainty in the value of the given bond 
strengths. 

In the same way, the difference of heat of 
formation of 2-methylpentane and n-pentane, 
2-methylhexane and n-hexane, 2-methylheptane 
and n-heptane, 2,4 dimethylhexane and 3-methyl- 
hexane give almost equally 4.78 kcal./mole, and 
the difference 4.78—3.67=1.11 kcal./mole can 
be interpreted as the excess energy because of 
the second methyl radical at the end of the 
paraffin chain, the average excess energy per 
methyl radical being in the case of two methyl 
radicals 1.11+0.31/2=0.71 kcal./mole. 

We have in the -paraffins beginning from 
hexane, 


(n—1)[C—C ]n-parafin + 2(n+1)[C —H 
(1) 


%E. J. Prosen, K. S. Pitzer, and F. D. Rossini, Bur. 
Stand. J. Research 34, 403 (1945). 


[C—C]n-paraftin + — H Jon, = L.+D(H2) 
+3.673 =276.63 kcal./mole, (2) 


and the two equations give 


4[C H Jcu,+2ACH; =L.+2D(He2) 
—(AHf°+(n—1)-3.673). (3) 


The value of the expression in the parenthesis is 
for the higher m-parafins beginning from hexane 
equally —12.61 kcal./mole, and this gives with 
L.=169.735 kcal./mole, D(H2)=103.22 kcal./ 
mole, ACH;=0.31 kcal./mole for [C—H ]cne, 
the mean bond energy of CH» which is supposed 
to be identical with the energy required to 
remove a hydrogen atom from a free methylene 
radical 
[C—H Jcu, = 97.04 kcal./mole. 


The mean paraffinic carbon-carbon bond strength 
which contains afso the slight interaction of the 
second next methylene radicals results to be 


[C—C]),-parafin = 82.55 kcal. /mole. 


For n=2 the analogous equation to equation 
(1) gives directly the energy required to rupture 
the bond of the two methyl radicals in ethane. 
With AHf’= — 16.517 kcal./mole, 


{C—C} ethane = 82.78 kcal. /mole. 


The cycle proposed by Van Artsdalen* gives with 
the 100.8+1 kcal./mole value of Kistiakowsky 
and Van Artsdalen* for the energy of rupture of 
the first hydrogen atom in methane at 0°K 


2CH, =2CH;+2H—201.6 kcal./mole 
2CH, + 15.46 kcal. /mole 
2H = H, +103.22 kcal./mole 
= 2CH; — 82.9 kcal./mole. 


The excellent agreement of the two values for 
{C—C} thane Shows again the reliability of the 
used value of the heat of sublimation of carbon. 

It is worth while to note also the relation 
between the mean [C—C],, bond energy 
in the higher x-paraffins and the mean C—C 
bond energy in diamond. This can be given by 
means of the recent thermochemical data inde- 
pendently of the value of the heat of sublimation 
of carbon. Equation (3) gives 


+6.31 kcal./mole, 


3tE. R. Van Artsdalen, J. Chem. Phys. 10, 653 (1942). 


fa 
ate 
the: 
nc- 
ond 
of 
rgy 
and 
the 
tro- 
ling 
the | 
the 
the 
oms 
ying 
sty 
oa 
ond 
ides 
rage 
d to 
po- 
ving 
yMic 
iter- 
point 
now 
yond | 
ding 
ibra- 
ig a 
a we 
ilues 
edge 


1024 J. 


and subtracting this equation from Eq. (2) 
results 


[C—C]n-parafin = L-/2+ACH;— 2.64 kcal. /mole. 


Since according to Wagman, Kilpatrick, Taylor, 
Pitzer and Rossini® 


Legraphite =0.577 kcal. / mole, 


[C—C ]siamona = Laiamona/2 
=L,./2—0.29 kcal./mole, 


[C C Jparaftin og +ACH; 
— 2.35 kcal./mole. 


With ACH;=0.31 kcal./mole this gives a differ- 
ence of 2.04 kcal. /mole between the mean carbon- 
carbon bond energies of the higher members of 
the paraffin series and of diamond. 


MONOOLEFIN HYDROCARBONS 


By means of the heat of formation values of 
Kilpatrick, Prosen, Pitzer and Rossini,** we can 
determine the energy values required to rupture 
the olefinic double bonds of monoolefins. With 
the same method as in the case of the paraffins 
we form the difference 


AHf° 1-heptene — -nexane = — 6.37 + 30.98 
= 24.61 kcal./mole. 


With the approximation that the energy differ- 
ence of m-hexane and the radical formed from 
1-heptene by removing the = CH, radical equals 
2(C—H ]Jcu2+ACHs, we get for the energy re- 
quired to rupture the double bond in 1-heptene 


{C =C } 1-heptene = L.+ACH;-— 24.61 
= 145.43 kcal./mole. 


By an analogous assumption the energy required 
to rupture the double bond in 2-methyl-1-butene, 
in which the remaining paraffin radical also 
contains two methyl radicals, results using 


+23.332 =21.03 kcal./mole, 


{ C=C } 2-methy]-1-butene =L,.—21.03 
= 148.70 kcal./mole, 


® Wagman, Kilpatrick, Taylor, Pitzer, and Rossini, 
Bur. Stand. J. Research 34, 143 (1945). 

% Kilpatrick, Prosen, Pitzer, and Rossini, Bur. Stand J. 
Research 36, 559 (1946). 


G. VALATIN 


which shows a stabilization owing to the inter- 
action of the double bond with the neighboring 
methyl radical. 


ACETYLENE HYDROCARBONS 


Using the data of Wagman, Kilpatrick, Pitzer 
and Rossini,** we obtain by the same method the 
energy required to separate the =CH radical in 


1-pentyne. 


pentyne — butane = 38.90 + 23.33 
= 62.23 kcal./mole 


{ C=C } 1-pentyne =L,.+2[C = H — D(H2) 
— 62.23 = 198.67 kcal. /mole. 


In acetylene 


{ Cc =C } acetylene H 
2L.+D(Hz2) —AHfracetyiene 


gives with AHf°acetyiene = 54.33 kcal. /mole 
{C=C} acetylene = 194.28 kcal./mole. 


As it is known from bond lengths and force 
constants, the C—H bond in acetylene is appre- 
ciably stronger than in the paraffins, this example 
showing clearly again that the sum of the 
energies required to rupture the different bonds 
may differ essentially from the heat of formation 
of the molecule corresponding to a dissociation 
into the constituent atoms. 
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A Tabulation of General Formulas for Inverse Kinetic Energy Matrix Elements 
in Acyclic Molecules 
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In terms of internal coordinates of the bond stretching and bending type, together with a 
torsional type of coordinate, general formulas for inverse kinetic energy matrix elements 
have been computed for all 33 possible acyclic configurations involving the given types of 
coordinates. Tabulations are given both for the general case and for three special cases in which 
all valence angles of the configuration are assumed to be 90°, 109°28’, and 120°, respectively. 


I. INTRODUCTION 


ENERAL methods for setting up the secular 
determinant for molecular vibrations have 
been given by Wilson!? and by Eliashevich? in a 
form which makes use of internal coordinates 
adapted to the valence structure of the mole- 
cule. Since the transformation from this basis to 
a basis adapted to the irreducible representations 
of the group of symmetry operations is now well 
known, and since methods for performing the 
numerical solutions of the separate factors of the 
secular determinant are being facilitated through 
the use of machines,*~’ the calculation of the 
elements of the inverse kinetic energy matrix 
seems likely to emerge as one of the more 
laborious steps in the process of computing the 
vibration frequencies of a molecule. In this con- 
nection, it might be mentioned that a recent 
paper by Wilson® provides a scheme whereby the 
potential energy coefficients can perhaps most 
conveniently be determined as unknowns to fit 
a given set of observed frequencies; this method 
similarly makes use of the inverse kinetic energy 
matrix. 

In view of these facts, a tabulation of general 
and special algebraic formulas for the inverse 
kinetic energy matrix elements has been prepared 

*This work was supported by the Office of Naval 
Research under Contract N6ori-88, Task Order No. 1. 

LE. B. Wilson, Jr., J. Chem. Phys. 7, 1047 (1939). 

2 E. B. Wilson, Jr., J. Chem. Phys. 9, 76 (1941). 

1940) Eliashevich, Comptes rendus U.R.S.S. 28, 605 
1G K, Carter and G. Kron, J. Chem. Phys. 14, 32 
aA ‘A. Frost and M. Tamres, J. Chem. Phys. 15, 383 

6 RH. Hughes and E. B. Wilson, Jr., Rev. Sci. Inst. 

18, 103 (1947). 


W. A. Adcock, Rev. Sci. Inst. 19, 181 (1948). 
SE. B. Wilson, Jr., J. Chem. Phys. 15, 736 (1947). 


in the hope of reducing the time consumed in 
calculating normal vibration frequencies. These 
formulas, which express the result of the trans- 
formation from Cartesian to valence-type internal 
coordinates of small displacements of the atoms 
from their equilibrium positions, may all be 
derived from the fundamental formula as given 
by Wilson,? 


N 


k=1 


in which g,, is a matrix element associated with 
the coordinates g and q’; Sq, is a vector repre- 
senting the contribution of the kth atom to the 
coordinate q, u; is the reciprocal mass of the kth 
atom, and the summation is extended (formally) 
over all N atoms of the molecule. 

Two types of internal coordinates have been 
principally used to date: the change in distance 
between two bonded atoms and the change in 
angle between two bonds having an atom in 
common. A third type of coordinate is suggested 
by the resistance, to torsion about, e.g., the 
carbon-carbon bond in ethane. As a general 
definition of a coordinate of a third type, we 
adopt the change in dihedral angle between 
planes of atoms (123) and (234), where neither 
(123) nor (234) are colinear and where bonds 
exist between 1 and 2, 2 and 3, 3 and 4. Sucha 
coordinate appears naturally in the C2. model of 
H,O2, for example. 


Il. NOTATION AND CLASSIFICATION OF THE 
INVERSE KINETIC ENERGY MATRIX 
ELEMENTS 


The three types of coordinates described above 
will be identified by the symbols 7, ¢, and + 
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INTERACTION OF 

STRETCH AND HgC-H, 

BENDING IN METHYL 
GROUP 


DIAGONAL TERM FOR 
ANY BENDING 


Fic. 1. Examples of notation. 


referring to stretching, bending, and twisting, 
respectively; thus g,, will represent a matrix 
element involving the interaction of a stretching 
and a torsional coordinate. Similarly, S,, will 
represent the contribution of the kth atom to a 
bending coordinate. The equilibrium configura- 
tion of the molecule will define the following 
set of symbols numerically: 


Qi;=Qj:=reciprocal of equilibrium distance 
between ith and jth atoms, 
disk = =%;) =equilibrium angle between 
bonds i—j and k—j, 
Tijkl= T ikjil Ti = equilibrium dihedral 
angle between planes of atoms (i;.) and 
(ikl): 


TaBLE I. Notations and numbers of distinct configura- 
tions for inverse kinetic energy matrix elements in acyclic 
structure. 


(0) 

(1) 

(2), (i) 

(2), (2) 

G), 

(6) 

(1), 

2)» (i)s (a) 

(2), (1), () 
G3), @), G), () 
(9) 

(i), () 

G), @), G) 
(3) (3) 


As will be seen below, the simplified notation in 
parentheses will be used, or the subscripts may 
be omitted entirely, wherever no ambiguity 
arises. The angles ¢ and 7 will have the following 
ranges: 0<@<a2; —t<7r<+m. The sign of 7 is 
determined through the following convention: 
viewing the configuration (ijkl) along the line 
j—k, with j nearer the observer than k, let 7 be 
positive if it is traced from the projections of 
i—j to 1—k in the clockwise sense. 

Wilson? has given expressions for and 
(the contribution of the kth atom to the stretch- 
ing and bending coordinates, respectively) in 
terms of the above quantities and of unit vectors, 
e:.1, directed from the kth to the /th atom: to these 
formulas, we now append appropriate ex- 
pressions for the torsional type of coordinate 
(S,x). It is assumed that the atoms are bonded 
together in numerical order. 


21, (2) 
{12}Si1, (3) 
(€12 X 


(cosp2€21 — €23) = Q21€12 X 
singe 


So1= 


singe 

— (Soit+Sys), 
So3= {13}So1, 

C12 X€23 


S.= ’ 
sings sings 


S., -| Qie 


sings 


€12 X C23 
—(costi4 ’ 


sings 


C12 X C23 
—cotds sint14Q23€23 X ( ). (8) 
singe 


(14) (23)} (9) 


and 
= { (14) (23) } 


In Eqs. (3), (6), (9), and (10) the expression 
in brackets signifies a permutation of atomic 
subscripts to be performed upon the previous 
expression as indicated. 

The classification of the matrix elements 
depends, in the first place, upon the fact that 
they can be non-vanishing only if the two coor- 


(10) 
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TABLE II. General formulas for inverse kinetic energy matrix elements in acyclic structures. 


Formula 


rr 

» 


(i) 
( 


(3) 
(5) 
(3) 


(3) 
(i) 
(i) 


(3) 


() 
(3) 


—singQesp2 
sing; CosTQi 3411 
— 
singsi4 
sings cotd; sin 
0 
sing 


sintosQ 213 
sing 


ind! 

Q? + 2343+ 2 CosQ12023) we 


+ Ke — €0S$123023 — 24) 
Sing23 


COS324 COSH123 COSH124 — COS"H123 — COS*h124+ 1 
314 


—cosrQ12[ 41 + (Qi2—cosd2Q23) ue 
— (sin725 sint34+COsT25 COST34 COSP1) 
[(singo14 COSH415 COST34— SiN215 COST35) COSH415 COST35 — SiNG214 COST34)Q15 

415 
C (4435 — (W4453 — 012015 + (Y4235 — 44435) 
, 
+ (W’253 4453) 
sint 
sings 
singse4 Sind*;3 cotds sinris 


12@1+ | sind®13(Q12 — Cosd123Q23) [sings24 Cosd°13(Q12 
Singi23 SIN124 singi23 


42 — —COSH3Q23) us ] 


+sing124(Q23 —Cosd123012) us 


sing: 


ind 
+ sind°34+cot¢: sinr35)Qi2— (Qi2—COsd123Q23) ue 


singi24 


1 


sing, 


[ (singars COSa216 SINT 34 —COS214 SiNGa16 (Sing214 COSH216COST 34 — COSH214 SING216 COST36) COSH214 


(sinr25 COST34 —COST25 SiNT34) 


—sin*po14 singaie sints4 sind® | sintseQis] Qis Mi 
singais 


C(singsie sinr4s— singers Sint24 COSH216) (Singais SINT 24— SINT4¢ COSH216) Qis Ge 
sings singais 
[Lsingars sind*s2(Q1s —Cos21sQi2) +singsis 
+cotds COSd*s2(Q1s +singsis a 


singers 
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TABLE II.—Continued. 


Type 


Formula 


ing ing 


end) 


sing, 


Ore 
ing 


7] 
singe 


V5 O12 


sing123 singi25 singies singie5 


(3) 


2 
+) Mit SiINT14 SiN T15Q23) we + +Cot? sinr14 SiNT15Q"23) 


+cotde sin715541) Jus 
COST14 7] Qs2 


wat [ cost14)Q32— 


sings Ising: 


Cots COsd2 sinr14 sinr35Q12Q32 


+cotdi25 Cotds sinri¢ ]+cotds sin714 Sin71¢ 


COST36 7 Qi2 


315 


(i) [ (cote, COST36-+COtd125 


sing; 


— (COS5® + COSP123 SINT14 


—cotd: Cotds COSP123 SINT14 SINT56 wo 


(3) cots sinr35 sint26Q15) 


sing 


cosd2 sint35 sin714Q23) = ] 
2 sing 


Cotdies SiNT35 SINT 46712) SINT35S 264 12 
+ CotHars Sin735 SiNT46Q712) SiNT465""35) O12 Jue 


35 Sint 25 COS@1 — COST35 COST 26) Ou. 
singe Sings 
Quiz 


sings17 


(3) [(cosrs sind®72 Cosd215) 


(2) COS617 — COS217 COSH416 
2 


cotd, sints7 
singai7 


cotde 
singers SiN417 singers 


— [cotds17(cosr35 sind®72 Cosda15) 


+cotds(cosrss5 Cost26+sin735 SiINT25 COSP215) 106} 
sing: 


Sind°¢7 — COSda16 SiNda14 sin S"75Q12 


cotd, sinrss sinra] 


0Sp214(COSH214 COSH617 — COSH416 COSH217) sind vin 


singare SiN417 


dinates involved possess one or more atoms in 
common. In the subsequent tabulations this 
number of common atoms is indicated by the 
superscript p; writing the subscripts always in 
such an order that a<a’, where a=number of 
atoms in first coordinate, a’=number of atoms 
in second coordinate, it is clear that 1<p<a. In 
the present paper, only acyclic configurations of 
atoms defining the superposition of coordinates 
q and q’ will be considered, i.e., the collection of 
atoms and bonds will be ‘‘trees” in the topological 
sense as defined by Cayley.® This being the case, 


®A. Cayley, Phil. Mag. 13, 172 (1857). 


it is possible to specify the distinct types of con- 
figurations uniquely through the following system 
of notation : imagine the set of common atoms to 
be drawn in a horizontal line, with the residual 
atoms of coordinate g drawn in lines with slopes 
of +45° upwards away from each end of the 
common set, and with the residual atoms of 
coordinate q’ similarly drawn downwards (see 
Fig. 1). Then a two by two matrix, C, whose 
elements are the numbers of residual atoms in 
the corresponding corners of the pattern just 
described, would define the configuration pre- 
cisely. However, only the first column of this 


§ 
g 


g 


1 
- 
g 
g 
g 
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matrix need be given, since Cyi+C,2.=a—p and 
Con t+Cy2=a'—p. Furthermore, in computing 
the number of distinct configurations, the 
equivalence of the atoms in each coordinate by 
pairs beginning at the ends and progressing 
towards the center must be considered (corre- 
sponding to the permutation symbols in Eqs. (2) 
to (10)). For example, if g is a stretching coor- 
dinate (type 7) and q’ is a bending coordinate 
(type ¢), and an interaction term with two 
common atoms is under consideration (g?,,), one 
could represent such an interaction by the 
matrices or which correspond, respec- 
tively, to cases in which the non-common atom 
of the ¢ coordinate is the first or third atom in 
the set defining ¢, reading from left to right in 
the standard configuration diagram described 
above. Because of the equivalence of the S 
vectors for the two atoms of coordinate r and 
for the first and third atoms of coordinate ¢, the 
two matrices represent equivalent interaction 
terms; in the tables, the first column of C is 
given where C is chosen in such a way that 
Cy >Cy, and if such that Cy1>Cyr. 
Explicit solutions of the combinatorial problem, 
giving the number, n”,, of distinct g? qq’, in 
terms of a, a’, and p may be readily obtained, 
but the results are merely summarized in Table I. 
The columns headed by gq, q’, p, and ?qq need 
no further definition; the final column gives the 


first column of the matrix C for each distinct 
2 qq (i.e., the number of residual atoms of coor- 
dinate g above a similar number for coordinate gq’ 
at the left-hand end of the pattern). 

Finally, it is necessary to establish a conven- 
tion for numbering the atoms in each of the 33 
possible configurations described above. The 
system followed throughout the tabulation of 
formulas given below in Section III is as follows: 
the numbering starts with the atoms common to 
both coordinates, reading from right to left along 
the horizontal line described above; the residual 
atoms are then numbered outwards from the set 
of common atoms, taking the upper left, lower 
left, upper right, and lower right corners of the 


_ diagram in that order. It is furthermore to be 


understood that C,,; and C, refer to the total 
number of residual atoms of each coordinate at 
the left-hand end of the diagram. This is illus- 
trated in Fig. 1 for several configurations; ex- 
amples of coordinates from well-known molecules 
are cited with each configuration diagram in the 
figure. 


III. TABULATION OF FORMULAS 


In the succeeding tables, it will be found that 
the algebraic expressions become rather cumber- 
some toward the end in the general case. While 
one purpose of the collection of formulas is to be 
as explicit as possible, some abbreviations are 


TABLE III. Special formulas: all ¢=90°. 


grr 
0 
— 


cost 
(i) 
0 
(3) 
1241 + + (0712+ ue 
() 
£1496(3) —sintes 
(3) 
Sor 


+0" 12(u1 + m2) [+613] 
(1) = +Qoqu2) [2654] 
(3) 
costes 
0 
(3) 
(2) 
Birr Q?12(ui + m2) +Q?s4(ust+ ps) 
(0) — +0087 1403403243 + 
0 
(7) —cosrs6Q12Q1 61 
(6) cosrasQ*12(u1 +2) 
(i) 0 
—COST35 cost26Q120 
@) 0 


| 
l- 
n 
of 
se 
in 
st 
is 
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TABLE IV. Special formulas: all ¢ = 109°28’. 


(24/24) 1241 + [ O28) [3012+ (1 +008715) O24] 
Brg — [613] 
cosrQrsy/3 (i) — (24/24) [ [18 sinrss)Q12 

(i) /3 + (3Q12+Qe3) (43'+42 sinrss)Qi2 Jue] [+834] 
sintQe3u2/3 (i) (2-4/3) { [9 sintssQi2+ COST34 

(3) (2 sint24+34)Q13 (3) (24/36) Q12[ [3Q16+ (1 +cosrs6)Q12] 

[+605 ] + +4- sinr34) sints6Q14] M1 [2846] 

+ + (30712 + 2012023) 
+[3Q%21+ O24) Q21 — 5Q23Q24 Jue] 

—(1/3) 41+ (3Q12+Q2s) ue] 
—(1/3)(3 sinres sinrs4—Cosres 

(i) 

(3) (24/24) [ + (F3!—2 sinrss) (3Q12+Q15) 1013] [+6425] 


Str (1/8) [9Q*1241 +90? [90712 +2(1+cosr) 723 Jus 


Jus] 


(1/16) {1801241 + [18012 (1+2 cosris+2 cosris)Q%23 
+[(1+2 cosris4+2 (2 cost14— 1)Q34+(2 costis— us} 
(3) — (1/24) {9Q12[ (1 +cos735)Qi2+3 cosrssQ15 Jui +9Qe23L (1 cosrisQs4 Jus 
+ [3(3Q12+ (1 +c0s714)Q23 |[3Q2s+ (1 +cos735)Qi2]— sint14 sint35Q12Q23] } 
— (1/96) {540% 121 + 1540212 +3012[(6(1 +c08714) +34 sinz16)Q25] 


+ [6(1 +cosr14) (1 +cosri¢) —10 sinr14 714) +sinrig—sin714) ]Q23Q25] ue} 


[+6415] 


(3) (1/48) { [9(2 COST 63 — 1)Q12+54 coste3Qi6 = [303Qi2+ (1 +cosri4)Q23 [3Q25+ (1 +cos756)Q12] 


sinr14[3Q25+ (1 +cos756)Q12 1023 F3! sints6l3Q12+ (1 +c0s714) 023 ]Q12—sin714 Sint56Q12Q23] we} 


[+6553] 


(5) (1/8) Qi2{ [9 cosrssQi2+[3 cosrss(1-+cosr26) —sinzss sinr25 


+ [9 cost35Q12+ [3 cosr35(1 +cosr14) sin735 sint14 ]Q23] pe} 


(i) (1/16) { [ [3Qis+ (1 +cosr753)Q12 J[3Q16+ (1 +cosr64)Q12 ]+sin7ss sint6sQ"12 


+3 (1+cosr53)Q12 34 sinzssQ12[ 3Q16+ (1 [+656] 
+[[3Q23+ (1+ cos735)Qe1 J[3Q21+ ]+sinzss 


(3) — (3/8) (3 cosrs5 cosr26+sin735 sint26)Qi12Q 


(3) (1/16) +sint35— 33 cosT35) Qirt+ [3 cosr35(2 COST 26 — 1) + sint35( +3!—2 sin T26) 1015] 


] 


(3) (1/24) { (1+cos775)Q14 (1+cosr6s)Q12 1014] 


—2 sintes sint7sQi2Qi4} ur [6°67] 


imperative. The following have been adopted: 
=W uj = — COSHE 
= (Qi;/singijx) — (Cost ix Cotd +-Cotdijx) 
and 
= Tik — Tjk- 
The transformation properties of the several 
types of g”,q° should be noted carefully. Under 


a reflection (generally an improper rotation), 
coordinates of the first two types are unchanged, 
whereas the third type, the twisting, suffers a 
change of sign. 

Thus, the blocks of the G Matrix corresponding 
to elements of the types g”,, and g4, for enantio- 
morphous configurations will be found to change 
sign. The eigenvalues are unaffected, of course, 


AS 
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since similar blocks in the F (potential energy) 
matrix also change sign. 


Table II contains the general expressions (the 


symbol ( is omitted wherever mw =t); 


in cases of extreme symmetry, it may occasionally 
be found preferable to develop the formulas by 
simplifying Eqs. (2) to (10) first and substituting 
in Eq. (1), using the table as a check. In cases 
practically devoid of symmetry, the general 
formulas serve both for direct calculation, and as 
a means of assessing ahead of time the magnitude 
of the numerical labor. 

In addition to the general case, the appropriate 
expressions for the g,q° are given in Tables II], 
IV, and V for the cases in which all ¢ assume the 
special values of 90°, 109°28’ (tetrahedral angle), 
and 120°. In Tables III and IV, some of the 


formulas will be found to contain a choice of 
sign. The proper sign is readily determined by 
taking the upper sign throughout if a 6*;; 
enclosed in brackets after the formula in question 
is positive or the lower sign if 6*;; is negative (the 
sign of 6*;; follows the same convention as that 
for 7;;). There is just one case, namely, g®,,(}) in 
Table IV, in which there occur two simultaneous 
choices of sign. As indicated in the table, upper 
signs are to be chosen for the coefficient of 1, if 
6°56 is positive, lower if 6°55 is negative, whereas 
6°34 provides a similar criterion for the choice of 
signs in the coefficients of po. 

Finally, it should be noted that the formulas 
for gigs(i), and g',,(3) are left blank in 
Tables III and V. This is because the associated 
configurations involve a tetravalent atom at 
which it is geometrically impossible to satisfy the 
condition that all bond angles be 90° or 120°. 


TABLE V. Special formulas: all ¢ = 120°. 


—p/2 () 
(34/2) costQismi 
(38/2) 
@re(2) (1/2) sinrQospe Sor 
(3) sinr2sQ (i) 
(1/2) (0) 


+ +0723 + 012023) g'or(3) 


(1/2) [2071241 + [0 12(2012 + Q23 + O24) — 24 luo] 
— (1/2) + (2012+ Q23) 
(1/2) (costs costs4—2 sinres 


(1/2) 


— (34/6) us] 

(34/6) sinr15(Qi2— Q23) Qoame 

— (34/6) + (2Q12+ Qos) 

(34/6)Q12{ [2 sinrss[2Qi12+ costssQ14] 


+2 sinrs4(2Q12+Qes) uo} 


— (3-4) (sinr25 costss+2 Costes sints4)Q12Q 


(34/6) Qr2[4 sint36Qie+ [2(1 +cosrse6) sint56 ]Q14] 


(3) (374) sinras(Qire— 


(2/3) [2Q*1241 + [2012+ (1+ cosr) 0723+ 2(1+cosr)Q12023 Jus 


Jus] 


(2/3) {2Q*1 241 + 2Q12(Q12 +Q23) u2— [2QssQ3s+ ime +cosri4)Q35+ (1 —cost14)Qs4 us} 
() — (1/6) {4[(1+cos735)Q12+2 cosrssQis costisQs4 
+ [4[(1 +cos735)Q*12+ (1+ cosr14) ]+ [8+ 2(1 —sinz14 sint35 wo} 
— (1/6) [2[2Q12+ (1 +c0s714)Q23 ][2Q12+ (i +cosr16)Q25] —sinri4 sint16Q23Q25] ue} 
(i) (1/6) {4[(1 —cost6)Q12—2 cosrseQ16 [ 2[2Q12+ (1 +c08714) 023 [2025+ (1 +cos7s6)Qi2] 


—sinri4 sin756Q12Q23] Me } 


(3) (1/3)Q12 { [ 2 +cosr26)Q15] sinr35 sint26Q15] Mi 


+[2 cosrss[2Q12+ sinr14Q23] 


G) — (2/3) [[2Q1sQ16+ Jui + [20230214 uo] 


(2/ 3) (sinr35 sint26+2 cosr35 cost 26)Q12Q15m1 


(3) (1/3) [4 cost3sQi7+ [2 —cosro¢) +sin735 sint25 1015] 
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IV. EXAMPLES OF APPLICATION TO — 
SPECIFIC MOLECULES 


1. (Cy Symmetry) 

Choose as internal coordinates the O—O and 
O-—H stretches (calling the coordinates P and Ri, 
respectively), the O—O—H bending (a;), the 
H—O-—O-—H twisting (8); the index 7 on R and 
a identifies the H atom involved. Let Qo-o = Po, 
Qo-n=Ro, ¢0-0-n =, TH-o-o-H Then the 
following elements of the G matrix associated 
with the indicated coordinates may be found 
from Table II. 


TABLE VI. 
Coordi- 
nate 
pair Type Formula 
PP 
PR; opto 
Pai —sinaoR oe 
Pp 2 0 
RiRi rr 
RiR2 0 
reo —sinaoP oo 
Ries sina COSBoP ope 
Rip —cosao sinBoP ono 
ajar; (2P%+ R%o—2 cosaoP oRo) po 
a1a2 —2 cosBo(Po—cosaoRo) Popo 
sinao 
BB [ +[R%0+2 cosao(1+cosBo) 
X (cosaoPo— Ro) Po 
Py Ru Rie Ra Roe Si Se 
Ru gre 0 grs 
Si gs 0 
Ti 
Bi 
621 522 €2 
P; g'Ps g'Ps Ps g'Pe Pc 
Ru g'Rs 0 0 0 
Si £88 £88 0 0 £Se 0 
0 0 fae 0 
Bi £B8 £85 0 0 0 
bu g's 25 gs 
€1 gle 
ou 
n 
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2. (Propane; Symmetry) 


Assume tetrahedral valence angles throughout; 
consider one H on each methyl group and the 
three C to be coplanar and of the form indicated 
in Fig. 2. Then all 7 are either +60° or 180°. 
With atoms numbered as indicated in Fig. 2, 
choose the following coordinates: 


TABLE VII. 

Type Symbol Number Detailed description 
Stretching P; 2 C;—C; 
Stretching Ri; 4 C;-—Hi; 
Stretching Si 2 C;— His 
Stretching T; 2 
Bending 4 His 
Bending Bi Hi—Ci—Hie 
Bending 1 H3i1—C3— Haz 
Bending bij 4 Hij—C:—Cs 
Bending 2 His—C;—Cs; 
Bending bij 4 H;;—C;—C; 
Bending n 1 C,—C.—C; 
Twisting Ti 2 


Count of these coordinates reveals three re- 
dundancies which correspond to the vanishing 
sums of six bending coordinates at each carbon 
atom. The special notation for the elements of 
the G matrix is given below, only the first row 
from each set of symmetrically equivalent coor- 
dinates being required. 


Ti T2 an a2 a@ a2 By Be Y 
gpr £Pa O O gpp gpPy 
0 0 Ra 0 gre O 0 
0 0 gsa 9 O gs O O 
gr er 0 0 0 0 

gla 8a 0 0 Sap 0 0 
28 0 0 
gy 
$12 S22 n T1 T2 
gat Sat 0 Zan Lar 
gar gar 0 0 
wer Zen fier Ber 
8n Snr Snr 
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The zeros in this matrix arise simply because 
the pairs of coordinates in question possess no 
common atoms. The matrix elements g'p,, gp,, 
without detailed calculation, since they represent 
interactions of one of the first two types of coor- 
dinates with the third type and the configura- 
tions of atoms defining them contain a plane of 
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symmetry (g'p, would vanish even without 
symmetry, being of type g*,,(?)). Many of the 
remaining elements are equal: formulas for them 
are obtained by substitution of appropriate 
values for r in Table IV. It is assumed that all 
C—H distances are equal; the reciprocal of such 
a distance is designated by Ro and that of the 
C—C distance by Po. Certain elements, fol- 


TABLE VIII. 
Type G elements in CHsCH:CH; Formula 
rr gip 2uc 
—pc/3 
Pp — (21/3) Rouc 
— (21/3)—P 
LTS — (23/3) Pouc 
(1) £Pa= £PB=&Py= Ra= (23/3) Rouc 
(23/3) Pouc 
&Rr(+) — (34/2) Powe 
(3) =g7r(+) — (34/2) Rowe 
8°46 £8 =8y (2/3) R*o(3uu+4uc) 
(14/3) 
8a = — (1/2) 
—(1/6)[3R% n+ (3Ro+5Po)(Ro—Po)uc] 
ad = Lac = £85 = Lyn — (1/6)[3R* +4Ro(Ro— Po)uc] 
—(1/6)(7Po—4Ro) Pouce 
() (2/3)(3Po+ Ro) Pouc 
gen (1/3) (7Po+Ro)Pouc 
— (1/6)(7P0+ Ro) Pouc 
ge (1/3) 
— (2/3) 
— (1/6) P?ouc 
= = fre (4/3)PoRouc 
Bat = San= Sat = Sys — (2/3)PoRouc 
(2/3)(Pot+Ro)Pouc 
(1/6) (Ro—5Po)Pouec 
Ber (1/6)(Po—5Ro)Pouc 
(1) Lyn — (8/3)PoRouc 
gtr — (2/3)(Pot+Ro)*ue 
Sad = — (4/3)(Pot+Ro)Rouc 
+) (64/8) R°o(3un+4yuc) 
— (64/8) (6Po+Ro)Pouc 
(1) (64/8) 
(64/8) (3Po—2Ro)Pouc 
64(6Po+ Ro) Pouc 
+) (64/2) PoRoue 
Str g', (9/8) [Roun + (2P%o+ Ro) uc 
(9/8)(2Ro— Po) Pouc 


it; 
he 
od 

2, 
n 
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ing 
on 
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0 

0 

gy 
T2 
py 
YT 
br 


1034 J. C. DECIUS 


lowed by a (+), would exhibit alternation of sign 
if the complete G matrix, rather than just the 
first row from each symmetrically equivalent set 
had been written out. 


Vv. CONCLUSIONS 


A few remarks on the limits of applicability of 
these formulas are necessary. In the first place, 
it remains to be seen whether the torsional type 
of coordinate introduced here will prove to be 
practical in the treatment of restricted rotation 
as a small vibration. This is not only a question 
of the magnitude of the potential barrier hinder- 
ing the rotation and the quality of the harmonic 
approximation in the trough of the potential 
energy versus angle of torsion, but also of 
whether or not the vibrational potential function 
is reasonably diagonal in terms of this type of 
force constant. It should be mentioned at this 
point that such torsional frequencies as are 
known at present are usually known or inferred 
to be quite low by comparison with the bending 


Fic. 2. Numbering of atoms 
in C3Hs. 


and stretching frequencies, at least in the case of 


hydrocarbons. Taylor and Pitzer!® have shown 
that an approximate treatment of such cases may 
advantageously be carried through in terms of 
the kinetic energy matrix (not its inverse). 

From a purely kinematical point of view, the 
set of coordinates chosen here is inadequate for 
the treatment of such planar molecules as BF; 
or CH,0 which possess no chain of bonded atoms 
with a length of four atoms; when a structure 
similar to that of such molecules forms a part of 
a larger molecule, the same coordinate deficiency 
will remain. Another type of deficiency of the 
tables occurs when a part of the molecule con- 
tains a chain of 3 or more colinear atoms, since 
one or more dihedral angles become indeter- 
minate. On the other hand, a large number of the 
G matrix elements for many cyclic molecules 
may be found in the tables, including, for ex- 
ample, all necessary in-plane interaction elements 
for benzene. 


10W. J. Taylor and K. S. Pitzer, J. Research Nat. Bur. 
Stand. 38, 1 (1947). 
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The vapor pressures of six alkali halides, NaBr, NaI, KI, CsCl, CsBr, and CsI, have been 
measured in the range from about 10~‘-mm to 0.1-mm pressure. The measurements were made 
by the diffusion method, using an oxygen-coated tungsten filament as detector. 


I. INTRODUCTION 


APOR pressures of solid alkali halides have 
been measured by a number of investi- 
gators. Mayer and Wintner' used Knudsen’s 
method of molecular effusion.? The amount of salt 
Q escaping into a vacuum from an oven at tem- 
perature TJ through a small hole of known area A 
in a definite time interval ¢ is related to the pres- 
sure P in the oven by 


P=(Q/At)(2eRT/M)}, 


where M is the molecular weight of the salt and 
R is the gas constant. The effused salt was 
condensed on a cold surface and the amount de- 
termined by quantitative chemical analysis. 
Niwa* and Dietz‘ used the same method, but 
measured the loss of weight of the oven. Dietz 
also used an electromagnetic balance as an 
absolute manometer. 

Zimm and Mayer® combined Knudsen’s effu- 
sion method with the positive ion method of 
Langmuir.® The effusing salt beam was inter- 
cepted by a tungsten filament whose tempera- 
ture and work function were sufficiently high to 
allow complete dissociation of the salt molecule 
and complete ionization of the metal atom, and 
the current due to the metal ions was measured. 
This method gives directly the rate of effusion 
and thus eliminates the measurement of the time 

* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

** Submitted in partial fulfillment of the requirements 
for the degree of Doctor of Philosophy at Columbia Uni- 
versity. 

1J. E. Mayer and I. H. Wintner, J. Chem. Phys. 6, 
301 (1938). 

2M. Knudsen, Ann. d. Physik 47, 697 (1915). 

3K. Niwa, J. Faculty Sci. Hokkaido Imp. Univ. Ser. 
III, 2, 201 (1938). 

‘V. Dietz, J. Chem. Phys. 4, 575 (1936). 

a : 4a) H. Zimm and J. E. Mayer, J. Chem. Phys. 12, 362 


6 J. B. Taylor and I. Langmuir, Phys. Rev. 51, 753 
(1937). 


interval. It is not necessary to hold the oven 
temperature constant for long periods of time. 
The sensitivity of the detector allows the vapor 
pressure measurements to be extended to lower 
temperatures. 

In this research, the method of Zimm and 
Mayer was used to measure the vapor pressures 
of sodium bromide, sodium iodide, potassium 
iodide, cesium chloride, cesium bromide, and 
cesium iodide. The heats and entropies of va- 
porization were calculated from these data. The 
entropies of the crystals at 800°K and the heats 
of vaporization at absolute zero were also 
calculated. 


II. PRINCIPLES OF THE METHOD 


The number of salt molecules escaping per 
second into a vacuum from an oven at known 
temperature through a hole whose diameter is 
smaller than the mean free path of the molecules 
is given by Knudsen’s formula 


n=SP/(2xmkT)}, (1) 


where S=area of hole, P=pressure inside oven, 
T =absolute temperature of oven, m=mass of 
molecules, and k=Boltzmann constant. The 
molecules leaving the oven are allowed to im- 
pinge upon a filament whose temperature and 
work function are sufficiently high so that every 
molecule hitting it is dissociated and every metal 
atom thus formed is ionized. If sufficient po- 
tential is put on the collector, the number of 
positive ions reaching it will equal the number 
of molecules hitting the filament. Consideration 
of the geometry of the apparatus allows the 
pressure in the oven to be expressed in terms of 
the positive ion current 4p. 


F Vio, (2) 


where r=radius of hole in oven, R=distance 
1035 
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between hole and filament, N»=Avogadro’s 
number, F=Faraday, and A=area of plane 
projection of filament. 

The molecular effusion method allows the 
calculation of the vapor pressure from the effu- 
sion rate, provided the molecular weight of the 
gas is known. Zimm and Mayer showed it to be 
unlikely that the degree of association of KCI is 
greater than 0.3 percent at 800°K. Because the 
properties of the other alkali halide vapors are 


similar, it is assumed that the amount of dimeri-" 


zation is negligible in all cases. The extent of 
dissociation of the molecules into atoms may be 
estimated from the dissociation energy of the 
molecule. In the case of CsI, which has the lowest 
dissociation energy of the alkali halides, the dis- 
sociation would be about 510-5 percent at the 
vapor pressure of 10-? mm at 800°K. 

The pressure of positive ions produced by the 
dissociation of the gas in the oven to ions is 
even less than the dissociation into atoms, and 
is consequently negligible. In the case of Csl, 
the extent of dissociation is about 5X10-® per- 
cent and corresponds to a pressure of cesium 
ions of 5X10-” mm at 800°K. Even considering 


K 


Fic. 1. Schematic diagram of the apparatus. 


that the ions leaving the oven orifice carrying a 
current directly to the large collecting plate 
without requiring collision with the filament so 
that the geometric factor of R?/A in Eq. (2) is 
much larger in this case, the low pressure negates 
any appreciable errors due to this cause. 


Ill. APPARATUS 


A diagram of the apparatus is shown in Fig. 1. 
The oven A in which the salt was heated was a 
4X1.5-cm polished silver cylinder whose wall 
thickness was 0.3 cm; the cover was 0.7 cm 
thick and screwed into the body of the oven to 
a depth of 0.4 cm. The hole £ through which 
the molecules escaped was 0.1057 cm in di- 
ameter, finished with knife-sharp edges. Another 
oven, used for additional runs on cesium bro- 
mide, had a hole 0.0105 cm in diameter. The 
oven was heated by a Nichrome coil B wound 
around a quartz tube C; both were enclosed by 
a silver housing D. Current for heating the coil 
was taken from a Solar regulating transformer. 
The silver-constantan thermocouple JN, used to 
measure the temperature of the oven, was in- 
serted into a hole bored in a silver plug which 
was then flattened to hold the couple securely ; 
the plug screwed deep into the cover of the oven. 
The cold junction O was kept in a bath of 
crushed ice and water; the e.m.f. of the couple 
was measured with a Leeds and Northrup 
Type K potentiometer P. The couple was stand- 
ardized at the freezing point of Bureau of Stand- 
ards samples of tin, zinc, and aluminum; when 
it was later checked at a single point, it showed 
no appreciable deviation from the initial cali- 
bration. The beam of escaping salt molecules, 
defined by the aperture in the collector plate H, 
fell on a tungsten filament F, 0.02 cm in di- 
ameter, which had been previously aged for 
one hour at 2600°K. The filament (13.5-14 cm 
long) was suspended vertically, using a heavy 
tungsten spring to keep it taut when heated. 
The length of the filament exposed to the mole- 
cules was approximately 2 cm. Current for 
heating the filament was supplied by the d.c. 
line K. The voltage drop across the filament 
could be measured by means of platinum poten- 
tial leads welded to its ends. The diameter of 
the filament was initially determined by weigh- 
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TABLE I. Values of (eg—J) for r=99. 


Temperature °K (eg —I) ev 
1200 0.547 
1400 0.638 
1600 0.729 
1800 0.820 


TABLE II. Ionization potentials of the alkali metals. 


Metal Ionization potential ev 
Sodium 5.138 
Potassium . 4.340 
Cesium 3.893 


ing a known length of aged wire. The diameter. 
and temperature of the filament during the run 
were calculated from the tables of Jones and 
Langmuir.’ The distance between the hole in 
the oven and the filament was approximately 
10.50 cm. It was redetermined whenever changes 
were made in the tube elements. The ions leaving 
the filament were collected by the collector plate 
H; stray ions were collected by the guard ring J. 
The plates were attached to the same potential 
source L but connection from the collector was 
through a Leeds and Northrup Type HS galvan- 
ometer M whose sensitivity was 2.3107! 
amp./mm. The beam leaving the oven could be 
interrupted by a magnetically operated shutter J. 

Thermocouple wires were introduced into the 
apparatus through capillary tubes sealed with 
Apiezon wax “W;” all other wires were intro- 
duced by means of tungsten to glass seals. The 
plates were made of hydrogen fired ‘‘A”’ nickel. 
Two auxiliary filaments for outgassing were of 
tungsten. The filament and oven assemblies were 
introduced into the glass envelope by means of 
ground glass joints sealed with Apiezon “W.” 
All stopcocks were greased with Apiezon ‘“‘N.”’ 

The system was evacuated by two single- 
stage mercury diffusion pumps in parallel, with 
mechanical backing. Mercury vapor was elimi- 
nated by use of dry ice traps. Residual gas pres- 
sure was measured with a Distillation Products 
VGI1A ionization gauge which had been cali- 
brated against a McLeod gauge with liquid 
nitrogen traps. 


7H. A. Jones and I. Langmuir, Gen. Elec. Rev. 310 
(1937), 


TABLE III. Dissociation energies of the alkali 
halide molecules. 


Salt Dissociation energy ev 
NaBr 3.85 
Nal 3.16 
KI 3.33 
CsCl 4.38 
CsBr 3.99 
CsI 3.32 


TABLE IV. Experimental conditions and necessary values 
of work function to satisfy r=99, g=99. 


eg for eg for 
= Max. = 
Filament at value at 
temp. (eg)r.99 NMX q=99 
Salt Ty (°K) - ev N X10-"4 at Ty ev 

NaBr 1400 5.78 0.447 —4.76 5.69 
Nal 1250 5.71 4.61 — 3.69 5.64 
KI 1600 5.07 ye | —3.78 4.29 
CsCl 1800 4.71 2.92 —3.73 4.41 
CsBr 1400 4.53 1.07 —4.10 4.83 
CsBr 1800 4.71 1.07 —4.05 3.96 
CsI 1400 4.53 4.92 —3.48 4.38 


IV. OPERATING CONDITIONS 
Scattering 
The relationship 
F jin (2) 


has been derived assuming there is no scattering 
of the salt beam. If the molecules are moving a 
distance R with uniform velocity through a 
scattering gas, the number J of beam molecules 
received per second by a detector will be 


where J) is the number of beam molecules re- 
ceived per second by the detector in the ab- 
sence of scattering gas and L the mean free path. 
The calculation of Z will be uncertain because 
the scattering cross sections are unknown. There- 
fore, the scattering correction J/I) has been 
evaluated experimentally by measuring the posi- 
tive ion current as a function of residual gas 
pressure (at constant oven and filament tempera- 
ture) and extrapolating to zero pressure. The 
value of ip used in Eq. (2) was calculated from 
the observed positive ion current 7 by applying 
the scattering correction I/Io. The pressure of 
residual gas was usually kept low enough so 
that I/Io was greater than 0.85. Within the 
working range, the curves of I/I versus pressure 
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1.0 


0.4 


Ong 


Fic. 2. Ion emission versus 
temperature. 


Filament Temperature (°K) 


for all the salts studied were within 5 percent of 
an average curve, and were reproducible within 
a few percent. 


B. Reaction on the Filament 
We will now consider the reaction 
MX—M++X-+e (in filament) (1) 


which occurs on the surface of the filament to 
determine whether the assumption of complete 
conversion of MX molecules to M* ions made 
in the derivation of Eq. (2) is valid. 


1. Equilibrium of Conversion Reaction 


The M atoms of the incoming MX molecules 
may leave the filament as MX molecules, M 
atoms, and M* ions only. The necessary condi- 
tions for complete conversion of MX to M+ are 
that the equilibrium ratios 


and mu+/NMx 


are both extremely large. The two together are 
a sufficient condition, assuming no kinetic ob- 
stacle to the establishment of equilibrium. These 
ratios are given by 


=3 exp(eg—T)/RT, (3) 


where I=ionization potential of the atom and 
eg=work function of the filament surface. This 
is the Saha-Langmuir equation® assuming zero 


8]. Langmuir and K. H. Kingdon, Proc. Roy. Soc. 
107A, 61 (1925). 


reflection coefficients for both atoms and ions. 
From the equilibrium constant of reaction (1) 
and the principle of conservation of matter one 
obtains: 


= ] 
(4’) 


where K =equilibrium constant for reaction (1), 
K=exp[(ee—D—I)/RT+ (Sm+°+Sx°+Sux°)/R 
—}], and w=reduced mass of MX. If the ex- 
tent of reaction is so great that nm*+>mmx, then 
the number of M+ ions leaving a unit area of 
filament per second, m+, approximately equals 
the number of MX molecules striking a unit 
area of filament per second, Vmx. 


(4) 


A necessary and sufficient condition for 99 per- 
cent .conversion of MX to Mt? is that both 
r>99 and q>99. 

Table I lists the quantity (eg—J) which satis- 
fies the relation r=99 at various temperatures. 
It can be seen that the work function necessary 
to ionize 99 percent of the metal atoms increases 
with increasing temperature. 

At a given pressure of salt molecules at the 
filament, the ratio g depends on the equilibrium 
constant for the over-all reaction K, which is a 
function of the ionization potential of the metal 
atom I, dissociation energy of the salt molecule 
D, work function of the filament eg, temperature 
T, and the standard state entropy change AS°. 
For any molecule D and J are fixed and AS° 
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TABLE V. Experimental data. TABLE V.—Continued. . 


Sodium bromide Sodium iodide Cesium iodide 
—log Pmm (1000/T) (°K) —log Pum (1000/T) (°K) —log Pum (1000/T) (°K) —log (1000/T) (°K) 
1.77 1.098 0.49 1.108 1.50 4.172 2.85 1.295 
1.91 1.114 1.11 1.160 1.53 1.173 2.85 1.297 
2.05 1.129 1.20 1.169 1.59 1.179 2.85 1.307 
2.32 1.153 1.49 1.190 1.61 1.177 2.87 1.296 
2.62 1.179 1.83 1.221 1.63 1.179 2.97 1.305 
2.67 1.182 2.03 1.244 1.65 1.183 3.01 1.306 
2.82 1.197 2.13 1.253 1.66 1.182 3.04 1.323 
2.87 1.201 2.17 1.255 1.71 1.186 3.07 1.321 
3.15 1.228 2.67 1.300 1.78 1.197 3.07 1.323 
3.73 1.278 2.75 1.311 1.88 1.203 3.08 1.323 
3.98 1.301 3.27 1.358 1.88 1.209 3.13 1.320 
4.20 1.319 1.89 1.206 3.24 1.339 
Potassium iodide Cesium chloride 0.59 1.118 
fen CR) (1080/7) (EK) 0.62 1.120 2.07 1.254 
0.90 1.114 0.89 1.122 0.70 1.128 2.29 1.265 
1.24 1.144 1.12 1.141 0.71 1.127 2.32 1.273 
1.37 1.152 1.16 1.146 0.96 1.147 2.58 1.305 
1.39 1.153 1.58 1.184 0.99 1.151 2.63 1.310 
1.46 1.162 1.60 1.185 1.32 1.182 2.73 1.315 
1.63 1.177 1.73 1.200 1.35 1.179 2.85 1.327 
1.67 1.179 1.85 . 1.209 1.50 1.198 3.12 1.353 
1.69 1.179 2.03 1.223 1.53 1.195 3.26 1.371 
1.95 1.207 9 37 1.256 1.70 1.214 3.54 1.398 
1.98 1.209 2.40 1.263 1.85 1.234 3.68 | 1.418 
2.00 1.204 2.61 1.281 2.07 1.245 3.85 1.437 
2.00 1.209 2.71 1.290 
2.02 ae 2.78 1.298 
varies slowly with temperature; K will therefore 
vary with and eg. If it were possible to in- 
263 1.262 3.33 1.354. crease the work function of the surface without 
2.63 1.264 3.43 1.363 limit, the equilibrium constant at any tempera- 
2.93 1.301 3.55 1.379 . ‘ 
294 1.297 3.58 1.381 ture would be sufficiently large to satisfy the 
2.94 1.302 3.76 1.397 condition g=99. The actual existence of an 
2.95 1.297 3.85 1.406 
309 1315 3.88 1409 upper limit to the work function means that 
there will be a minimum filament temperature 
3.55 1.363 4.14 1.435 below which the conversion of MX to M+ will be 
3.81 1.383 4.25 1.449 incomplete. 
3.84 1.386 
Cesium bromide 0.8- 
—log pmm (1000/T) (°K) —log Pmm (1000/T) (°K) 
0.92 1.111 1.92 1.202 
1.00 1.120 1.95 1.210 Filament temperature 
1.02 1.121 1.95 1.210 
23 1.141 1.97 1.215 
1.21 1.141 2.00 1.209 
pes 1.145 2.00 1.218 
1.23 1.145 2.04 1.220 ) 
1.24 1.144 2.13 1.230 
1.24 1.146 2.23 1.236 0. 
1.26 1.144 2.25 1.237 
1.26 1.148 2.25 1.239 
1.27 1.153 2.37 1.243 4 
1.31 1.150 2.38 1.254 
1.38 1.158 2.39 1.255 oO 
1.39 1.161 2.42 1.266 
1.40 1.159 2.48 1.255 
1.40 1.159 2.50 1.261 
1.44 1.162 2.55 1.271 
1.44 1.165 2.66 1.277 09 + 
1.47 1.164 2.69 1.277 600 700 800 ee 1000 
1.166 1.281 ture 
149 2:70 1281 Pilenens tenpevature 
1.49 1.170 2.84 1.291 Fic, 3. Adsorption of cesium on tungsten (de Boer). 
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1.0_ 

Sodium Bromide |, 
t Py, V8. 1000( 
+ 
*g + Mayer and Wintner 

2.0L 
Pam 

Fic. 4. Vapor pressure ver- 
sus 1000/T°K (NaBr). 
4.0. 
1.0 1.1 1.2 1.3 


1000( °k 
) 


Ionization potentials listed in Table II are 
taken from Herzberg.’ Dissociation energies 
tabulated in Table III are from Herzberg!® or 
Mayer and Helmholz." 

The experimental filament temperatures for 
each salt are given in Table IV. Column 3 of this 
table lists (e¢),-99, the work function necessary 
for 99 percent ionization at 7;; in column 6 are 
tabulated (eg) -99, the minimum work function 
satisfying the relation g=99 at 7; and the high- 
est pressure of incoming MX molecules. The 
larger of these two determines the lowest work 
function sufficient for 99 percent conversion. 

9G. Herzberg, Atomic Spectra and Atomic Structure 
(Dover Publications, New York, 1944). 

10G. Herzberg, Molecular Spectra and Molecular Struc- 
ture (Prentice-Hall, Inc., New York, 1939). 


J. E. Mayer and L. Helmholtz, Zeits. f. Physik 75, 
19 (1932), 


Under experimental conditions it was possible 
to maintain a work function of about six volts, 
which is higher than any of the required work 
functions. Thus, equilibrium favors conversion 
of every molecule which hits the filament into 
the appropriate ion. 


2. Kinetics of the Conversion Reaction 


The following assumptions are made: 

a. The rate of reaction is so fast that the re- 
action goes virtually to completion. There are 
three possibilities : 


(i) The rate is very slow; no conversion. 
(ii) There is partial conversion. 
(iii) The reaction is very rapid; complete 
conversion, 
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The existence of a positive ion current eliminates 
the first possibility. If the conversion were in- 
complete, in all probability it would have a 
temperature coefficient. Only a reaction which 
had gone virtually to completion would not, 
since a change from 99.0 to 99.9 percent con- 
version is a tenfold change without change in 
result. A plot of the experimental ion current 
divided by the maximum ion current versus 
filament temperature for KI is shown in curve J 
of Fig. 2. The existence of a plateau indicates 
that the reaction is actually very fast in this 
temperature interval. 

b. The rate of removal of products is rapid 
enough to prevent accumulation of metal ions 
on the surface, with consequent reduction of the 
work function. When the equilibrium requires 


incomplete conversion, one may calculate the 
number of metal ions formed at any given tem- 
perature and work function (i.e., for any value 
of K and r) for any value of the pressure of in- 
coming molecules. One arrives at an expression 
for mm*: 


][r/(A+r)]=0. (5) 
Thus a=7y+/Nmx, the fraction of MX mole- 
cules which are converted to M* ions under 
these conditions, may be calculated. a for KI 
has been evaluated as a function of temperature, 
assuming a constant work function eg=5.15 
volts, the work function required for r>99 at 


1800°K. The oxygen coating strips at high tem- 


peratures, lowering the work function. Under 
the conditions of these experiments, stripping 


Sodium Iodide ° 
~10g Phy VS. K) 
1.0L 
Pon 
2.01. 
Fic. 5. Vapor pressure versus 
1000/T°K (Nal). 
320 |_ 
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did not occur under 2000°K. The results of this 
calculation are plotted in curve II of Fig. 2. 

A comparison of the experimental curve I and 
the theoretical curve II in Fig. 2 shows that in 
the case of KI the shift in the equilibrium does 
not account for the sharp drop in ion current at 
low filament temperatures. This phenomenon is 
probably due to the fact that desorption of 
metal ions is slow at low temperatures. Since 
the adsorption energy for metal ions on an oxi- 
dized tungsten filament is probably very large, 
the layer of metal ions will disappear in a very 
small temperature range above some critical 
temperature, and the work function of the 
filament will increase abruptly. This is illustrated 
by Fig. 3 taken from de Boer! which shows the 

123. H. de Boer, Electron Emission and Adsorption 


(Cambridge University Press, Teddington, 


variation of 6, the fraction of tungsten surface 
covered by adsorbed cesium, with increasing 
filament temperature, when the cesium is de- 
sorbed as ions. According to Taylor and Lang- 
muir,® with a given pressure of cesium and with 
a sufficiently strong accelerating field for ions to 
overcome space charge there is a critical tem- 
perature T., above which every cesium atom 
that strikes the filament escapes as an ion; 
below 7,, ion emission is very small. The plateau 
in curve I of Fig. 2 implies a high enough work 
function and a rapid removal of metal ions in 
this temperature range. 

Therefore it is possible to choose the experi- 
mental filament conditions so that the assump- 
tion of complete conversion of MX to Mt on 
the filament is valid. In these experiments, 
precautions were taken ‘to have sufficient oxygen 


if 


e 
+ No 
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present (using a capillary leak when necessary) 
to coat the filament and increase its work func- 
tion. For each salt, the positive ion current at 
some constant oven temperature was measured 
as a function of filament temperature, and the 
vapor pressure measurements were made in the 
temperature range of the plateau. 

c. In order to measure the positive ion cur- 
rent, a potential of —45 to —67.5 volts was put 
on the plates to overcome ion space charge. To 
eliminate the effects of stray currents, the gal- 
vanometer deflection was first measured with 
the shutter closed, and the additional deflection 
when the shutter was opened was assumed to be 
due entirely to the positive ions. The deflections 
with closed shutter rarely exceeded ten percent 
of the ion deflection. 


V. MATERIALS 


The sodium and potassium salts were prepared 
by recrystallizing C.P. grade salts from distilled 
water and drying them in a vacuum oven. The 
cesium salts were prepared by adding a small 
excess of the appropriate redistilled C.P. acid 
to C.P. cesium carbonate (Maywood Chemical 
Works) which had been recrystallized from 
absolute alcohol, precipitating the salt, and 
drying it in a vacuum oven. 

The presence of a highly volatile impurity in 
the salt would give a very high apparent vapor 
pressure, but its concentration would change 
during the course of the run, and the resulting 
pressure would drift with time. The results of 
measurements on the same sample made on 
different days agreed with each other within the 
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experimental error, and low temperature data 
could be repeated after keeping the sample at 
high temperatures for several hours. 


VI. PRIMARY RESULTS 


The experimental data are tabulated in Table 
V and plotted in Figs. 4 and 9. The data are well 
fitted by the equation: 
logPmm = C—A(1000/T) 

+5/2 log(1000/7). (6) 
The term 5/2 log(1000/7) arises from an esti- 
mate of the difference of heat capacities between 
the gas and crystal. The exact value of the AC, 
term is unimportant, since it represents only a 
small correction term. An average value of 
AC,=-—5/2R for the temperature range 500°K 
to 1000°K has been used. 

The constants in Eq. (6) have been evaluated 


LoL 


3-0}. 


4,0 


by the method of least squares and are tabulated 
in Table VI. The curves in Figs. 4 to 9 have 
been drawn using the constants in Eq. (6). 


VII. SOURCES OF ERROR 


The sources of error in p are: 


A. Geometry: +2.9 percent; distance from 
hole to filament: R+0.5 percent; radius 
of hole: r+1.1 percent; length of fila- 
ment: 11.0 percent; diameter of fila- 
ment: d+1.5 percent. 

B. Ion current: +2.8 percent; measure- 
ment of galvanometer deflection: +2 
percent ; incomplete reaction: +0.3 per- 
cent; scattering correction: +2 percent. 

. Temperature: +0.2 percent; maximum 
error: +2 percent; probable error: 0.2 
percent. 
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The probable error in p is +4.0 percent. There- 
fore, the probable error in Inp is +0.040. 


VIII. DISCUSSION OF RESULTS 


Cesium chloride exjsts in two crystalline modi- 
fications. The low temperature form is a body- 
centered lattice; at 733°K x-ray data!" show 
a transition to a face-centered lattice. It is ex- 
pected that the energy difference between these 
two cubic forms is small. From x-ray measure- 
ments, the heat of transition AH; is estimated as 
1.8 kcal." Thermal data give a AH; value of 1.34 
kcal.!® There is no observable break in the vapor 
pressure curve of cesium chloride (see Fig. 7). 


8G, Wagner and L. Lippert, Zeits. f. physik. Chemie 
B31, 263 (1936). 

“C. D. West, Zeits. f. Krist. 88, 94 (1934). 

15 International Critical Tables 5, 207. 


1.3 1.4 
1000( °x) 
T 


The probable error in the slope of the line in- 
cluding only those points below the transition 
temperature is 3.1 percent. The accuracy of 
these measurements would only be sufficient to 
say that the heat of transition probably does 
not exceed 3.1 percent of AHyap or 1.5 kcal. 
Additional data on the vapor pressure of 
cesium bromide which have not been used in 


TABLE VI. Constants to be used in Eq. (6). 


Standard 
deviation 
Salt A of A Cc 

NaBr 12.003 0.083 11.347 
Nal 11.895 0.14 12.525 
KI 11.455 0.097 11.660 
CsCl 11.076 0.062 11.346 
CsBr 11.421 0.079 11.676 
CsI 11.193 0.066 11.726 
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TABLE VII. Heat content in kcal./mole. 


TABLE VIII. Entropies in cal./degree. 


wo* 
Hegas wave Herystal 
Salt AHs00 T =800 no. T =800 AHo 


NaBr 50.93+0.35 6.79 133.9 9.38 53.52 
Nal 50.43+0.91 6.83 117.0 9.65 53.25 
KI 48.42+0.46 6.92 98.0 9.65 51.15 
CsCl 46.69+0.26 6.90 98.0 9.76 49.55 
CsBr 48.26+0.33 6.97 74.6 10.21 51.50 
CsI 47.23+0.28 7.01 (61)* 9.82 50.04 


Sgas° —In Pmm Sgas Scrystal 
Salt T =800 T =800 T =800 ASs00 T =800 


NaBr 79.60 7.86 95.21 63.66 31.55 
Nal 81.33 4.84 90.94 63.04 27.90 
KI 83.55 5.57 94.60 60.53 34.07 
CsCl 83.17 5.20 93.50 58.36 35.14 
CsBr 85.71 5.43 96.49 60.35 36.14 
CsI 87.53 4.65 96.81 59.04 37.77 


* Estimated. 
*R. B. Barnes, Zeits. f. Physik 75, 723 (1932). 


the evaluation of the constants in Eq. (6) are 
plotted in Fig. 10. These peculiar results are not 
explanable. Two ovens were used in these meas- 
urements (Oven I with a hole diameter of ap- 
proximately 0.1 cm and Oven II with a hole di- 
ameter of approximately 0.01 cm) and the 
data were taken at two filament temperatures 


(1400°K and 1800°K) using many different sam- 
ples of salt. 

The portion of curve I labeled AB corre- 
sponds to Fig. 8 and was obtained in both ovens 
with fresh samples which had never been heated 
above 7, (approximately 875°K). The material 
looked powdered before and after heating. BC 
shows the rapid rise in vapor pressure when the 
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salt was heated above 7, in Oven I. The data 
then fall on curve II when the salt is cooled and 
remain on curve II throughout all subsequent 
heating and cooling. This material now appears 
sintered, but there is no indication that the 
sample had melted (Tmn.p.=909.2°K). No dif- 
ference in structure was found from x-ray powder 
pictures of the original powdered salt, crushed 
sintered material, and small (but uncrushed) 
fragments of sintered material. Previous x-ray 
investigators had looked for a transition in 
cesium bromide but had never found one. 
Microchemical analysis showed these materials 
to be the same within the error of the determina- 
tion (original sample: 37.55 percent Br; sintered 
material: 37.98 percent Br). BD was obtained 
when fresh samples were heated above 7, in 
Oven II and is a continuation of AB. This ma- 
terial was still powdered after heating. DE 
shows the extremely rapid rise at 900°K (an 
increase in pressure from 0.1 to 7 mm in 2.5°) 
when the sample was heated in Oven II. On 
cooling the data could be represented by curve 
II and the material was sintered. There is no 
evidence that the salt, which leaves the oven at 
a very great rate in the region of anomalous 
pressures, has condensed on the outer surface 
of the oven or radiation shields. 

The data of Mayer and Wintner for sodium 
bromide are plotted in Fig. 4. In Fig. 6 the re- 
sults of Mayer and Zimm for potassium iodide 
are given. Their entropy of vaporization agrees 
with ours to 0.5 percent, but the vapor pressures 
obtained in these measurements are 1.55 times 
those of Zimm and Mayer. Results of measure- 
ments on two samples of potassium iodide, one 
prepared in this laboratory and the other in the 
phosphor research laboratory at the Polytechnic 
Institute of Brooklyn, were the same within 
experimental error. 


TABLE IX. Entropies in cal./degree. 


Seryat(T =800) Serystal Sgas°? 
Salt —Seryst(T =298) T =298 T =298 
NaBr 12.76 18.79 70.85 
Nal 13.15 14.74 72.61 
KI 12.92 21.15 74.80 
CsCl 13.10 22.04 74.35 
CsBr 13.73 22.41 76.87 
CsI 12.80 24.97 78.72 


IX, DERIVED RESULTS 


The measurements of this research permit the 
calculation of AF 7° for the vaporization of the 
salt (from the vapor pressure at any one tem- 
perature) and the AH; of vaporization from the 
relation : 


d Inp/d(1/T) = —AHp/R. 


From these one may calculate the entropy dif- 
ference between vapor and crystal. The heats of 
vaporization at 800°K are listed in Table VII; 
the entropies of vaporization at 800°K are given 
in Table VIII. 

The entropies of the crystal and gas may be 
calculated from statistical equations!® but with 
uncertain accuracy. The equilibrium inter- 
nuclear distance in the gaseous. molecule is 
known in all cases.!” Since the frequencies of the 
gaseous molecule are known in some cases,® and 
a reasonably good estimate of the frequencies 
has been made in others,'® the accuracy with 
which the vapor entropies may be obtained is 
good. The calculation of the crystal entropies, 
however, is more uncertain. In order to get good 
agreement between crystal entropies obtained 
from heat capacity data and those calculated 
using the Debye function, it is necessary to use 
the reststrahlen frequency to determine the pa- 
rameter 0p. This is theoretically unsound. The 
use of this as an empirical relationship is limited 
because few of these reflection frequencies are 
known. The Nernst-Lindemann equation has a 
sounder theoretical basis, but does not give as 
good agreement with experiment. Therefore, the 
entropy of the crystal was obtained by using 
the experimental entropy difference between the 
gas and crystal and the calculated vapor entropy. 
The values of the entropies of the gas and crystal 
at 800°K are tabulated in Table VIII. Because 
heat capacity data are lacking above room 
temperature for most of these crystals, ex- 
trapolation of the crystal entropy from 800°K 
to 298°K has been made, using the heat capacity 
equations estimated by Kelley.'® In the case of 
sodium bromide, a long extrapolation of the 


16 J. E. Mayer and M. G. Mayer, Statistical Mechanics 
(John Wiley and Sons, Inc., New York, 1940). 

17E. J. W. Verwey and J. H. Boer, Rec. trav. chem. 59, 
633 (1940). 

18 W. Gordy, J. Chem. Phys. 14, 305 (1946). 

19K, K. Kelley, U. S. Bur. Mines Bull., 371 (1934). 
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TABLE Al. Quantities used in calculating gaseous entropies: I.=yuro?, Be=h/8x°I.c, 


Inter- Equi- 


Reduced Fre- Anhar- atomic* librium>’ Moment of 
mass quency monicity distance distance inertia 
We correction 71500 ro Ie =Be/we 
Salt grams WeXe cm cm BeX102 uT =hew*/k oT =Bohc/k yX10 


Sodium bromide 2.965 315» 1.15 2.64 2.59 1.989 14.07 450.4 0.2019 4.47 
Sodium iodide 3.232 286° 0.75» 2.90 2.84 2.607 10.74 409.4 0.1543 3.76 


Potassium iodide 4.962 212 0.7 3.23 3.16 4.955 5.649 302.9 0.08110 2.66 
Cesium chloride 4.643 226° (0.7)* 3.06 3.00 4.179 6.698 323.8 0.09620 2.96 
Cesium bromide 8.278 168° (0.5)* 3.14 3.07 7.802 3.588 240.3 0.05154 2.13 
Cesium iodide 10.77 137¢ (0.4)* 3.41 3.33 11.94 2.344 195.7 0.03367 ° 1.71 


* Estimated. 


aL. R. Maxwell, S. B. Hendricks, and V. M. Mosley, Phys. Rev. 52, 968 (1937). 


b See reference 10. 
© See reference 18. 


TABLE A2. S,°=translational entropy at 1 mm pres- 
sure, S,=vibrational entropy, 5S,;=rotational entropy, 
S-=correction to the entropy due to anharmonicity and 
stretching. 


T =800°K 
‘Salt Se°/R S/R Se/R S°/R 
NaBr 29.144 1.588 9.286 0.058 40.076 
Nal 29.698 1.681 9.556 0.019 40.954 
KI 29.860 1.977 10.198 0.032 42.067 


CsCl 29.881 1.912 10.028 0.060 41.881 
CsBr 30.233 2.207 10.652 0.067 43.159 
CsI 30.533 2.410 11.077 0.079 44.099 


T =298°K 
Salt Sk°/R S/R Se/R S°/R 
NaBr 26.673 0.688 8.298 0.015 35.674 
Nal 27.234 0.757 8.568 0.003 36.562 
KI 27.397 1.019 9.256 0.010 37.662 


CsCl 27.403 0.964 9.040 0.029 37.436 
CsBr 27.760 1.242 9.664 0.040 38.706 
CsI 28.059 1.439 10.090 0.048 39.636 


experimental heat capacity equation was neces- 
sary. The entropy of the crystal and the calcu- 
lated entropy of the gas at 298°K are given in 
Table IX. 

Because sufficient heat capacity data are un- 
available, it was necessary to calculate the heat 
content of the crystal. The heat content of the 
crystal at 7 = 298 was calculated in the following 
way: 


J 2 
Hr—Ho 
0 


where the last term arises from the Lindemann- 


Magnus approximation for C,—C,. The Nernst- 
Lindemann heat capacity equation was used to 
evaluate C,: 


Cy = (3nR/2)f2(92/2T), 


where fz = Einstein function, 0g =hvo/k, vo = opti- 
cal absorption frequency, and »=number of 
atoms per molecule of crystal. AH for the in- 
terval T=298 to T=800 was evaluated using 
the heat capacity equations of Kelley. The values 
of H.=H,—H are tabulated in Table VII. The 
difference between the heat content of the crystal 
and the gas at 800°K added to the experimental 
heat of vaporization at 800°K gives AH, the 
heat of vaporization at absolute zero. These 
results are given in Table VII. 
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APPENDIX 


The entropy of the gaseous alkali halide mole- 
cules at 298°K and 800°K have been evaluated 
by the usual statistical mechanical equations.'® 
The quantities used in these calculations are 
listed in Table A.1. In Table A.2 are given the 
various contributions to the entropy and the 
total standard state entropy (pressure 1mm). 
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Previous calculations of integrals which enter computations by the method of antisymmetric 
molecular orbitals are reviewed, several new integrals are evaluated, and a master table of 


integrals is given. 


N a paper entitled ““Antisymmetric molecular 
orbitals. I]. Calculation of integrals,’ Sklar 
and Lyddane! summarized the evaluation of 


some of the integrals over atomic orbitals which ~ 


occur in quantitative calculations by the method 
of antisymmetric molecular orbitals. In the 
present paper, the work of Sklar and Lyddane 
and previous investigators will be critically 
reviewed, and certain new integrals will be 
evaluated. 

The integrals to be considered involve 2p 
atomic orbitals on two identical nuclei: A and B. 
If these nuclei are taken as origins for the 
Cartesian coordinate systems (%a,Va,Za) and 
[Xa=Xv, 2's along internuclear 
axis ], and the polar coordinate systems (a,0q,¢) 
and (b,6,,¢), respectively, the three orthonormal, 
real, hydrogen-like 2 wave functions for electron 
von A and B may be written? in atomic units’ as 


a(v) =(Z5/32m)*y. exp(—Za,/2) 
= (Z5/327)!a, exp(—Za,/2) cos¢,, 
a,(v) exp(—Za,/2) (1) 
= (Z5/32m)!a, exp(—Za,/2) sinOa sing,, 
a,(v) = (Z5/327)'z, exp(— Za,/2) 
= (Z5/322)!a, exp(—Za,/2) 


and 


* This paper is based partly on a thesis submitted by 
Robert G. Parr to the Graduate School of the University 
of Minnesota, August, 1947, in partial fulfillment of the 
requirements for the degree of Doctor of Philosophy, and 
partly on research supported by the Office of Naval 
Research, U. S. Navy, through Contract N5ori-147 with 
the University of Minnesota. 

** Present address: Department of Chemistry, Carnegie 
Institute of Technology, Pittsburgh, Pennsylvania. 

1A. L. Sklar and R. H. Lyddane, J. Chem. Phys. 7, 
374 (1939). 

*Eyring, Walter, and Kimball, Quantum Chemistry 
(John Wiley and Sons, Inc., New York, 1944), pp. 89-90. 

3 Distances in units of the first Bohr radius a9=0.5292A; 
energies in units of e?/2a9=2.1792 X10-" ergs = 13.602 ev, 
where e is the charge on an electron. 
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b(v) = exp(—Zb,/2) 
= (Z*/327)'b, exp(—Zb,/2) sin®,, cos¢,, 
b.(v) = (Z°/327)'x, exp(—Zb,/2) 
= (Z5/32m)'b, exp(—Zb,/2) sind, sing,, (2) 
b,(v) = exp(—Zb,/2) 
= (Z5/327)'b, exp(—Zb,/2) cosy; 
where Z is an effective nuclear charge and the 
subscript v denotes coordinates of electron ». 
The integrals are of several types. First, there 
are the overlap integrals‘ 


(3) 
where a is one of {a,a,,a,} and 8 one of {6,b,,b,}. 
Secondly, there are the Coulomb integrals 


= f (4) 


the hybrid Coulomb-exchange integrals 
(2181;8283) 


f (5) 


and the exchange integrals 
(a181;a282) 


f (6) 


where the a; are members of the set {a,a,,a,}, 
the 8; are members of the set {b,b,b,}, and rie is 
the distance between electrons 1 and 2. For zero 
internuclear distance, these integrals reduce to 

4 dv, denotes element of volume in coordinates of elec- 


tron v; unless otherwise specified, integrations are over all 
space. 
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mononuclear integrals of the form 


= f (7) 


Finally, there are the Coulomb penetration inte- 
grals 


(A:68) = f (8) 
and the exchange penetration integrals 
(A:ap) = — f (9) 


where H,*(1) is the spherically symmetric po- 
tential energy operator for electron 1 in the field 
of neutral atom A. 

The various integrals are given in Table I as 
function of the screening constant Z and the 
parameter 


p=ZR, (10) 


where R is the internuclear distance; their 
evaluation is described in detail below. 


R. G. PARR AND B. 


L. CRAWFORD, JR. 


OVERLAP INTEGRALS 


The only non-vanishing integrals of the form 
(3) are and To evaluate these, 
one may introduce elliptical coordinates 
h=(a+b)/R, pw=(a—b)/R, dasabove. (11) 
One then finds® 
Sav =exp(—p/2)[1+ (p/2) 

+ (p?/10)+ (9°/120)], (12) 


and 


=exp(—p/2)[1+ (p/2)+ (p?/20) 
— — (p*/240) J. (13) 


Values of Sa, and Sa,b, calculated from these 
formulas are given in Table I. 
MONONUCLEAR INTEGRALS 


Sklar and Lyddane gave to the integral 
(aa;aa) = f (14) 


which represents the mutual potential energy 
of two 2p electrons in the same 2 orbital on 
the same atom, the value 0.3881Z. Actually, 
Z—(aa;aa) = 0.3914,* as will now be shown. 


TABLE I. Values of integrals (atomic units®). 


Integral p 0 3e 4c 5e 6 7 8 8.37 9 10 12 14.5 16.7 20 
Z“(aa;bb) 0.3914 0.3476 0.3151 0.2890 0.2616 0.2367 0.2163 «0.20868 =0.1972 0.1806 0.155 0.131 0.115 0.097 
Z-\(aa;babs) 0.3492 = 0.3164 0.2909 «0.2710 (0.2484 0.2273 0.2097 0.20368 0.1928 0.1774 
Z-\(aas;bbs) 0.0211 0.0156 0.0121 0.0090 0.0066 0.0047 =: =0.0022 0.0016 


Z"\(aa;brbr) 0.3492 0.338 0.324 0.305 0.285 0.262 


Z"\(aar;bbr) 0.0211 0.0056 —0.0009 —0.0044 —0.0064 —0.0067 —0.0062 —0.0053 —0.0042 
Z-\(araribrbr) 0.3914 0.359 0.346 0.332 0.315 0.296 0.278 0.252 0.231 
Z-\(ab;bb) 0.3914 : 0.1538 0.0766 0.0448 0.0087 
Z"\(ab;bebs) 0.3492 0.1399 0.0704 0.0415 0.0082 
Z-\(abs;bbs) 0.0211 0.0069 0.0031 0.0016 0.0003 
Z-\(ab;ab) 0.3914 0.073 0.022 0.0082 0.0031 0.0006 0.0001 0.0000 
Z-\(ab;asbs) 0.3492 0.066 0.020> 0.0074» 
Z"\(abs;abs) 0.0211 0.005 0.005 0.0035 0.0025 0.0012> 0.0004 
Z-"\(ab;;abr) 0.0211 0.043 0.030 0.039 0.028 0.021 0.015 0.010 0.006 
Z\(abrjarb) —0.0211 —0.048 —0.030 —0.0389 —0.028 —0.021 —0.015 —0.010 —0.006 
Z-1(A:bb)4 0.5859 0.0773 0.01976 0.00702 0.00182 0.00031 0.046 0.0546 
Z-\(A:ab)4 0.5859 0.1200 0.04364 0.02090 
Sab 1.0000 0.46800 0.25995 0.16396 0.08874 0.03883 0.01798 
Sayby 1.0000 —0.15356 —0.08946 

Interpolated. 

tn dightly to the trensilen to the 0 values. Such ould be i 

© Values in these columns probably y to sm e ition e limiting p= ues. ustment w accord with the state- 
522 of reference 5) that his values may be somewhat in error for small p. 

=carbon. 
© See reference 3. 


5 J. H. Bartlett, Phys. Rev. 37, 507 (1931). 


®In a letter to one of the authors (RGP) dated August 30, 1947, Dr. Lyddane agrees to this value. 
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INTEGRALS IN ORBITAL CALCULATIONS 


In polar coordinates (15) reads the integral over 6; and @a: is*® 
(aa;aa) = (Z5/32n)? f exp(—Za,a;" f 00, 
1 2 
Xsin?6a; (2/112) -| f (1—x?)P,,!™! (dr | 
Xexp(—Zad2)a2? sin*Oaz cos*podve. (15) 
Now? (16/9 if m=0 and n=0, 


16/225 if m=0 and n=2, 

=4 256/25 if |m| =2 and m=2, (22) 
0 if m=0 and 10 or 2, 

0 if |m| =2 and n#2. 


(2/ri2) =2 x (n= |m|) Cont 
- (cos6a;) 


xX P,.'"! explim(¢1—¢2) ], (16) 


where Hence 


a<=min.{@1,d2}, (17) (aa;aa) =2(24/32)+ (16/9) f exp(—Za) 
and the quantities 
(18) (4/4) (2/41) (256/25)] f exp(~ Za) 
are the associated Legendre functions. Also, 


Xexp( [ (16/225) 


dv, =a,’ sin ( 19) exp( ~ 
Therefore 
(aa;aa) =2(25/32n)? (n—|m|)! = (Z/7200) | 25 f exp(—u)u‘gs, 1(u)du 
n=) m=n 0 


f COS": +4 f sud, (23) 
0 


Xexp[im(¢1— $2) where 
x f (cos@a1) 25,1(u) = (1/u) f exp(— V) VidV 
0 
a d a d a f 
X P,,'™! (cos@a2)d0a;d exp(—Za;) +f exp(— V) 
The integral over ¢; and is 87, 3(u) = (1/u*) f exp(— V) 
0 
J cos’, exp[im(d1—¢2) +42 f exp(— V) VdV. 
2 But 
-| f cos’, 
0 8s, 1(u) = (1/u) {24—exp(—u)[u5+ 6u? 
+18u+24]}, (25) 
gr, a(t) = 144 —exp(—u)[u8+ 6u! 
if |m| =2, (21) + 24u'+ 72u?+ 144u-+ 144 }}, 
0 otherwise ; 5See, for example, Whittaker and Watson, Modern 


Analysis (The Macmillan Company, New York, 1943), 


7 Reference 3, Appendix V. p. 324. 
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so that 


f exp(—u)u'gs, 1(u)du 837/8, 
(26) 
f 3(u)du =405/8, 
and 
(aa;aa) = (4509/11520)Z =0.3914Z. (27) 


The integrals (aa;a,a,) =(aa;a,a,) and (aa,;aa,) 
=(aa,;aa,) can be found similarly. These differ 
from (aa;aa) only in the integral over ¢; and ¢»; 
in the case of (aa;a,a,) there enters 


f cos’; exp — $2) 
-| f cosméndés | 
0 


x| f 
0 


mr if m=0, 
— 79/4 if |m| =2, (28) 


0 otherwise ; 


in the case of (aa,;aa,) there enters 


f sing: Cosd, Sings Cosh: exp[im (di — $2) 


2 
=| (1/2) f sin2¢, 
0 


if |m| =2, 
= (29) 
0 otherwise. 


Thus 


(aa;a,a,) =(Z /7200)|25 f exp(u)u'gs, 1(u)du 
0 
—2 4 7.3 d 30 
J exp(— 1) u4gz, 3(u) u| (30) 


(aa,;aa,) = (Z/7200) | exp(— 3(u)du 
0 


L. CRAWFORD, JR. 


and 
(aa;a,a,) = (4023/11520)Z =0.3492Z, 
(aa,;aa,) = (243/11520)Z =0.0211Z. 


(27) and (31) establish most of the limiting p=0 
values which are listed in Table I. 


COULOMB AND EXCHANGE INTEGRALS 


The non-zero Coulomb integrals are (aa;bb), 
(aa;b.b,), (aa;b,b,), (aa,;bb.), (aa,;bb,) and 
(a,a,;b,b,); the non-zero exchange integrals are 
(ab;ab), (ab,;ab.) =(ab,;a.b), (ab,;ab,), (ab;a.b.), 
(ab;a,b,), (ab,;a,b) and (a,b,;a,b,). These may be 
expressed as follows in terms of the integrals 7, 
of Bartlett = 


(aa;bb) =i1+(1/2)%7, (abj;ab) (1/2)ts; 
(aa;b,b.) (1/2)27, (ab;a.b,) =1,— (1/2)is 
(aa,;bb,) = (1/2)tz, 

(ab,;ab.) = (ab,;a,b) =(1/2)is; (32) 
=%3, (ab;arb+) =%4; 
(aa,;bb,) =%5, (ab,;ab,)=%6, (ab,;a,b) =111; 
=%9, = 110. 
The first pair of these formulas was given by 
Goeppert-Mayer and Sklar.’ 

Table | gives values of the Coulomb inte- 
grals (a1a2;8:82) and those exchange integrals 
(a181;a282) which are multiples if 7», 73, and 411, 
calculated from Bartlett’s values for the 7,.5 1%" 
Also given in Table I are Sklar and Lyddane’s 
values for the Coulomb integral (aa;bb), for 
10<p< 20, and the exchange integral (ab;ab), for 
0<p< 20. Bartlett did not tabulate explicitly, 
but Sklar and Lyddane appear to have dug it out 
of his paper, allowing them to calculate (ab;ab) 
They found values for large p 
developing approximation formulas which are 
valid for large p and by extrapolating Bartlett’s 
tables for small quantities.” 

Bartlett failed to tabulate 7, and i19, so that 
the exchange integrals (ab;a,b,) and (a,b,;a,b,) 
are not included in Table I. 


HYBRID COULOMB-EXCHANGE INTEGRALS 


The hybrid Coulomb-exchange integrals are 
somewhat troublesome. The only ones which 


9M. Goeppert-Mayer and A. L. Sklar, J. Chem. Phys. 
6, 645 (1938). 

” The signs of is, is, and i; are wrong in reference 5; this 
was pointed out by W. H. Furry and J. H. Bartlett in 
Phys. Rev. 38, 211 (1932). 

i This procedure was first carried out for (aa;bb) and 
— by Goeppert-Mayer and Sklar (reference 9%) 
corrected by Sklar and Lyddane (reference 1)]. 


(31) 


| 
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will be considered here are 


(ab;bb) = faye) 
(ab;b,b,) = f (33) 


(ab,;bb,) = f a(1)b(2) 


These will be evaluated by the method which 
was applied by Sklar and Lyddane to the first 
of them—expansion of (2/ri2) in spherical har- 
monics followed by successive integrations over 
$1, ¥2, and the remaining coordinates of elec- 
tron 1. 

Substituting into 


(ab:bb) = f exp[ —(Z/2)(a1-+b1) 


Xsin@a; sin@o; (2/r12) 

Xexp(—Zb2)be? (34) 
the expansion of (2/r12) in spherical harmonics 
about nucleus B (Eq. (16) with a replaced by 3) 


and the elements of volume dv;= (dv;/d¢:)d¢i, 
= be? sinBbodbod one finds 


(ab;bb) =2(Z5/32m)? (n—|ml)! 


M=—n 


Xexp[im($1— $2) 


x 
x f 


xf exp — (Z/2) (ai +01) sin@a, 
Xsin6o,P,,'"! (35) 
By virtue of (21), (22), and the relations 
0 


(36) 
exp ( be) 2 
0 


=(1/Z*)g;, 3(Zb:), 


(35) becomes 
(ab;bb) = (2/24) (Z*°/32)? 
f 
XsinGa; sin; - { (4/3) gs, 3(Z01) 
— (1/15) P2?(cos61) |gz, 3(Zb1) } 
-(dv;/d¢1), (37) 


which is in agreement with Sklar and Lyddane’s 
Eq. (7). Similarly, making use of (28) and (29), 
one obtains 


(ab;b.b.) = (2/Z*)(Z°/32)? 
x 
XsinGa; sino; { (4/3) gs, 
—[(4/15) P2(cos@o;) 
+ (1/15) P2?(cos601) }gz, 3(Zb1) } 
(38) 


and 

(ab,;bb,) = (2/24) (Z°/32)? 
(2/2) ai Jan 
XsinOa; - { (1/15) 
X P2?(cos61) gz, 3(Zb1) } (dvi/dd1). (39) 


As independent coordinates of electron 1, it is 
convenient to take 


a=Za, and B=Zb,. (40) 


With p=ZR and (37)- 
(39) then become 


Z-*(ab;bb) = (1/7680p") —(1/2)(a+8)) 
X {5gs, 1(8) — [1+ (9/166"p) (a*+B*+ 
-adapdp, (41) 
=(1/7680p%) [expt —(1/2)(a+8)) 
 {5g5,1(8) (a! 
— (3/8B?p*) (a?B? + a*p? 3(8)} 
-adaBdB, (42) 
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and 
(ab,;bb.) = (2) (ab;bb) — Z— (ab (43) 


Equation (41) is equivalent to Sklar and 
Lyddane’s Eq. (8). 

According to Sklar and Lyddane, one must 
now indulge in a ‘‘considerable amount of 
rearrangement and integration by parts.” In 
terms of the integrals 


= (1/2) exp(o/2) 


Xa*B'dadp =I(t,s), 
J(s,t) = (3/2) exp(p/2) 


x f expl—(1/2)(a+36) 
K(s) = (1/2) exp(o/2) f 
x Lexp(8) 


A,=exp(2p) f exp(—2a)a'da 
=f (44) 
0 


B,=exp(—p) f exp(a)a‘da 
0 
f 
Cmexpo) exp(—a)a‘da 
=f 
0 
f 


= f 
0 


there result!” 


12288p* exp(p/2)Z—(ab;bb) 
de), 
12288p* =48> 1° — > 
+9 
12288 exp(p/2)Z—(ab,;bb,) 1°* — 
t+ Ls +L"), 


2 Details of the reduction of Eqs. (41)-(43) to Eqs. (45) 
will be furnished upon request. 


(45) 


L. CRAWFORD, JR. 
where the sums >> are given by: 


= (162 — 8p") (5,0) + (— 108 + 16p?) 
-I(3,2) + (12p?+ 16p*) (3,0) 
+ (27 —8p?)I(1,4) + ( — 60p?+ 
-I(1,2) + (66p*—8p°) (1,0), 
= (54 —8p*)I(5,0) + (—36+16p?) 
- (3,2) + (—60p?+ 16p*)I(3,0) 
+ (9—8p?)I(1,4) + (12?+ 169%) 
-I(1,2) ; 


(46) 


=3I(9,1) +18I(9,0) —127(7,3) 
—72J(7,2) —(288+12p?)J(7,1) 
— (864+72p?)J(7,0) +18(5,5) 
+1087(5,4) + (432+ 12 2) (5,3) 
+ (1296 +72?) J(5,2) +(2592+320p? 
+18p*)J(5,1) +(2592+-1120p? 
+108p*)J(5,0) —12J(3,7) —72J(3,6) 
— (288 — 12?) J(3,5) — (864 —72p?) 
- J (3,4) — (1728 — 2249? — 12p*) J(3,3) 
— (1728 —352p?—72p*)J(3,2) 
+ (192p2+ 224p4— 12) J(3,1) 
+ (192p?+352p!— 72p*) J (3,0) 
+3J(1,9) +18J(1,8) + (72 —12p?) 
J (1,7) + (216 —72p?)J (1,6) + (432 
—256p?+ 18") J(1,5) + (432 —608p? 
+108p*)J(1,4) — (960p? — 368p! 
+12p*)J(1,3) — (960p?— 784p!+72p°) 
- J(1,2) + (1056p4— 256p°+3p8)J(1,1) 
+ (47) 
¥2*=J(9,1) +6F(9,0) —4(7,3) —24J(7,2) 
—(96+4p?)J(7,1) — (288+ 24p?) 
-J(7,0) +6 (5,5) +36I(5,4) + (144 
+4p?)J(5,3) + (432+ 24p?)J(5,2) 
+ (864+ 128,?+- 6p) J(5,1) + (864 
+544p?+36p*)J(5,0) —4J(3,7) 
—24J (3,6) — (96 —4p?) J (3,5) — (288 
— 24%) J(3,4) — (576 — 32p?— 4p!) 
- (3,3) — (576+ 224p? — 24p') J(3,2) 
— (960p? — 32p!+-4p*) J(3,1) — (960p? 
+224p'+24p*) J(3,0) +J(1,9) 
+6J(1,8)+(24—4p?)J(1,7) 
+(72—24p?) J(1,6) + (144-649? 
+6p*)J(1,5) + (144 —32p?+ 36p*) 
-J(1,4) +(192p?+ 80p* — J(1,3) 
+ (192p?—80pt— 24°) J(1,2) 
— (672p*+-64p° — J(1,1) 
— (672p'+32p* —6p*) J (1,0) ; 


x 
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3= —3K(9)+(288+12p?)K(7) 
— (416p?+ 18p*) K (S) — — 12p°) 
-K(3)+(160p*—3p*)K(1), (48) 
—K(9)+(96+4p*)K (7) — (2245? 
+6p*)K (5S) + (160p*+4p°)K (3) 
— (32p%+p*)K(1) ; 


= (443/840) p® — (303/140) p® 

— (8276/105)p? —1086p° 

—7208p'— 27045? —69708p? 

— 1886895 /16)p —(791505/8), 

= (443/2520)®— (101/140) 

+ (676/105) — 

—1208p'—6615p*—21252p? 

— (628965 /16)p — (263835/8); 


(49) 


(264 —54p 

—96p?) Cs — (1800p — — 168 *) C; 
+ (4864p? — 336»? — 144p*) C, 
— (6656p? — 288p* —48p*) C; 

+ (4736p* —96p*) C2 — 1408p°Ci, 

(88 — 18p—32p?) 
— (600p — 64p? — 56p*) C; 
+ (1536p? — 112p?—48p*)C, 
— (1792p — 96p*— 16p*)C; 
+ (896p* —32p*) C,—128p°C;; 


(50) 


— (264+ 162p — 288p*)(A6+De) 
+ (1800p +576p? — 504p*) (A;+Ds) 
(4864p?+ 1008p? — 432p')(A 
+ (6656p? + — 144°) (A3+Ds) 
— (4736p*+ 288°) (A2+Dz) 
— (88+54p —96p?) 
+ (600p + 192p? — 168p*)(As+Ds) 
— (1536p? + 336p* — 144p*)(A4+Du) 
+ (1792p? + 288p' — 48p*)(A;+Ds) 
— (896p*+96p°) (A2+Dz2) 


=(1/2)(Le—- Le’). (52) 


The first of Eqs. (45) does not agree with 
Sklar and Lyddane’s Eq. (9) for Z—!(ab;bb). The 
following changes are necessary in the latter to 
bring it into accord with the former :" 


(51) 


3 The ).’s defined in the present paper are multiples of 
the ones employed by Sklar and Lyddane. 


(1) Change (1/2)>°; to —¥3. 


(2) In change p®J(1,1) — (p°—p*)J(1,0) to 
(1/16)p*J(1,1) —[(38/3)p*®— (3/8) p* JJ(1,0). 


(3) In change —(107/9)p*— (4195/8) 
— (14279/6)p?— (637065/128)p — (277815/64) to 
(901/9)p*+ (3005/8) *+ (5809/6)? 

+ (209655 /128)p+(87945/64). 


(4) In change (14/3)p*(A4+D,) to 
14p*(A4+D,). 


Some of these errors are doubtless misprints." 

Values of Z-(ab;bb), and 
Z—(ab,;bb,), calculated from (45) and the tables 
of Sklar and Lyddane for the IJ(s,t), J(s,t), and 
K(s)," are given in Table I. These integrals fall 
nicely between the corresponding Coulomb and 
exchange integrals, as they should. 


PENETRATION INTEGRALS 
The Coulomb penetration integral 


(A :bb) = — 
[A=carbon], (53) 


and the exchange penetration integral 


f 
[A=carbon], (54) 


have been evaluated by Sklar and Lyddane [in 
Sklar and Lyddane’s notation, (A:bb) =Q22(p), 
(A :ab) =Q12(p) ], but the numerical values quoted 
by those authors do not agree with the formulas 
given by them. 

If the 1s electrons on carbon A are treated as 
if on the nucleus,!® the potential H,*(1) on 
electron 1 due to the carbon atom A is (in c.g.s. 


“In a letter of one of the authors (reference 6), Dr. 
Lyddane states: “The errors . . . in (ab;bb)/Z ... are 
all typographical with one exception; the value of >). 
is wrong in our paper.” 

16 For p=14.5, [(5,0) =2.4864 10°, not 2.486410" as 
Sklar and Lyddane state. 

16 This method of evaluation was outlined by M. 
Goeppert-Mayer and A. L. Sklar (reference 9) (see also 
Sklar and Lyddane, reference 2). 
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units) 
H,*(1) = — (4e2/a:) +(1/ay) f éo(a)a%da 
0 


+ J (55) 


where a; is the distance between electron 1 and 
nucleus A, and 


o(a) | Poo |?+ | + | 
= fC (56) 
+ | ](dv/arda). 
But 
f o(a)atda= f 


+ | (57) 


so that in c.g.s. units 
(a—ao(a)ada (58) 
a1 
and in atomic units* 


H,*(1) = —(2/a,) f (59) 
Here? 
o(a) = (2*/4) (60) 
Hence (in atomic units) 


Ha*(1) = —(1/2a1) exp(—Za;)[16+12(Za1) 
+4(Za;)?+(Za)*]. (61) 
Substitution into (53) of (61) and the first of 


Eqs. (2) and integration by the method applied 
above to the overlap integrals yield the formula 


Z-(A:bb) = (1/1920) exp(—p)[(4/7)p° 
+ (67/7) p*+86p? +-461p?+1125(p+1) ] 
[A=carbon ]. (62) 


L. CRAWFORD, JR. 


Similarly, integration of (54) by the method 
applied above to the hybrid Coulomb-exchange 
integrals gives 


:ab) = (1/384p*) exp(—p/2) 
X { —J(5,3) —4J(5,2) —12J(5,1) —16J(5,0) 
+2/(3,5)+8J(3,4) +24J(3,3) +32J(3,2) 
—J(1,7) —4J(1,6) —12J(1,5) —16J(1,4) 
+p°(2J(3,3) +8J(3,2) +24J(3,1) +327(3,0) 
+2J(1,5)+8/J(1,4) +24J(1,3) +32J(1,2) ] 
+16/J(1,0)]} [A=carbon]. (63) 


For p=0, each of these integrals reduces to 


Z—(A = 1125/1920 =0.5859 
[A =carbon]. (64) 


Equation (62) differs from Sklar and Lyd- 
dane’s Eq. (18) by a factor of 2 and the appear- 
ance of 86p* in place of 46p*; Eq. (63) differs 
from Sklar and Lyddane’s Eq. (19) by a factor 
of 2. Equations (62) and (63) reproduce exactly 
Sklar and Lyddane’s numerical values, however, 
indicating that those authors employed the 
correct equations in their computations.’7—” 
These values are included in Table I. 

Similar methods could be applied to the 
calculation of (A:b,b,) and (A:a,b,). 
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17 Dr. Lyddane has confirmed this. 

18 V. Griffing, J. Chem. Phys. 15, 421 (1947), noted that 
the values given in the first row of Sklar and Lyddane’s 
Table VI are larger than those calculated from their 
Eq. (18), but laid the blame on the calculations rather 
than on the equation. 

19 Other misprints in the Sklar-Lyddane paper include. 

In Eq. (6) for Z read 25, for bz read b:*. 

In Eq. (8), for (5/8) read (5/8). 

In Eq. (15), for (1+e7?) read (1+e*?). 

In Eq. (16), for (1/2)sA,—1 read (1/2)sA 5-1. 
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The thermal dissociation of hydrogen on hot tungsten filaments has been used in the inves- 
tigation of the reaction of atomic hydrogen with acetylene. A quantitative method has been 
developed for the detection of atomic hydrogen at pressures of the order of 10-* mm. The 
catalytic effect of acetylene on the recombination of atomic hydrogen has been studied and the 
rate constant of the rate-determining step calculated to be 8.6X10-“ molecules“ cm sec.~! 
at 25°C. In contrast with results obtained by the discharge tube method, small amounts of 


hydrogenated products are formed. 


INTRODUCTION 


HE reaction of atomic hydrogen with 

acetylene has been studied by a number 
of investigators using the discharge tube method 
of Wood and Bonhoeffer.'-> The nature of the 
discharge is such that a flow system must be used 
in which the total pressure is of the order of 0.5 
mm. In most experiments the pressure of atomic 
hydrogen was from 0.1 to 0.2 mm, and under 
these conditions it was found that while small 
amounts of acetylene greatly accelerated the 
consumption of hydrogen atoms, no products 
other than hydrogen and acetylene were de- 
tected. Similar results were obtained with atomic 
deuterium,‘ but in this case there was a rapid 
exchange. The collision yield of deuterization 
was estimated to correspond to an activation 
energy of less than 5 kcal. per mole and two 
possible mechanisms were suggested to explain 
the results: 


H+C.H:=C2Hs, (1) 
(2) 
(3) 
H+C.H (4) 


*Present address, Division of Chemistry, National 
Research Laboratories, Ottawa, Ontario. 

1K. F. Bonhoeffer and P. Harteck, Zeits. f. physik. 
Chemie 139, 64 (1928). 

2H. v. Wartenberg and G. Schultze, Zeits. f. physik. 
Chemie B2, 1 (1929). 

3K. H. Geib and P. Harteck, Ber. 66, 1315 (1933). 

‘K. H. Geib and E. W. R. Steacie, Zeits. f. physik. 
Chemie B29, 215 (1935); Trans. Roy. Soc. (Canada) III, 
29, 91 (1935). 

5D. J. Le Roy and E. W. R. Steacie, J. Chem. Phys. 12, 
369 (1944), 


Using \the method of photosensitization to 
produce atomic hydrogen, Le Roy and Steacie*® 
found as products, ethane, ethylene, butane, and 
a polymer of the approximate composition 
(C.H,),. Since the pressure of atomic hydrogen 
was lower than that obtained by the discharge 
tube method by several orders of magnitude, 
they concluded that hydrogenated products 
could be produced only in the presence of large 
H::H ratios. On the basis of Mechanism I vinyl 
radicals could then partake in the reaction 


(5) 
or 
H2=C2Hs. (6) 


If Mechanism II obtains the corresponding reac- 
tion would be 


CH+H2=CHs, (7) 


followed by either (5) or (6). Vinyl radicals would 
play an important part in either case. 

In order to obtain a clearer insight into the 
elementary reactions involved it seemed desirable 
to develop a method capable of yielding pressures 
of atomic hydrogen which would overlap those 
previously used. In addition it was essential to 
devise some means of measuring these pressures. 


EXPERIMENTAL DETAILS 


Despite Langmuir’s classical work on the 
formation of atomic hydrogen by thermal dis- 
sociation on tungsten filaments this method has 
only recently been applied to chemical inves- 
tigations.* Relatively wide variations in condi- 


6J. P. W. Houtman, K. Van Den Berg and P. M. 
Heertjes, Rec. trav. Chim. 62, 672 (1943). 
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T 
S A 
} 
s 
Fic. 1. Apparatus for the production and detection of 


atomic hydrogen. 


tions are made possible by controlling the tem- 


perature of the filament and the pressure and 


flow rate of the hydrogen. 

The part of the system used to produce and 
detect atomic hydrogen is shown in Fig. 1. The 
tungsten filament, 7, 0.004 inch in diameter and 
from one and a half to two inches long, was spot- 
welded to conventional nickel-tungsten lead-in 
wires. The filament assembly was sealed to a 


standard 24/40 ground joint lubricated with 


Apiezon N grease; the filament chamber and 
ground joint were surrounded by a water jacket. 
The current was controlled manually and the 
filament temperature measured to +50°C by 
means of an optical pyrometer. Purified hydrogen 
entered at Hz, and the mixture of atomic and 
molecular hydrogen then passed into the reaction 
tube, R, through a jet approximately 8 mm in 
diameter. The opening in the acetylene jet, A, 
was 40 mm down stream from: the hydrogen jet. 
The reaction tube itself was 18 mm in diameter 
and about 25 cm long. 

The detector, D, was mounted on a standard 
24/40 ground joint, and could be placed at either 
one of two positions with the use of an extension 
joint. The detector wire was of platinum, from 
0.0020 to 0.0025 inch in diameter and approxi- 
mately 25 mm in length. It was spot-welded to 
nickel-tungsten lead-in wires. The nickel was 
shielded from contact with hydrogen by soft 
glass skirts, S, which reached down to the 
tungsten-nonex seal at the one end and at the 
other end were sealed to the fine platinum wire 
at the nickel-platinum weld. The small portion of 
the platinum wire covered by the soft glass did 


not introduce a serious error. All glass surfaces 
exposed to hydrogen, from the top of the tung- 
sten filament chamber to somewhat past the 
detector, were given a liberal coating of meta- 
phosphoric acid to reduce recombination on the 
walls. 

The detector wire, its leads, and a calibrated 
500 ma meter formed one arm of a Wheatstone 
bridge ; the other three arms consisted of a 1000 
ohm Leeds and Northrup standard resistance and 
two decade boxes. A G.M. table galvanometer 
with a sensitivity of approximately 0.4 micro 
ampere per mm served as the null detector. The 
bridge current was supplied by a storage battery 
and could be varied by means of a rheostat. 

Considerable time was devoted to a study of 
the characteristics of the detector before doing 
any experiments with acetylene. With the tung- 
sten filament off, but with hydrogen flowing 
through the system at the desired pressure, the 
bridge current was adjusted until the resistance 
of the detector wire corresponded to the desired 
temperature, say 400°C. Allowance was made for 
the resistance of the meter and leads. When the 
tungsten filament was turned on, the tempera- 
ture of the detector wire would increase due to 
the recombination of atomic hydrogen on its 
surface. The bridge current was then lowered 
until the bridge was again in balance. With the 
detector wire again at 400°C, heat losses thereby 
being eliminated, the decrease in the wattage 
supplied to the detector is a measure of the 
number of hydrogen atoms recombining on its 
surface per second. 

If the recombination coefficient is assumed to 
be unity’ the number of hydrogen atoms striking 
a unit area per second can be calculated in terms 
of the wattage decrease, Wy. It is equal to 


26.023 Wr 
4.183 X103,680XA 


where A is the area of the wire in cm?. The area 
of the wire was determined at the end of each 
series of experiments by cutting it off closely, 
measuring its length, and weighing it on a micro 
balance. From the kinetic theory expression for 
the number of particles striking a unit area per 
second it was then possible to calculate the 


7W. V. Smith, J. Chem. Phys. 11, 110 (1943). 
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pressure of atomic hydrogen, Py. The resulting 
equation is: 


Py=1483X mm. (a) 


It was assumed that the average velocity of the 
hydrogen atoms striking the detector wire was 
determined by the temperature of the walls of 
the reaction tube. Experiments using nitrogen 
instead of hydrogen showed that no heat was 
carried from the tungsten filament to the 
detector. 

The commercial electrolytic hydrogen used in 
this work was purified continuously by passing 
it over hot platinized asbestos and then through 
a trap cooled with liquid air, before admitting it 
to the reaction system. Acetylene was freed from 
acetone and other impurities by passing it 
through two traps cooled with a mixture of dry 
ice and acetone and then subjecting it to a series 
of bulb-to-bulb distillations. The middle fraction 
was stored in a two liter bulb. No impurities 
could be detected by analysis. 

In carrying out an experiment using acetylene 
a sample was transferred from the storage bulb 
to a calibrated gas burette. The volume of the 
gas in the burette was adjusted by a mercury 
piston to give a pressure somewhat less than 
atmospheric. This pressure was read on a U tube 
manometer having an electrical contact in the 
vacuum arm. As the acetylene left the burette 
and entered the reaction vessel through an 
orifice the electrical contact in the manometer 
operated a relay which caused a current to pass 
through a solution of sodium sulphate in an elec- 
trolysis cell. The hydrogen and oxygen produced 
raised the mercury piston and kept the pressure 
in the burette constant to within a fraction of a 
mm. The gases issuing from the reaction tube 
passed through several traps cooled with liquid 
air; the liquid air around the last trap was 
pumped continuously, reducing its temperature 
to a point where all gases except methane and 
hydrogen would be frozen out. It would seem 
highly unlikely that any methane would be 
produced in this reaction. The hydrogen was 
pumped off through a mercury diffusion pump 
and a Pressovac pump; it could be either col- 
lected and measured over a pneumatic trough or 
passed into a fume duct. The steady-state 
pressure in the system was adjusted by two 


Hoke high vacuum valves (No. 318), one on each 
side of the reactor. 

At the end of an experiment all residual 
hydrogen was pumped off and the traps warmed 
to room temperature. The gases were condensed 
into a bulb at the bottom of a second gas burette. 
After isolating this burette from the rest of the 
system, the liquid air was removed and the gas 
compressed into the calibrated portion by raising 
a mercury levelling bulb. After measuring its 
volume and pressure it was transferred, by means 
of a Toepler pump, to one of the small test tubes 
used with the Blacet-Leighton gas analysis 
apparatus. In later experiments a large part of 
the unreacted acetylene was first removed by a 
preliminary treatment with alkaline mercuric 
iodide under vacuum conditions. In the actual 
analyses acetylene was removed with alkaline 
mercuric cyanide beads.* Ethylene was deter- 
mined by the hydroxy mercurial method,’ par- 
affins by combustion.*® 


EXPERIMENTAL RESULTS 


Production and Determination of Atomic 
Hydrogen 


The effect of the temperature of the detector 
wire on the wattage developed by recombination 
was investigated over a temperature range of 
approximately 125°C to 460°C, for various 
pressures of hydrogen. The whole range could 
not be covered at a single hydrogen pressure, 
since to heat the detector to 460°C under a high 
pressure of hydrogen would require a current in 
excess of the capacity of the bridge. A large 
number of experiments were made, of which the 
two recorded in Table I are typical. The re- 
sistance of the detector wire, R, is given rather 
than its temperature. Correction was made for 
the resistance of the leads and the meter. The 
limiting value of R at low currents was 1.33 
ohms. Wy was calculated from the expression 


Wu=R(I0?— In’), (b) 


where J) and Jy are the detector currents in the 
absence and in the presence of atomic hydrogen. 


8 F. E. Blacet, A. L. Sellers and W. J. Blaedel, Ind. Eng. 
Chem. (Anal. Ed.) 12, 356 (1940). 

®R. Pyke, A. Kahn and D. J. Le Roy, Ind. Eng. Chem. 
(Anal. Ed.) 19, 65 (1947). 
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It is clear from these experiments that Wg is 
independent of the temperature of the detector 
wire over quite a wide range; this is not proof, 
however, that Wy is a true measure of the pres- 
sure of atomic hydrogen. Experiments in which 
purified nitrogen was used instead of hydrogen 
showed no change in R when the tungsten fila- 
ment was turned on, showing that the effect was 
not due to hot gas. Two additional experimental 
tests were applied, one to study the effect of the 
temperature of the tungsten filament on Wu, the 
other to study the effect of hydrogen pressure. 

If it can be assumed that thermodynamic 
equilibrium is obtained at the tungsten filament, 
and that this equilibrium is “frozen” until the 
gas reaches the detector, then the decrease in 
standard free energy for the reaction H,=2H is 
given by 


~AF° = RT InP (c) 


where Py is the pressure of atomic hydrogen at 
the detector, PH: is the pressure of molecular 
hydrogen, R is the gas constant and T is the ab- 
solute temperature of the tungsten filament. Pu: 
is so close to the total pressure that the difference 
can be neglected. If the total pressure is kept 
constant, and if the standard heat of formation 
of atomic hydrogen is approximately independent 


_ of temperature over the range, there should be a 


linear relationship between the logarithm of Py 
and the reciprocal of the absolute temperature 
of the tungsten filament; the slope should be 


given by the expression: 


d logPx/d(1/T) =d logWu/d(1/T) 
= —AH°/(4.606R). (d) 


The results of two sets of experiments, done 
several months apart, are shown in Fig. 2. Dif- 
ferent detector wires and hydrogen pressures were 
used in the two cases and hence the intercepts of 
the two curves would not be expected to be the 
same. The values of AH® obtained from the 
lower and upper curves are 76 and 89 kcal., re- 
spectively, the latter value being the more 
recent. 

A more definitive test was made by computing 
the equilibrium constant, K,, for the reaction 
H.=2H at a fixed tungsten temperature of 
1700°C. The data are given in Table II. Giauque’” 
obtained a value of 3X10-* for K, at 727°C. In 
view of the uncertainty in our value of the tem- 
perature the agreement is probably as good as 
can be expected. 


The Catalytic Effect of Acetylene on the Recom- 
bination of Atomic Hydrogen 


-The effect on the concentration of atomic 
hydrogen of the addition of small quantities of 
acetylene is shown in Table III. The total 
pressure in the reaction tube, essentially the 
pressure of molecular hydrogen, is given in 
column one. The acetylene pressures in column 
three were calculated from the measured flow 
rate of hydrogen and the rate at which the 
acetylene left the first burette. The reaction time 


TABLE I. Effect of detector temperature on Wy.* 


Hydrogen pressure =0.547 mm. 
Flow rate** =25 cc/min. 
=75 mm. 
H 


Hydrogen pressure =8.56 mm. 
Flow rate—330 cc/min. 
H 


‘To Wu 
amps. amps. ohms watts 
0.269 0.213 3.237 0.0875 
0.259 0.200 3.096 0.0839 
0.249 0.180 2.894 0.0857 
0.239 0.154 2.672 0.0892 
0.229 0.134 2.404 0.0828 
0.220 0.099 2.186 0.0844 
0.215 0.091 2.112 0.0800 
0.210 0.068 2.039 0.0805 

0.0843+3% 


Io Tu 
amps. amps. ohms watts 
0.498 0.337 2.627 0.354 
0.477 0.296 2.421 0.339 
0.455 0.228 2.259 0.350 
0.434 0.168 2.146 0.344 
0.424 0.129 2.093 0.341 
0.415 0.042 2.020 0.345 
0.405 0.000 1.996 0.336 
0.344+1.5% 


* Tungsten filament temperature =1650+50°C; R =1.33 ohms at 25°C. 
** Flow rates are exp: 


10 W. F. Giauque, J. Am. Chem. Soc. 52, 4808 (1930). 


as cc of hydrogen, measured at 25°C and one atmosphere pressure, per minute. 
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given in the fourth column is the time required 
for the reactants to pass from the acetylene jet 
to the detector. At the beginning of an experi- 
ment the apparatus would be run for some time 
before adding the acetylene. This served to 
establish a constant flow rate of hydrogen and to 
“burn up” any impurities on the detector. The 
final steady value of the wattage due to recom- 
bination on the detector is given in column five. 
The values found after adding acetylene to the 
gas stream are given in column six. Pseudo-first 
order rate constants were calculated from the 
expression 


In(Wu®/Wu) =In((H]o/[H]) (e) 


The values of k are given in the last column of 
Table III. Despite the fluctuations in k it is 
clearly independent of the pressure of molecular 
hydrogen. It would appear that the rate equation 
for the removal of atomic hydrogen is of the form 


d(H Vdt= (f) 


It should be noted, however, that the range of 
acetylene pressures used was not large. 


The Formation of Products 


All attempts to detect products other than 
hydrogen and acetylene under the conditions of 
Table III were unsuccessful. However, con- 
densible products were found when the acetylene 
pressure was increased to a value sufficient to 
reduce the pressure of atomic hydrogen at the 
detector to zero. Reaction times of as much as 
0.28 seconds were used but in no case was the 
percentage of products greater than about 34 
percent, based on acetylene. Increasing the tem- 
perature of the-reaction tube to 100°C had little 
effect. Although the analyses were subject to 
considerable experimental error, evidence was 
found for the presence of both an olefin and a 
paraffin. The latter appeared to be ethane. Ex- 
periments performed with the tungsten filament 
at room temperature showed no measurable 
amount of impurity in the original gases. In 
general the amount of acetylene recovered 
agreed with the amount introduced within two 
percent, 


Ne 


46 5-4 6-2 70 


+70" 


Fic. 2. The effect of the temperature of the tungsten 
filament on the logarithm of the wattage produced by 
recombination of the atomic hydrogen on the platinum 
detector. The upper and lower curves correspond to heats 
of recombination of 89 and 76 kcal., respectively. 


DISCUSSION 


The results given in Table III can be inter- 
preted on the basis of either Mechanism I or 
Mechanism II. It has not been pointed out 
previously that reactions (1) and (2) alone, or 
(3) and (4) alone, do not fulfill the requirement 
that the ultimate products be stable. For ex- 
ample, if (3) and (4) were the only reactions 
occuring the concentration of ethynyl radicals 
would build up to a maximum and remain at that 
value indefinitely. The maximum value would 
be given by 


d[C3H ]/dt = H ] = 0, 
or 
(C2H (g) 


Some reaction consuming ethyny] radicals would 
have to occur even if no products could be 
detected. This might be reaction (7) or (8), 


2C2H = (8) 


The analogous reactions for Mechanism I would 
be (5), (6), (9) or (10), 


2C.H3 = C,4Hg, (9) 
2C;H;= (10) 


Some evidence for (10) has been found in another 
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TABLE II. Apparent oe constant for the reaction 


Pres- : 
sure Wu Pu Kp Flow rate 
mm watts mm X108 mm X106 cc/min. 
1.15 0.149 4.62 18.5 40 
5.64 .194 6.02 6.41 200 
8.75 .234 7.25 6.02 360 
12.9 .283 8.77 5.97 460 


* Tungsten filament temperature =1700+50°C; distance, tungsten 
filament to detector =73 mm. : 


TABLE III. Recombination of atomic hydrogen in the 
presence of acetylene.* 


tance 
Pres- jet-de- CsHe k 
sure tector mm T WH? Wu molec.“ 
mm mm X108 sec. watts watts cm sec.~! 


1.00 52 1.41 0.0246 0.0727 0.0620 1.42x10—% 
3.62 52 1.47 0.0250 0.0732 0.0539 2.68 
6.14 52 1.61 0.0231 0.0680 0.0507 2 
10.80 52 1.35 0.0220 0.0641 0.0553 1 
1.98 52 1.68 0.0307 0.0953 0.0789 1. 
1.02 110 1.86 0.0447 0.0916 0.0560 1. 
4.59 °110 2.38 0.0555 0.0675 0.0394 1 
9.36 110 2.10 0.0519 0.0647 0.0368 1 
14.60 110 1.81 0.0496 0.0620 0.0403 1 


* Reaction tube at room temperature; tungsten filament temperature 
=1650+50°C; cross section of reaction tube =2.26 cm?. 


connection." To avoid confusion it will be 
assumed that Mechanism II obtains; any con- 
clusions drawn can be expressed in terms of 
Mechanism I if the need should arise. 

Unless atomic hydrogen is regenerated in some 


| way its rate of decay will be determined by reac- 


tions (3) and (4), and at any time, ¢, its concen- 
tration will lie within the limits, 


Jo |t >In [H] > ]o 
}t.  (h) 


If (3) were faster than (4) the upper limit of (h) 
would be approached, the concentration of 
ethyny]l radicals would become quite appreciable, 


J. W. Hodgins, A. W. Tickner, and D. J. Le Roy, in 
press. 


and a considerable amount of product would be 
expected. For this reason we believe that the rate 
constant, k, of Table III is probably very close 
to 2k3(2k; for Mechanism I). 

It is concluded that at 25°C the rate constant 
of the rate-determining step for the recombina- 
tion of atomic hydrogen catalyzed by acetylene 
is 4k, or 8.6X10-" molecules“ cm sec.—!. It has 
been customary in work of this type to express 
experimental results in terms of an activation 
energy, even when the data are confined to a 
single temperature. Although such a procedure is 
admittedly dangerous, the values obtained may 
be of some use in comparing data from various 
sources. The usual value assigned to the steric 
factor is 0.1. An estimate of 2.72A for the mean 
collision diameter was obtained from Sponer’s 
values for the C—C and C—H bond distances in 
acetylene” and Bonhoeffer and Harteck’s value 
of 2.14A for the diameter of the H atom." The 
value of the activation energy of reaction (3) 
(or (1)) computed in this way is 3.9 kcal. per 
mole. 

Considerable work remains to be done on this 
reaction. Not only is it impossible to distinguish 
between the two mechanisms, but the range of 
experimental conditions which we have been able 
to use does not eliminate the possibility that the 
catalytic recombination may be at least partially 
heterogeneous. It should also be noted that the 
concentration of atomic hydrogen was deter- 
mined at only one point in the reaction vessel 
during an experiment and hence no correction 
could be made for possible recombination on the 
wall, even in the absence of acetylene. 

The authors wish to express their thanks to the 
Advisory Committee on Scientific Research of 
the University of Toronto for a grant which sup- 
ported this research. 

2 A, Farkas and H. W. Melville, Experimental Methods 
3 S. Reactions (MacMillan and Company, Ltd., London, 


18K. F, Bonhoeffer and P. Harteck, Photochemie (Stein- 
kopff, Leipzig, 1933). 
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It is shown that if a crystal is approximated as a har- 
monic oscillator, only a very small number of vibrational 
modes may be active in the infra-red or Raman spectra. 
These modes are all totally symmetric with respect to 
translation. They must also have proper symmetry with re- 
spect to the factor group of the space group which we have 
called the unit cell group. The active frequencies may be 
calculated from a single unit cell. It is shown further that 
the number and symmetry type of vibrational modes may 
be correctly calculated from the local symmetry (site 
group) of the center of gravity of each kind of molecule. 
In the harmonic approximation the spectrum should 
consist of very sharp lines. The wide bands obtained 
experimentally are probably due to combination lines with 
lattice vibrations and indicate anharmonic coupling be- 
tween molecular and lattice modes. Such combinations 
may Cause apparent violations of selection rules for internal 


vibrations. The symmetry basis for obtaining selection 
rules in this case is given. The importance of obtaining 
spectra over a wide temperature range and particularly at 
very low temperatures is emphasized. In cases where all 
internal modes have a frequency great compared to the 
frequency of lattice modes the problem is examined from 
the point of view of the Born-Oppenheimer approximation. 
It is shown that for some symmetries perturbations linear 
in the lattice displacement may occur. In these cases 
considerable splitting of degenerate molecular vibration 
modes is possible. Groups for which first-order splitting is 
not expected are listed. Symmetry species of other groups 
are classified into those which may and which may not be 
split in this approximation, and the lattice modes re- 
sponsible for splitting are classified as rotations or trans- 
lations and by symmetry type. 


I. INTRODUCTION 
HE molecules of a crystal, unlike those of 


a gas, have a relatively fixed orientation” 


in space. If the structure of the crystal is known 
from x-ray data, the use of polarized radiation 
in infra-red or Raman investigations aids enor- 
mously in making vibrational assignments. On 
the other hand, if the vibrational assignments of 
the free molecule are known, it may be possible 
to use the spectroscopic data to indicate crystal 
structure, particularly in hydrogen containing 
compounds for which it is exceedingly difficult 
to obtain x-ray data. It may also be possible to 
obtain data regarding the crystalline field. 
Unfortunately, the spectra obtained from 
molecular crystals are complicated and their 
interpretation has not been clear. It has been 
known for a long time that the vibration fre- 
quencies of molecules in crystals are not very 
different from the corresponding frequencies in 
the gaseous state. When the problem is examined 
in greater detail, however, a number of significant 
differences appear. There is a great change in 
the shape of the envelope or the nature of the 


* The work reported here was supported in part by the 
Office of Naval Research under Contract N6ori-88, Task 1. 


fine structure of vibration bands. Degenerate 
vibrations are split, even in cases such as calcite, 
where the equilibrium crystal symmetry as well 
as the ion symmetry is consistent with degener- 
acy. Both degenerate vibrations are split by 
about 50 cm! in this case. Finally, selection 
rules are altered and apparently ‘‘forbidden’’ 
transitions appear. 

This problem has been discussed by a number 
of authors, but it seems desirable to develop the 
problem of the dynamics of molecules in crystals 
systematically, with particular emphasis on their 
rigorous basis in symmetry. We shall give a 
quite general discussion of the crystal in the 
harmonic oscillator approximation, noting that 
the predicted spectra are very simple. We shall 
conclude that the harmonic approximation is 
inadequate to account for important experi- 
mental effects. Since the coupling between mole- 
cules in a crystal is usually weak compared to 
the intramolecular forces, we shall introduce a 
general crystalline field as a perturbation and 
examine the more important consequences. These 
arise chiefly through anharmonic coupling be- 
tween lattice modes of the crystal and internal 
modes of vibration of the molecules, 
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Il. SYMMETRY PROPERTIES 


Let us consider the symmetry properties of a 
crystal which we shall assume to be infinite. 
The theory of space groups? shows that a lattice 
structure must have the symmetry of one of the 
230 space groups. The space group of a crystal 
is that group (in a mathematical sense) of 
operations which carries each atom into an 
identical atom. These symmetry operations in- 
clude the translations which generate the lattice 
plus the symmetry operations familiar in point 
groups. The unit cell is the smallest unit in which 
no atoms are equivalent under simple transla- 
tions. However, some of them may, in general, 
be equivalent under any of the other symmetry 
operations. The crystal lattice is constructed 
from the unit cell by the translations which 
carry any unit cell into any other. Similarly, the 
complete set of symmetry elements of the crystal 
is obtained from those of the unit cell by these 
same translations. The choice of unit cell is not 
unique but this makes no difference in the 
following discussion. 

The space group is then an infinite group. 
However, if we define the translations which 
carry a point in a unit cell into the equivalent 
point in another cell as identity, we define a 
finite group which is a factor group of the space 
group.” With our convention regarding transla- 
tions, this group describes the symmetry of a 
unit cell, and we shall call it the wmit cell group. 
In general, it is mot a point group although the 
possible factor groups are isomorphous with the 
32 point groups which may occur in crystal- 
lography. 

If we define the operations of the unit cell 
group as (U), those of the infinite group of 
translations by (7), and those of the space 
group by (S$), we can then write: 


(S) =(U) X(7). (1) 


It is a theorem of group theory that irreducible 
representations of a factor group, (U), are also 
irreducible representations of the entire group, 
(S). In these particular irreducible representa- 


1A. Schoenflies, Krystallsysteme und Krystallstruktur 
(B. F. Teubner, Leipzig, 1891). 

?W. H. Zachariasen, Theory of X-Ray Diffraction in 
Crystals (John Wiley and Sons, Inc., New York, 1945). 
The theory of space groups is concisely presented in 
Chapter 2, 


tions the elements of the invariant subgroup, 
(T), are all given by identity. Furthermore, these 
are the only irreducible representations having 
this property. Consequently, the set of irreduc- 
ible representations of (.S) which are totally sym- 
metric with respect to translation are given by 


(2) 


No simple structure can be assigned to the in- 
finity of remaining irreducible representations, 
but in the detail to which we shall go we will 
find them unnecessary.** 

It is possible to simplify still further. About 
any point in the crystal, there is a local symmetry 
which, for most points, will consist only of the 
identity operation, E. If, however, the point in 
question is located on one or more symmetry 
elements, the corresponding operations leave 
that point invariant. The symmetry operations 
belonging to the elements passing through the 
point form a point group which describes the 
symmetry of the crystal viewed from the point 
in question. Following the nomenclature of 
Halford’ we shall call this group which describes 
the local symmetry a site group. A point of 
given local symmetry may be invariant or have 
one, two, or three degrees of freedom, depending 
on whether it lies at the intersection of two or 
more elements (not both planes), on an axis, on 
a plane, or on no symmetry elements. Since the 
remaining symmetry operations of the space 
group (or unit cell group) carry the point into an 
equivalent point (which has, therefore, the same 
local symmetry about it), there is a definite 
number of equivalent points of each type. A 
convenient table giving all of the possible local 
symmetries (site groups) of each type for any 
space group, together with the number of 
equivalent points of each symmetry, has been 
given by P. Nigglit and more recently by Hal- 
ford* who does not, however, distinguish between 
points of the various types. These site groups 
are the subgroups of the unit cell group or, 


** Note added in proof. The author has been informed by 
Mr. H. Winston and Professor R. S. Halford that they will 
discuss the irreducible representations of space groups in 
detail in a forthcoming publication. 

3R. S. Halford, J. Chem. Phys. 14, 8 (1946). 

4P. Niggli, Geometrische Kristallographie des Diskon- 
tinuums (Gebriider Bornstraeger, Leipzig, 1919), Table IV, 
pp. 404-411, 
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alternatively, factor groups of the subgroups of 
the space group. 

Consequently, all of the operations of the unit 
cell group may be obtained by taking all possible 
products of operations contained in the site 
groups, noting, of course, that any operation 
which carries a point into the equivalent point 
in another unit cell must be treated as identity. 
It is always possible to write linear relations 
between the irreducible representations of the 
unit cell group and those of any site group. 
Irreducible representations which are so related 


may be called corresponding irreducible repre- | 


sentations. The correspondence may be deter- 
mined in any specific case by comparing the 
character tables of the site group and the unit 
cell group. (The character table of the unit cell 
group can always be obtained from that of the 
isomorphous point group containing the same 
subgroups.) In this way, correlation charts can 
be set up between the various site groups and 
the unit cell group. These relations are discussed 
analytically in Section III. Bethe® has described 
the correlation between spherical symmetry and 
the local symmetries which may occur in crystals, 
and Mulliken® has discussed the correlations 
between the point groups in some detail. 

To illustrate these ideas, let us examine the 
space group D3,°, which describes the crystals of 
calcite and a-Al,O;3. From Niggli we find that 
the unit cell contains two equivalent points of 
symmetry D3, two of symmetry C3;(.S¢), and six 
of symmetry C,(S2). These points are fixed 
uniquely. In addition, it contains an infinite 
number of sets of four equivalent points of 
symmetry C3, which are arbitrary in that one 
member of each set may be selected at any point 
on the C; axis (i.e., it has one degree of freedom), 
and sets of six equivalent points of symmetry C2, 
one of which is any point on a C, axis. Finally, 


TABLE I. 
E 2C3i(2S6) 2€3 3C2*t 3C2 
Aw 1 1 1 1 1 1 
Ai 1 —1 1 1 
Ag 1 1 1 1 
1 —1 1 —1 
E, Z -1 —1 2 0 0 
Ey 2 1 —1 —2 0 0 


5H. Bethe, Ann. d. Physik 3, 133 (1929). 
®R. S. Mulliken, Phys. Rev. 43, 279 (1933). 


it contains sets of twelve points about which 
there is no symmetry at all. The character tables 
for these possible local symmetry groups are as 
follows:. 


D3\|E Cs| E C3 C3? Se Seo 
1-1 1-1 -1 
-1 0 E 1 w wo 1 w w? 
7 w w 1 w 

1 w w —1 —w 

Eu 1 w —w —w 

C:|\EC: CF |Z i C:|E Cs 
A,|1 1 1 
E 1 w wo Bii -1 

1 w 


In these tables w =exp(27i/3). 

The possible elements of D3a® are obtained by 
taking all of the possible products . between 
elements of the subgroups, noting, of course, that 
the points about which the point symmetries are 
taken are not coincident in space, and that 
lattice translation is regarded as identity. It is 
not necessary to consider any groups which are 
subgroups of other groups, so that only D; and 
C3; need be considered. In this case, the only 
unusual operation is C2:7 where the axis and 
center of symmetry do not coincide. In crystal- 
lography this element is called a glide plane and 
the corresponding operation o, (glide). Table I 
the character table for the unit cell group of D3a°®. 

The group D3¢*° has the same character table 
and is isomorphous with the point group D3u. 
The correlations between the irreducible repre- 
sentations of the various site groups and the 
unit cell group are given in Table II. (When 
two lines originate at the same £ it is only 
pseudodegenerate and is composed of two non- 
degenerate irreducible representations whose 
characters form complex conjugate pairs.) 

A few very general statements about the trans- 
formation of certain quantities by a translation 
will suffice for our purposes. Since translation 
leaves any vector invariant, it is necessarily 
true that, in addition to vectors, dot and cross 
products of vectors as well as rotations are in- 
variant. Furthermore, all orders of derivatives 
are invariant under translation. 
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Now we shall be lademncnnd in matrix elements 
of the form 


where M is any operator and, in particular, the 
Hamiltonian operator H, the dipole moment 
operator wu, or the polarizability tensor P, and 
the w’s are suitable wave functions for the 
system. Such an integral can exist only if the 
reducible representation of its integrand contains 
the totally symmetric representation of the 
space group. From Eq. (2) we can then write 
for any of the specific operators mentioned 


(4) 
Consequently, we require that 
contains (5) 


where A,‘ is the totally symmetric representa- 
tion of the space group. This then demands that 


contains (6) 


since the product of two irreducible representa- 
tions cannot be totally symmetric unless they 
are identical (or complex conjugates). In par- 
ticular, such integrals can always exist in the 
special case where 


TY(M) contains (7) 


Similarly, in connection with the mixing of 
states by perturbations, we shall encounter 
integrals of the form 


f Mab dr, (8) 


and by the same reasoning we require that 
contains PS) PS(¥,). (9) 


TABLE II. 


Such integrals may always exist in the special 
case where 


contains 


Ill. THE HARMONIC APPROXIMATION 


If the potential function for the system in- 
cludes only quadratic terms, it is known from 
the theory of small vibrations that it is possible 
to find coordinates (which are linear combina- 
tions of the original coordinates) such that the 
kinetic and potential energy can be put in the 
form 


(11) 
2 V= (12) 


In this case the motion of the system can be 
resolved into the motions of a set of independent 
harmonic oscillators of angular frequency ; in 
the normal coordinates Q;. The wave function 
for the crystal is then 


= (13) 


when the functions, ~,:, are given explicitly by 


Wns(Qs) =expl — 
X Anil (14) 


where H,,; is the Hermite polynomial of degree n. 
The energy levels in the harmonic approximation 
are given by 


(15) 


where 
Vi=i. 


If the Hamiltonian, H, and hence JT and V, are 
to be invariant under all of the operations of the 
space group of the crystal, the coordinates Q; 
must transform as the irreducible representa- 
tions of the group. The functions H,,; are odd or 
even functions of Q; according to whether 1; is 
odd or even. If Q; belongs to a non-degenerate 
irreducible representation, the corresponding Wns 
transforms either as Q; or the totally symmetric 
representation according as m; is odd or even. 

The selection rules in either the infra-red or 
Raman spectra are given, respectively, by the 
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non-vanishing of the components of the integrals 


Mr'n= f iWWnid Qi, 
and (16) 


The situation is the same in this respect as in 
free molecules except that the distinction be- 
tween the components, e.g., uz, wy, and wz, has 
physical significance. Hence, the selection rules 
may be different for different orientations of the 
crystal and different polarization of the radiation. 
As is customary in vibrational theory, we 
shall consider only linear terms in the expansion 
of v and P. In that case such integrals can 
exist in the harmonic approximation only if 
An;=+1, and the observed frequency is the 
same for all values of . From Eq. (4) we know 
that v and P have components only in the normal 
coordinates belonging to those irreducible repre- 
sentations of S which are also irreducible repre- 
sentations of U. Since the integrals of Eq. (16) 
involve only one normal coordinate, they can 
exist only for normal coordinates whose sym- 
metry type is some I”. That is, only those 
normal coordinates of the crystal which are 
totally symmetric with respect to translation 
may be active, and then only if the transition is 
allowed under U. The selection rules for funda- 
mentals are thus given exactly by Eq. (7). 
Explicitly, we may write any totally sym- 
metric coordinate of the crystal in the form 


where the sum is over all the unit cells and the 
coordinate in brackets might be called a normal 


coordinate of the unit cell. Physically, this. 


motion is that in which equivalent atoms of all 
unit cells are in phase. Such modes are relatively 
easy to compute and ‘have been used as an 
approximation.? We see that this approach is 
rigorous for the calculation of observed funda- 
mental modes. 

An immediate consequence of this one-to-one 
correspondence between the irreducible repre- 


7S. Bhagavantam and T. Venkatarayadu, Proc. Ind. 
Acad. Sci. A9, 224 (1939). 


sentations of the unit cell group and those 
irreducible representations of the complete space 
group which can be active as fundamentals in 
either the infra-red or Raman spectrum is the 
rule: in the harmonic approximation the total 
number and symmetry type of the observed fre- 
quencies in either the Raman or infra-red spectrum 
must be equal to the number and symmetry type 
determined by the molecules in a single unit cell. 

Since the number of atoms in a unit cell is 
finite and, in most cases of interest, small, we 
conclude that if the harmonic approximation is 
adequate, the infra-red and Raman spectra 
should consist of a few very sharp lines of breadth 
equal to the natural line breadth. This follows 
since we have already included all of the degrees 
of freedom of the system in our discussion and 
there are no others (such as rotational degrees 
of freedom in the gas phase) to couple with a 
vibration to yield a band rather than a line. 

In the simple case of cubic crystals, such as 
NaCl and MgO, we should then expect a single 
sharp line in the infra-red spectrum. Absorption 
measurements on these materials*® at room 
temperatures, as well as on a large number of 
other cubic crystals by other investigators have 
shown bands more than 100 cm-! wide (1000 
cm-! in the case of MgO) with a great deal of 
fine structure. This has already been ascribed 
by the authors to anharmonic coupling between 
normal modes and discussed theoretically by 
Born and Blackman” and by Seitz.® 

However, we can go much further than this. 
A considerable line breadth is found in all infra- 
red spectra of crystals," and in the relatively 
few cases which have been examined under high 
resolution such as, for example, NH,Cl and 
NH,Br®” and the 7y band of calcite,!* a wealth 
of fine structure is evident. Similarly, diffuse 
fundamentals are frequently found in the Raman 
spectra although, in general, the Raman lines 
are much sharper than the corresponding lines 


a 931} B. Barnes and M. Czerny, Zeits. f. Physik 72, 447 

°R. B. Barnes, R. R. Brattain, and F. Seitz, Phys. Rev. 
48, 582 (1935). 

10M. Born and M. Blackman, Zeits. f. Physik 82, 551 
(1933); M. Blackman, Zeits. f. Physik 86, 421 (1933). 

11 See, for instance, C. Schaeffer and F. Matossi, Das 
Ultrarote Spektrum (Verlag-Julius Springer, Berlin, 1930). 

#2 C, Beck, J. Chem. Phys. 12, 71 (1944). 

18S, Silverman, Phys. Rev. 39, 72 (1932). 
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in the infra-red spectrum. We conclude that the 
available data points to the inadequacy of the 
harmonic approximation. 

Nevertheless, we should like to develop it a 
little further as a basis for higher approxima- 
tions. It has been pointed out by Halford* that 
the local symmetries may be used to determine 
selection rules if the molecules other than the 
one considered are kept in their equilibrium 
position. Let us investigate the relation between 
the rigorous (in harmonic approximation) selec- 
tion rules obtained under the unit cell group 
and those obtained using the local symmetries. 

Suppose that the reducible representation of 
a physical property C, such as any of the compo- 
nents of uw or P, has been put in the form 


TC) (18) 
k 


Then, from Eqs. (7) and (16), this property 
may be observed with those frequencies, vi, 
whose corresponding normal coordinate, Q;, be- 
longs to an irreducible representation, T;”, for 
which ¢,” is not zero. If we now represent the 
irreducible representations of U, considered as 
vectors whose components are the corresponding 
characters x,7(R), in terms of the irreducible 
representations, I';’, of the rth site group L’, 


we write 
(19) 


where 
be? = = (R)xi7(R)/h’. (20) 


This is just an analytical statement of the 
correlation chart. The 0,;7 are non-vanishing 
when a line connects [;* and I,2. Since this 
transformation will have applied to the normal 
coordinate and the physical property, they will 
both transform as any I,’ for which 3;;:740 and 
hence would be predicted to be active under L’. 
That is, every normal coordinate which may be 
active under the unit cell group must also be 
active under any of the local symmetry groups. 

Let us now examine the converse situation. 
Suppose that we consider the crystal from the 
point of view of Halford.* Then we divide the 
atoms of the crystal among molecules M’, with 
local symmetry L’. From this point of view we 
consider that all molecules, M*, of a symmetri- 
cally equivalent set have the same frequencies, 


so that in order to count frequencies we need 
consider only a single molecule of each set. If a 
reducible representation of a molecule, M’, is 
I'"(M"*), the number of distinct frequencies, 1,’, 
occurring in I,” is given by 


ni? = (21) 


where h’ is the order of L". These frequencies 
may be active in the infra-red or Raman spec- 
trum, respectively, if the reducible representation 
of v or P contains I',;’. According to Halford the 
selection rules for the crystal are obtained if 
Eq. (21) is written for each symmetrically 
non-equivalent kind of molecule. 

A difficulty arises in applying these selection 
rules to the separate components of v or P as 
we should like to do when considering polarized 
radiation. The reason is that if the x, y, and z 
coordinates are fixed in space (and hence in the 
crystal), the various symmetrically equivalent 
molecules, M’, of a given set-are not, in general, 
similarly oriented with respect to these axes. 
Consequently, although Eq. (21) will be shown 
to yield the proper number of active frequencies 
in either the infra-red or Raman spectrum, it 
must be modified if the selection rules with 
regard to the separate components of v and P 
are to be obtained when the components are 
designated by space-fixed axes. 

This is readily done if we remember that the 
reference axes of the groups U and S are fixed in 
space. Suppose that according to Eq. (18) a 
given component of v or P transforms as TZ. 
The number of times each I,’ occurs in I,” is 
given by Eq. (19). Using Eq. (21) we then find 
that the total number of distinct frequencies 
occurring in T',2 (and which are therefore active 
and have the polarization of the given component 
of v or P) resulting from molecules of type M’ is 


‘= 


The selection rules in this form are plainly 
equivalent to those expressed in Eq. (21). In 
the first case the reducible representation of v 
or P is reduced directly into the irreducible 
representations of L* while in the second the 
same reduction is effected in two steps. The 
total number of distinct frequencies occurring in 
T,2 due to all molecules is then, from the 
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Halford point of view, 
; nv= (22) 


We shall now show that, in the harmonic 
approximation, these selection rules are correct. 
We make no assumptions whatever about the 
degree to which the motions of the molecules are 
coupled. It has previously been shown that the 
selection rules for the crystal may be obtained 
from the unit cell group and the molecules within 
a single unit cell. The unit cell must contain 
h” /h" molecules of type since there are hY/h’ 
equivalent site groups, L’ (conjugate subgroups). 
We shall formally designate this set of symmetri- 
cally equivalent molecules by 


A reducible representation of the entire unit cell 
may be designated by 


The number of distinct frequencies, m2, occur- 
ring in some T’,” is 


Since none of the sets with different r are 
symmetrically equivalent, we can write 


Furthermore, the character of any transforma- 
tion of a set of symmetrically equivalent mole- 
cules, ((h"/h")M’), is zero if the transformation 
is not included in L’, since all such operations 
interchange members of the set (i.e., any opera- 
tion of U which leaves the center of gravity of 
a molecule, M’, invariant implies a corresponding 
symmetry element passing through the center 
of gravity. The operations of such elements are 
included in L’ by definition). Noting the defini- 
tion of a dot product given in Eq. (20), we have 


= (25) 
i 
Equation (23) then becomes, after substituting 
Eqs. (24) and (25) 
= (26) 


This is just the expression obtained in Eq. (22) 
by considering isolated molecules with the appro- 
priate local symmetry. We had previously noted 
that the total number of frequencies active in 
either the infra-red or Raman spectrum must be 
the same whether computed from Eq. (21) as in 
Halford’s paper or from Eq. (22). We can now 
state the following conclusion which is generally 
valid in the harmonic approximation ; the number 
of frequencies active in either the infra-red or 
Raman spectrum calculated from the local sym- 
metries of the molecules (including each non- 
equivalent type of molecule once) is exactly equal 
to the number calculated on the basis of the unit 
cell, which is equal to the allowed number for the 
entire crystal. It is not possible to calculate the 
polarization properties from local symmetry 
alone, but using the coefficients 5;;” they can be 
obtained from Eq. (26). 

Since the group L’ is a subgroup of U, an irre- 
ducible representation of U cannot contain an irre- 
ducible representation of L’ of higher degeneracy. 
Consequently, no degenerate vibration obtained 
in the local symmetry approach can ever be split 
when the entire crystal is considered. However, 
we have previously pointed out that several 
irreducible representations of U may contain the 
same irreducible representation of L’. They 
differ physically in the phase in which equivalent 
molecules are coupled. If the polarizability occurs 
in [,”, the dipole moment occurs in I¥, and 
both };;7 and b;;7 exist, then a mode which was 
active in both the infra-red and Raman spectra 
under I,’ might appear with different frequencies 
in the two spectra when the entire unit cell is 
considered. The difference in frequency is a 
measure of the effect of coupling between 
molecules. 

In order to see the problem more explicitly, 
let us designate the displacement of the center 
of gravity of each ion or molecule as well as the 
displacement of the angles which define the 
orientation of a set of axes fixed to each molecule 
(with respect to axes fixed in the entire crystal) 
by the coordinates s;;. Here the superscript 
refers to the kind of molecule and the subscripts 
to the specific molecule and coordinate (e.g., the 
ith coordinate of the jth molecule of the rth 
kind). These coordinates we shall call lattice 
coordinates, and they are analogous to the 
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coordinates which define the translation and 
rotation of a free molecule. Six of them are 
required to describe each non-linear ion or 
molecule (five of them are required in the linear 
case and three for a monatomic ion). 

Let us describe the distortion of the configura- 
tion, relative to the center of gravity and body 
axes, of the rth kind of molecule by the internal 
coordinates, m;;". These are internal coordinates 
in the same sense as in the usual molecular 
vibration problems. There are then 3n—6 in- 
ternal coordinates describing a non-linear mole- 
cule or ion containing atoms. We may use an 
equivalent set of coordinates, both lattice and 
internal, for every symmetrically equivalent 
molecule of a given kind. 

In these coordinates, the potential energy is 
given in the harmonic oscillator approximation 
by the expression 


nm + Vine (27) 


The first term is the potential function of the 
lattice if there is no interaction with the internal 
vibrations of the molecules (or, alternatively, if 
the molecules are considered rigid bodies). The 
second term is the sum of the potential functions, 
V,’, for the internal motion- of the individual 
molecules when there is no interaction between 
the internal motions of the various molecules 
nor between the internal motions and the lattice 
vibrations. The remaining terms are small; the 
third represents the coupling of the internal 
modes of different molecules, and the fourth 
represents interaction between lattice and in- 
ternal vibrations. 

In terms of these same coordinates, the 
kinetic energy is of the form 


= Tas’ + (28) 
ra 


where &;," and m,;" are the velocities correspond- 
ing to ss” and m,* and the coefficients are 
functions of the coordinates and masses. The 
first term is the translational and rotational 
energy of the molecules considered as rigid 


bodies, the second is the internal kinetic energy 
of the molecules and is analogous to that 
encountered in any molecular vibration problem, 
and the third represents the interaction between 
the internal angular momentum of each molecule 
and the overall angular momentum. Hence, only 
three of the six coordinates, s, of each molecule 
appear in this term. 

It is apparent empirically that the interaction 
terms, Vinm: and Vs, are small. If they were 
neglected entirely the problem would be sepa- 
rable, both classically and quantum-mechani- 
cally, into the problems of the lattice vibrations 
and the indtvidual molecules. The lattice vibra- 
tion problem may be handled by methods de- 
veloped by Born, Blackman, and others.'* We 
shall not pursue it in detail except to note some 
general behavior. It is convenient to replace the 
original translations and rotations, s;;", of a given 
molecule by linear combinations, g;"(I'."), which 
transform as irreducible representations of L’. 
In the case of these lattice coordinates the 
interaction terms with the lattice coordinates of 
other molecules are strong, so that these coordi- 
nates are not even approximately normal co- 
ordinates. Using the correlation chart between 
Lr and U, or Eq. (20), we can form h/h’ 
orthogonal linear combinations from the equiva- 
lent coordinates, of the symmetri- 
cally equivalent molecules having symmetry L’. 
These coordinates are symmetry coordinates of 
the unit cell group and transform as some Ij’. 
We shall denote them by g,¥(T,”). If the original 
set of coordinates was orthogonal, there is no 
ambiguity in this procedure, even in the case of 
degenerate representations. 

We have previously shown that the lattice 
modes which are active in the infra-red and 
Raman spectra are totally symmetric with re- 
spect to translation. The sum over all unit cells 
of gi¥(Tj”) is a totally symmetric (translation) 
symmetry coordinate of the crystal and starting 
from those IY which may be active under the 
unit cell group we can immediately calculate the 
active normal lattice vibrations. Actually, the 
same calculation is made more simply by noting 
that in these totally symmetric modes corre- 
sponding coordinates in all unit cells are equal 


See, for example, M. Born, Atomtheorie’ des festen 
Zustandes (B. G. Teubner, Berlin, 1923). 
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in any coordinate system. The original kinetic 
and potential energies may therefore be expressed 
in this case in terms of the coordinates of one 
unit cell. The secular determinant arising with 
these coordinates gq: (I'j”) is factored as far as it 
can be from symmetry considerations. Its solu- 
tion yields ‘“‘normal coordinates of the unit cell’’ 
which are linear combinations of g:7(T,%2) be- 
longing to the same irreducible representation, 
T,.”. The frequencies are correct for the entire 
crystal. 

If Vinm: and V,», are neglected, we must also 
solve the internal vibration problem once for 
each kind of molecule. The usual methods for a 
molecule whose Hamiltonian function is 


Ht=V'+T" (29) 


may be used.'® (We have omitted the subscript 
since the identity of the molecule is immaterial.) 
This problem is essentially the same as that of 
the free molecule except that H* must have the 
local symmetry, LZ’, of the molecule in the 
crystal. In general, this symmetry will be lower 
than that of the free molecule, the group L’ 
being a subgroup of the symmetry group of the 
free molecule. The solutions to this problem are 
a set of frequencies, v;,”, and normal coordinates, 
which must transform as We have 
previously shown that the selection rules in 
terms of these coordinates are rigorous. The 
frequencies are very nearly correct if the neg- 
lected terms are small. In the absence of specific 
intramolecular coupling such as hydrogen bond- 
ing this approximation is quite satisfactory. 

As a matter of fact, since T’ is identical with, 
and V’ is only slightly different from the case of 
a free molecule, even these frequencies will 
usually be very nearly the same as those in the 
free molecule, although degeneracies may be 
split. Similarly, although selection rules are 
relaxed by the lowered symmetry, the new 
frequencies need not appear with appreciable 
intensity. If the normal coordinates and the 
corresponding frequencies of the free molecule 
are known,.the change in frequencies may be 
calculated by the usual perturbation methods. 

From the sets of equivalent Q;"(I'.") of mole- 
cules in a single unit cell we can, following Eq. 


15 E. B. Wilson, ds J: Chem. Phys. 7, 1047 (1939); 9, 
76 (1941); 15, 736 (1947). 


(20), construct 4”/h’ orthogonal symmetry co- 
ordinates, S,;“(T;”), of the unit cell. In the limit 
of vanishing interaction terms these are true 
“normal coordinates’ of the unit cell. If the 
cross terms are included, rewritten in terms of 
the coordinates and they will 
consist only of products of coordinates of the 
same species, ',.”. The exact normal coordinates 
will then contain small admixtures of other 
symmetry coordinates of the same _ species. 
Nevertheless, because the coupling is weak the 
normal modes of the unit cell will still be divided 
into a set whose frequency is very nearly the 
same as the Q;’(T.") which we shall continue to 
call internal modes, and a low frequency set 
whose frequencies are nearly those of the g;"(T;”) 
which we shall still term lattice modes. 

The normal coordinates of the unit cell have 
frequencies identical with those modes of the 
entire crystal which are totally symmetrical 
with respect to translation (and which may be 
observed in the infra-red or Raman spectra). 
The calculation of the other crystal modes 
(which are inactive, however, as fundamentals 
in either the infra-red or Raman spectrum) is 
enormously difficult. Because of the weakness of 
the interactions we may anticipate that all of 
the high frequency crystal modes arising from a 
single internal mode will lie within a fairly 
narrow frequency range. 

A comparison of the observed Raman spec- 
trum with that expected in the three cases: the 
hypothetical free molecule, the molecule using 
local symmetry, and the molecule using unit 
cell group symmetry, has been made by Couture” 
for single crystals containing CO3;-, NO;-, and 

To summarize: in the harmonic approximation 
the spectrum of a crystal containing molecules 
or complex ions should consist of a small number 
of very sharp lines. The selection rules for those 
frequencies which correspond to internal modes 
of the molecules are obtained in this approxima- 
tion by considering the individual molecules, 
using their actual local symmetry in the crystal. 
The observed frequencies may be calculated 
exactly by considering only one unit cell, with 
the aid of the unit cell group which is the factor 


16L. Couture, Ann. de physique 2, 5 (1947), ser. 12. 
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group of the space group obtained when all 
translations are considered an identity operation. 

When the harmonic approximation is valid 
there is some advantage in the study of molecules 
in the crystalline state. The extreme sharpness 
of the lines in the spectrum would reduce the 
overlapping of bands as compared to the gas 
spectrum. The additional information which may 
be gained from polarization measurements, both 
in the Raman and infra-red spectrum, because 
of the fixed orientation of the molecules in space 
simplifies the analysis of the observed spectrum. 
This approximation may be expected to become 
increasingly valid as the temperature of the 
specimen is lowered, and the sharpening of the 
lines has indeed been observed.!” 

These advantages must, of course, be balanced 
against the disadvantage of reduced symmetry 
and the complications which are introduced by 
anharmonic interaction with the lattice vibra- 
tions of the crystal. 


IV. THE EFFECT OF ANHARMONICITY 


Just as in the case of free molecules, the effect 
of anharmonic terms in the potential function is 
to introduce overtone and combination fre- 
quencies into the spectrum. As far as the internal 
vibrations of the molecules in the crystal are 
concerned, we do not expect the situation to be 
particularly different from that of free molecules. 
An important qualitative difference may arise 
in the spectrum of the crystal, however, in that 
we may also expect combination levels to occur 
between internal and lattice vibrations. Since 
the lattice modes of most crystals are of fre- 
quency lower than 200 cm~, such combinations 
do not appear as separate bands, but convert the 
sharp line predicted in the harmonic approxima- 
tion to a band which, depending on the temper- 
ature, may or may not be centered about the 
fundamental. Secondly, because of the enormous 
number of lattice modes, all of which yield 
combination lines in the same narrow region of 
the spectrum, the total intensity due to combi- 
nation frequencies may be very considerable. 
In the case of a fundamental which is forbidden 
but whose combination bands with the lattice 


17 See, for example, W. H. Avery and J. R. Morrison, 
i: App. Phys. 18, 960 (1947); and R. S. Krishnan, Proc. 
nd. d. Sci., A26, 932 (1947); A27, 321 (1948). 


are allowed, the result may be to introduce a 
spectrum in the immediate vicinity of the for- 
bidden fundamental. Such a case may occur in 
the doubly degenerate, 5.7-4 band of the NH,*+ 
ion in NH,Cl and NH,Br which, although for- 
bidden, appears with considerable intensity and 
a great deal of structure. 

The theory of the effect of anharmonicity in 
crystals is identical in principle with that for 
free molecules.'* If V’* represents the anharmonic 
part of the potential function, we may get inter- 
action of the two states ¢;* and ¢,° of the crystal 
if the integral 


Hj = f (30) 


exists. Since the potential, V’*, must be totally 
symmetric, the product ¢;°¢;* must be so too. 
Consequently, both states must be of the same 
symmetry species with respect to the group, 
the unit cell group and the site groups. The 
symmetry discussion can always be reduced, 
then, to the site groups where all lattice coordi- 
nates transform as simple translations or rota- 
tions. The situation is not as simple for quanti- 
tative calculations, however, because the lattice 
coordinates of a given molecule are closely 
coupled to those of other molecules. We cannot, 
on that account, write any meaningful lattice 
wave functions for a single molecule. In Section 
VII we shall discuss an approach which simplifies 
this discussion, both mathematically and in that 
it provides a simple physical picture, so we shall 
not develop this approach further. 

It is worth pointing out, however, that the 
symmetry properties of the states, wni(T;°), 
which enter the ¢’s, according to Eq. (13), are 
not complicated in the case of internal modes, 
since at room temperature we are interested 
only in the ground state and the first excited 
level. However, the low frequency lattice modes 
may have quite complex symmetries if they are 
degenerate, since even at room temperature 
quite high values of the quantum numbers may 
be important. The symmetry of excited degen- 
erate states has been given by Tisza.!® 


18G. Herzberg, Infrared and Raman Spectra (D. Van 
Nostrand Company, Inc., New York, 1945). 
1” L. Tisza, Zeits. f. Physik 82, 48 (1933). 
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V. SELECTION RULES FOR THE 
ANHARMONIC CRYSTAL 


Let us consider first the selection rules in- 
volving only one normal coordinate. In the 
general case we can use Eq. (16), but without 
the restriction that An= +1. Then from Eq. (7) 
we conclude that all such transitions are allowed 
which are allowed under (U) (or, from our 
previous discussion, under L’. However, a further 
set may be allowed by the more general Eq. (6), 
but only if the transition is also allowed under U, 
and L’, since the normal coordinates under any 
I'S may be written as a linear combination of Q.”. 
Consequently, the selection rules may be deter- 
mined from U alone, but instead of a sharp line 
we may expect a band. In particular, since all 
first overtones are allowed in the Raman spec- 
trum, we expect that they may be diffuse. 

In the case of combination bands in which 
changes in two quantum numbers may occur, 
the selection rules depend on the non-vanishing 
of the integrals 


m= f (31) 


or 


p= f (32) 


The most important cases are certainly those 
for which 7’=i+1 and j’=j+1. Just as before, 
we can say that if these transitions are to occur, 
they must be allowed under U, and L’, according 
to Eq. (10). If they are allowed under U, they 
may also be allowed under Eq. (9), in which 
case the line will be broadened into a band. In 
the case in which y; and y; both refer to internal 
vibrations, these transitions give rise to the usual 
combination bands. If y; refers to an internal 
vibration and y; refers to a lattice vibration, 
these transitions have the effect of adding struc- 
ture to the simple transition in y;, if it is allowed, 
or of producing an apparent violation of selection 
rules if the simple transition is forbidden. The 
situation is complicated, of course, because at 
room temperature the crystal is not normally in 
a ground state as regards lattice vibrations but 
may be in a rather high level. Consequently, a 
multiplicity of lines may be expected due to 
anharmonicity. 


We can proceed similarly for higher combina- 
tions, and the number of allowed frequencies 
increases rapidly, but in the absence of more 
empirical information, it would be of little value 
to go further. 


VI. THE EFFECT OF TEMPERATURE 


If a crystal is studied at sufficiently low 
temperatures, both the lattice and molecules 
may be regarded as in their ground state. The 
infra-red or Raman spectrum obtained under 
these circumstances is as simple as it can be for 
the crystal since all anharmonic effects are 
minimized. If they make no contribution, each 
line should be exceedingly sharp and the selection 
rules derived from the harmonic approximation 
should be valid. Even if anharmonicity is present, 
difference bands are not allowed, so that all 
combination bands must appear on the high 
frequency shoulder of the fundamental. 

As the temperature is raised, we’ may expect 
the molecules to remain largely in the ground 
state. The population of upper lattice states will 
increase. Consequently, difference bands will 
appear on the low frequency shoulder and at 
sufficiently high temperatures the band should 
become symmetrical. Because of the increasing 
effect of anharmonicity as high lattice modes are 
excited, we may expect the line width to increase 
further and the intensity of the combination lines 
to rise. 

Secondly, forbidden lines may appear. The 
intensity of such lines may be expected to 
increase as the temperature is raised. The 
existence of a line whose total intensity depends 
strongly on temperature would argue for large 
anharmonic terms, since in the harmonic ap- 
proximation the intensity should be independent 
of temperature. This effect may make it possible 
to identify ‘‘forbidden”’ lines. 

It seems probable that the study of the 
temperature dependence of spectra will be of 
critical importance in the'elucidation of crystal 
dynamics. 

VII. ANOTHER POINT OF VIEW 


We should like to reduce the discussion of the 
effect of lattice vibrations on the internal modes 
of a molecular crystal to a discussion involving 
only individual molecules, in the same fashion as 
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we found possible in the harmonic approxima- 
tion. This was not possible from the point of 
view of ordinary anharmonic theory. The reason 
is that although the problem could easily be set 
up in terms of the vibrational wave functions of 
individual molecules, Wn;:” (since internal modes 
are very weakly coupled), the lattice wave 
' functions involve coordinates throughout the 
crystal and are not conveniently discussed from 
a local point of view except as regards their 
symmetry properties. 

That some of the interaction, at least, can be 
treated from a local point of view seems apparent 
physically in the following way. If the frequency 
of all lattice vibrations is sufficiently low com- 
pared to the frequency of the internal vibrations 
of the molecules, the lattice coordinates may be 


considered to be essentially constant during 


several cycles of the internal motion.. However, 
the lattice coordinates would affect this motion 
in that they would determine the symmetry 
and quantitative magnitude of the potential field 
in which the molecular vibration takes place. 
In particular, when all lattice coordinates have 
the value zero, the crystal is in its equilibrium 
configuration and the symmetry of the Hamil- 
tonian of each molecule (H’) is L’. For any other 
value of the lattice coordinates, this symmetry 
will, in general, be lowered. In many cases, this 
effect may be visualized qualitatively. If we 
imagine, for example, a tetrahedral molecule 
such as NH, at the center of a cubic unit cell 
with the N—H bonds directed along threefold 
axes of the cube toward four cube corners, the 
local symmetry of the NH,* ion (if similarly 
oriented in all unit cells) is Tg. Such a crystal 
possesses only two types of translational lattice 
modes, those in which the motion of the ions at 
the corners of the cube is parallel to that of the 
NH," ion and those in which it is antiparallel. 
The first mode is triply degenerate and is 
equivalent to a translation of the entire unit cell. 
Depending on the relative phase of the motion 
in different unit cells, it gives rise to the various 
acoustical modes of vibration of the crystal. 
Particularly in the case of long wave-length 
modes of this type, we expect no perturbation of 
the local potential field. The other gives rise to 
the optical modes of the crystal and is observed 
as the Reststrahlen frequency (in the case of 
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NH.,Cl it has a frequency of 183 cm~). This 
mode is also triply degenerate and it may be 
considered to be composed of the three vibrations 
perpendicular to the cube faces. If, in first 
approximation, these are considered separately, 
it is apparent that such a displacement of the 
center of gravity of the NH,* ion reduces the 
symmetry from Tz to C2. In this case, both of 
the triply degenerate modes of the NH,* ion 
might be expected to be split and the inactive 
double degenerate mode might become active. 
The third lattice mode is the triply degenerate 
torsional oscillation of NH,*+ about the same 
three axes. In the same approximation, it might 
reduce the symmetry to C». 

Conversely, if the frequency ratio is great, 
the internal vibrations would effectively be 
averaged so that in calculating lattice frequencies 
we need include only the mean molecular con- 
figuration. Since this mean configuration would, 
in general, vary with the lattice coordinates, the 
internal molecular energy (and frequencies) 
would be a function of these parameters. This 
function then constitutes a further potential 
energy term in the computation of lattice modes. 

This type of problem was solved long ago by 
Born and Oppenheimer” for the case of electronic 
and nuclear motion. Their approach has been 
widely applied to the coupling of electronic and 
vibrational states in molecules." A great deal of 
this work can be applied without change to our 
problem if the wave functions are appropriately 
renamed. 

Let us write the complete potential function 
in the form 


Vins") + (33) 


where V;* is that part including only lattice 
coordinates, V;’ is that part including only 
internal coordinates of the jth molecule of local 
symmetry L’, Vjms;" includes those terms con- 
taining both internal and lattice coordinates of 
the jth molecule, Vinm* represents (as in Eq. 
(27)) cross terms between internal coordinates 


20M. Born and R. Oppenheimer, Ann. d. Physik 84, 457 
(1927). Another mathematical technique which should be 
applicable to this problem has recently been’ given by R. 
(i9 a and L. H. Thomas, J. Chem. Phys. 16, 310 

21H. Sponer and E. Teller, Rev. Mod. Phys. 13, 85 
(1941). 
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of different molecules, and V,,,’* represents any 
terms not previously included. Following our 
discussion of the harmonic case, we begin by 
ignoring the last two terms which are certainly 
small. The kinetic energy is given classically by 
Eq. (28). Since we are looking for relatively large 
effects, we shall also ignore the last term of Eq. 
(28) which is known to be small. The Schrédinger 
equation is then of the form 


[Hy Vet) (34) 


where the H’s are the corresponding Hamiltonian 
operators and W7*(m,s) is the total wave function 
of the crystal. If all of the frequencies associated 
with lattice coordinates are small compared to 
any of the frequencies associated with internal 
coordinates, we may, following Born and Oppen- 
heimer, set 


Vr*(m,s) ‘Wm*(m,S), (35) 


where ¥»°(m,s) is a solution of the equation 


in which the lattice coordinates, s, are regarded 
as parameters, and y;‘(s) is a solution of the 
equation 


(37) 


Since Eq. (36) is separable, we have also 


En(s) = in’(s) (38) 
Wm*(m,S) T]¥m"(m,s) (39) 


Vims’ —Ejm(S) Wim"(m,s)=0. (40) 


If the functions y;*(s) and Wm*(m,s) are de- 
generate, Eq. (35) must be replaced by a linear 
combination of similar terms. 

This is discussed by Born and Oppenheimer 
and by Sponer and Teller. The quantitative 
validity of this approximation will be discussed 
at a later date. If the frequency ratio is very 
great, it is certainly quite valid, but in cases in 
which some of the internal molecular modes are 
of very low frequency, it is certainly not appli- 
cable. Even in these cases, a similar approach 
can still be used in discussing the perturbation 
of the high frequency molecular modes. 

We are particularly interested in the solutions 
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of Eq. (40) which yield the molecular energy 
levels (and frequencies) as a function of the 
lattice coordinates. We need solve this equation 
only once for each type of molecule present so 
that we will omit the subscript 7. The quantity 
H* may be taken to be identical with that in 
Eq. (29) since, for our purposes, the harmonic 
approximation is adequate for the internal vibra- 
tions. The quantity V,,” is in the nature of a 
perturbation, and we shall take it to include 
harmonic coupling terms as well as anharmonic 
terms of any functional form whatever, subject 
only to the general symmetry restriction that it 
must have the symmetry properties of the 
group L’. 

We shall use as the basis functions for a per- 
turbation treatment the functions y,:{Q;(T,") ] 
which we shall designate hereafter simply as 
Wniz”(m). These functions are solutions of the 
equation 

(Ht — Enix” =0, (41) 


and have the form given by Eq. (14) if the 
coordinates Q;([,") are the normal coordinates 
of the unperturbed molecule. The corresponding 
frequencies are just those of the local symmetry 
approximation which we have previously dis- 
cussed. 

According to the usual methods of perturba- 
tion theory a state Wnix”(m) may be perturbed 
by those states Wn’i7.-"(m) for which the integrals 


exist. In the case of small amplitude lattice 


TABLE III. The splitting of degenerate molecular vibrations 
by crystal lattice vibrations. 


' Irreducible Lattice mode responsible 
Site group representation for splitting 
C3 E T(E); R(E) 
Cov E T(E); R(E) 
Car, Ey’, Ey” T(Ey’) 

3 E RED R(E) 
Dsa E,, Ex R E,) 

E T(B) 

Ssv(Dea) E T(B)) 
E Not split 
F T(F), R(F) 
Ta E Not split 
Fi, Fe T (Fo) 
Th Not split 
o Fu R(F,) 
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vibrations we may expand V,,,” in terms of the 
symmetry coordinates q,’(I',"), which describe 
the rotations and translations of the molecule in 
question. We thus obtain 


Vins’ = LL (43) 


If the amplitudes are small the linear terms are 
most important. Equation (42) then becomes 


H' 
x f (44) 


The states may perturb each other if any of the 
integrals exist. Since 
must be totally symmetric, each a,j; which is a 
function of the internal coordinates, must have 
the same symmetry, I,’, as qg;"(I',"). The integrals 
can thus exist only if contains 

The perturbation will only be important if 
Enix” — Eni? is small and will be most important 
in the case of degenerate states. The perturbation 
will split the degenerate level in such a case. 
Since the symmetry of the product of the wave 
functions is the (I;,”)?, the degeneracy may be 
split if (T,")? contains T',’. We conclude that in 
this approximation a degenerate mode of vibra- 
tion of a molecule in a crystal may be split if the 
square of the symmetry species of the degenerate 
level contains the symmetry species of transla- 
tion or rotation. We assume, however, that the 
molecule is stable with respect to totally sym- 
metric lattice vibrations. 

A very similar discussion has been given by 
Jahn and Teller” for the perturbation of degen- 
erate electronic states by molecular vibrations. 
From their tables, it is immediately possible to 
predict which degenerate irreducible representa- 
tions may be split. In particular, it is found that 
in this approximation degenerate modes of the 
following local symmetry groups are not split by 
lattice vibrations: C4, Cav, Car, Da, Dan, Co, Cov, 
Cos, De, Den, O, and O,. The symmetry types of 


2H. A. Jahn and E. Teller, Proc. . Soc. A161, 220 
(1937). 


vibrational modes which may be split are tabu- 
lated in Table III. The modes responsible for 
the splitting are designated T or R depending-on 
whether they are translational or rotational in 
origin, and their symmetry species is indicated 
in parenthesis. On this basis we would predict 
that the degenerate vibrations of CO;= in calcite, 
dolomite, etc., may be split. In fact, both de- 
generate modes are split by approximately 50 
cm (see reference 11, pp. 338 ff.). Similarly, it 
is predicted that the triply degenerate modes of 
NH,* may be split in NH,Cl but that the 
doubly degenerate mode may not. This has 
been found to be the case, the two triply degen- 
erate modes being split 35 cm— and 50 cm™, 
respectively.” 

The only other case in which effects of con- 
siderable magnitude might occur is that in which 
a degenerate mode of the free molecule is slightly 
split by the equilibrium symmetry of the crystal. 
In this case the two states interact if the product 
of their symmetry species contains the symmetry 
of either rotation or translation. The conse- 
quence, in all such cases, will be to increase the 
separation between the two levels. 

It seems fruitless to carry the discussion 
further at this time. The general nature of the 
transitions is the same as those discussed by 
Sponer and Teller. These same authors discuss 
selection rules in this approximation. It is only 
necessary to note that in their discussion we 
must substitute molecular vibrations for elec- 
tronic states and lattice vibrations for molecular 
vibrations to make it applicable to our case. 

It is even reasonable to enunciate a principle 
similar to the Franck-Condon principle, namely, 
that the position and velocity of the center of 
gravity of the molecules will not change during 
a transition between internal vibrational states. 
Finally, we should emphasize again the impor- 
tance of temperature. The effects we have just 
discussed yield energy terms which are linear in 
the lattice vibration amplitudes. Hence, they 
too are minimized at low temperatures. 


2 R. Pohlman, Zeits. f. Physik 79, 394 (1932). 
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The heat capacity, entropy, internal energy, and free energy are calculated to 2500°C for 
tritium and tritium hydride. The distribution of tritium between the ortho and para states 
is considered. The equilibrium constant for the formation of tritium hydride from hydrogen 
and tritium is calculated from 50°K to 2500°K. The dissociation of tritium and tritium hydride 
into the atoms is discussed. The bearing of the radioactivity of tritium and tritium hydride 
on the attainment of thermodynamic equilibrium is briefly considered. 
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CCURATE calculations of the thermo- 
dynamic properties of hydrogen have been 
carried out by Giauque,! and similar calculations 
have been made for deuterium and deuterium 
hydride by Johnston and Long.? The present 
paper presents corresponding work for tritium 
and tritium hydride. The values refer to a mole 
of the ideal gas. 

The general methods of such calculations are 
well known. In the present calculations partition 
functions and their derivatives have been calcu- 
lated by direct summation over rotational and 
vibrational states. The energies of T; and HT 
are representable in the form, 


Ey, =w,(v+3) —X w.(v+})? 
+6-(v+3)* JJ(J+1) 
JJ?(J+1)? 
(1) 


The constants are listed in Table I. They were 


TABLE I. Molecular constants of Tz and HT. 


We 3608.3 2553.8 

X eve 87.585 43.872 

40.5747 20.3243 

Qe 1.67048 0.59222 
0.007636 


0.0655 X 10-8 


F, 1.5378 X 0.19328 X 


1 Giauque, J. Am. Chem. Soc. 52, 4826 (1930); see also 
C. O. Davis and H. L. Johnston, J. Am. Chem. Soc. 56, 
1045 (1934). 

2H. L. Johnston and E. A. Long, J. Chem. Phys. 2, 
389, 710 (1934). 


ish and H. A. Bethe, Rev. Mod. Phys. 9, 373 
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computed from the data for ordinary hydrogen 
as given by Jeppeson,* who expresses the energy 
in a form slightly different from Eq. (1), through 
the theory of isotopic diatomic molecules.t The 
excellent agreement of the constants so calcu- 
lated with the observations of Jeppeson® in the 
case of HD lends confidence to the present 
procedure. The agreement in the case of D,*® 
is not quite so good. It is felt that, in general, 
the calculations are correct to the number of 
significant figures given except at the highest 
temperatures where the large values of v and, 
particularly, of J needed for the convergence of 
the summations lie outside the range employed 
by Jeppeson in fitting the constants of Eq. (1). 
The constants in Table I were calculated 
taking the atomic weights (physical scale) of 
hydrogen and tritium to be 1.008130 and 
3.017050.7 The physical constants are otherwise 


TABLE II. Distribution of molecular tritium among ortho 
and para states. 


% ortho 
0.000 


% ortho 
66.803 


10 2.823 50 71.272 
15 16.429 75 74.475 
20 33.814 100 74.928 
25 47.446 125 74.990 
30 56.694 150 74.999 

62.783 75.000 


3C. R. Jeppeson, Phys. Rev. 44, 165 (1933). 

4 Jevons, Report on Band Spectra of Diatomic Molecules 
(Cambridge University Press, Teddington, 1932). 

5C. R. Jeppeson, Phys. Rev. 45, 480 (1934). 

6 C. R. Jeppeson, Phys. Rev. 49, 797 (1937). 

7J. Mattauch, Phys. Rev. 57, 1155 (1940); M. S. 
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TABLE III. Heat capacities and entropies (including nuclear spin) of T: and HT in calories per degree per mole. 


So Cp® So 
p-Te2 e-Te n-Te2 HT HT 


0.000 
10 4.967 4.967 6.770 4,967 18.860 
15 4.967 4.981 8.970 4.971 20.874 
20 4.967 5.103 8.688 5.001 22.304 
25 4.967 5.447 7.743 5.087 23.414 
30 4.998 5.994 7.040 5.247 24.321 
35 5.055 6.609 6.630 5.444 25.095 
40 5.152 7.156 6.420 5.653 25.776 
50 5.458 7.989 6.390 6.091 26.956 
75 6.332 7.661 6.656 6.664 29.342 
100 6.784 7.145 6.874 6.874 31.237 
125 6.925 7.008 6.945 6.945 32.769 
150 6.958 6.975 6.963 6.963 34.035 


6.966 


0.000 0.000 0.000 
14.495 14.873 18.886 4.977 16.043 
16.510 18.115 20.901 5.142 18.083 
17.956 20.697 22.334 5.607 19.619 
19.126 22.535 23.460 6.195 20.934 
20.166 23.881 24.400 6.676 22.108 
21.136 24.933 25.222 6.972 23.162 
22.056 25.804 25.963 7.110 24.104 
23.732 27.220 27.268 7.141 25.698 
26.827 29.845 29.830 7.023 28.568 
29.001 31.795 31.795 6.987 30.582 
30.577 33.338 33.338 6.977 32.140 
31.851 34.607 34.607 6.974 33.412 


34.487 


C. 0 
°K oT o-Ts HT HT 


So 
HT 


6.969 
225 6.971 37.432 6.974 36.239 
250 6.973 38.166 6.975 36.974 
275 6.975 38.831 6.976 37.639 
298.1 6.978 39.394 6.978 38.202 
300 6.978 39.438 6.978 38.246 
400 7.007 41.449 6.986 40.254 
500 7.077 43.019 7.002 41.815 


7.196 


46.430 7.181 45.136 
47.325 7.285 45.987 
7.401 46.760 
7.708 48.445 
7.993 49.877 
2000 8.779 53.916 8.442 52.242 
2500 9.011 55.901 8.758 54.159 


those of the International Critical, Tables® in 
order that these calculations may best be utilized 
in connection with the large body of similar 
calculations now existing. The experiments of 
Bloch, et al.° indicate the expected one-half unit 
of spin for the tritium nucleus. 


——o 


Fic. 1. Fraction of tritium molecules in odd rotational 
states as a function of temperature. 


8 International Critical Tables, Volume’I, p. 16. 
® Bloch, Packard, Graves, and Spence, Phys. Rev. 71, 
373(L), 551(L) (1947). 


THE EQUILIBRIUM DISTRIBUTION OF MOLECU- 
LAR TRITIUM AMONG ORTHO AND 
PARA STATES 


Like ordinary hydrogen, molecular tritium 
exists in ortho and para forms which occupy the 
odd and even rotational states, respectively, and 
which have the spin weights 3 and 1, respectively. 
The fraction of the molecules which are ortho 
attains the limiting value of 2 at moderate 
temperatures as may be seen from Table II and 
from Fig. 1. The limiting value is attained at 
lower temperatures than is the case for hydrogen 
and deuterium. 

It is well known that both hydrogen and 
deuterium gases show a very slow attainment of 
ortho-para equilibrium in the absence of cata- 
lysts at room temperature. It seems very likely 
that the radioactivity of tritium, giving rise to 
ionization, dissociation, and subsequent recombi- 
nation, would provide an internal mechanism 
whereby pure tritium could attain ortho-para 
equilibrium at a much more rapid rate than pure 
hydrogen or deuterium. Rough calculations indi- 
cate that the half-time for conversion by this 
mechanism should be of the order of a few 


days. 


Cp® 
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TRITIUM AND TRITIUM HYDRIDE 
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Fic. 2. Internal portion of the low temperature heat 
capacities at constant volume of ortho, para, and equi- 
librium tritium. 


THE HEAT CAPACITY AND ENTROPY 


The heat capacities at constant pressure of 
ortho, para, equilibrium, and normal tritium 
(room temperature ortho-para distribution) are 
given in Table III along with the corresponding 
total entropies which include the nuclear spin 
contribution of R1In4. Values are also given for 
HT. By 175°K the four types of T2 have attained 
identical heat capacities, and the entropies of 
normal and equilibrium tritium have become 


identical. Moreover, above this temperature, 


TABLE IV. E°—£,° in calories per degree per mole for 
T2 and HT (excluding trans ca. 


e-T2 n-Te HT 


85.47 0.00 
85.47 0.01 
85.48 0.35 
85.55 2.27 
6.94 

14.36 

23.72 

34.12 


Fic. 3. Internal portion of the low temperature heat 
Capacities at constant volume of equilibrium tritium, 
normal tritium, and tritium hydride. 


the entropies of para and of ortho tritium are 
less than the entropy of equilibrium or of normal 
tritium by very nearly R1ln4=2.754 e.u. and 
R 1n4/3 =0.572 e.u., respectively. This simplifies 
Table III. Figures 2 and 3 show the internal 
portion of Cy/R at low temperatures. The 
characteristic features are owing to rotation. 

It is well known that calorimetric measure- 


TABLE V. —[(F°—E °)/T] for e-T2, HT, and T in calories 
per degree per mole (spin effects included). 


e-T2 HT 


11.075 31.280 30.164 
13.092 31.937 30.814 
14.539 32.493 31.365 
15.689 3 32.537 31.408 
16.662 34.525 33.381 
17.517 36.073 34.918 
18.284 37.341 36.176 
19.612 38.420 37.245 
22.153 39.361 38.172 
24.021 40.197 38.994 
25.494 40.952 39.733 
26.711 42.572 41.312 
27.747 43.928 42.623 
28.649 46.129 44.744 
29.448 47.891 46.439 


1206.1 11670.0 
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‘ 
d T.°K eT: HT T 
26.208 
Te) 15 11.898 26.679 
20 13.792 27.080 
o-T2 p-T2 25 15.368 27.111 
d 0 113.96 0.00 0.00 
10 113.96 0.00 3.22 
at 40 18.736 30.554 
15 113.96 0.02 18.74 50 20.296 31.320 
25 114.01 1.75 55.00 100 25.016 32.563 
30 114.06 5.46 67.03 125 26.532 33.092 
0 50 118.28 47.35 97.91 100.56 55.89 200 29.745 36.535 
100 182.68 176.60 181.16 181.17 159.45 ; ° 
ly 125 230.25 228.94 229.92 229.92 209.81 oo 
to 150 279.70 279.45 279.64 279.64 260.01 
5 175 329.60 329.55 329.59 329.59 310.16 TABLE VI. Equilibrium constant for H2+T2=2HT. 
200 379.62 379.61 379.62 379.62 360.30 
m Te HT Te HT K 
ra 225 429.69 410.46 700 1436.2 1379.7 50 0.242 700 
re 250 479.81 460.64 800 1682.3 1596.6 100 0.947 800 
, 275 529.98 - 510.85 900 1944.7 1823.1 200 1.948 900 
li- 298.1 576.40 557.27 1000 2223.3 2060.5 298.1 2.57 1000 
A 300 580.21 561.09 1250 2983 2707 300 2.58 1250 
11S 400 782.44 762.49 1500 3816 3429 400 2.99 1500 
Ww 500 989.5 95.0 2000 5641 5061 500 3.26 2000 
600 2500 7608 6881 600 3.45 2500 
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TABLE VII. Dissociation constants of isotopic hydrogen 
molecules in atmospheres. 


TABLE VIII. Percentage dissociation of isotopic hydrogen 
molecules at a total pressure of one atmosphere. 


1000 1500 2000 2500 


T, °K 1000 1500 2000 2500 
He 3.60X10- 2.96K10-* 692xX10-* He 12.7X 10-8 9.5 10-4 1.31 
HD 2.48 1.57 1.36 3.28 HD 7.9 6.3 5.8 91 
D, 3.71 2.70 2.49 6.17 D2 9.6 8.2 7.8 1.24 
HT 2.26 1.49 1.32 3.21 HT 7.5 6.1 5.7 90 
Tz 3.01 2.44 2.35 5.94 Te 8.7 7.8 7.7 1.22 


ments performed on Hg: in the absence of a 
catalyst give rise to an entropy for the gas which 
is lower than the computed value by R1n4 
+3/4R1n3. These two terms correspond to 
nuclear spin degeneracy and to failure of ortho- 
para equilibration to occur in the time of experi- 
ment, respectively. This ‘‘discrepancy”’ is well 
understood and has been considered in detail by 
Giauque.! A similar situation exists for D,” 

Experimental determination of the entropy of 
tritium gas by the customary calorimetric tech- 
niques would be a difficult matter because of the 
fairly large amount of energy steadily being 
released by nuclear disintegration. Granting, 
however, that the calorimetric measurements 
could be made, we may speculate on the ‘“‘dis- 
crepancy” which would exist between the calori- 
metric and computed entropies for T: gas. The 
term R1In4 would remain even if calorimetric 
measurements on solid tritium could be extended 
to a small fraction a degree absolute. In spite of 
the fact that tritium may attain ortho-para 
equilibrium more rapidly than hydrogen, the 
additional discrepancy might be a sizeable frac- 
tion of the value 3/4R1n3 fully realized by 
hydrogen if the tritium were rapidly cooled to a 
low temperature and calorimetric measurements 
commenced. 


THE INTERNAL ENERGY AND THE 
FREE ENERGY 


The internal energy and the free energy, both 
given with reference to the state in which 
y= J=0 are given in Tables IV and V. Table IV 
does not include the contribution of 3/2RT 
for (E°—Ep°)transiation. The entries of Table V 
include the nuclear spin contribution of R 1In4 
for T, and HT and of R In4 for atomic T, arising 
from the half-units of nuclear and electron spin. 


10 Clusius and Bartholomé, Zeits. physik. Chemie B30, 
237 (1935). 


THE EQUILIBRIUM OF HYDROGEN AND TRITIUM 
WITH TRITIUM HYDRIDE 


Values of [(F°—E °)/T] for He, Tz and HT 
may be employed to calculate the equilibrium 
constant, K, for the reaction Hz+T.2=2HT from 
the relationship: 


—RinK= s(——) (2) 
T 


The values of —[(F°—£,°)/T] for Hz were 
taken from the tables of Giauque! at 298.1°K 
and at higher temperatures, while the values of 
Urey and Rittenberg," corrected to R=1.9869 
cal. deg.—! mole, were used below this temper- 
ature. The quantity AE»° is obtained from the 
zero point energies 6193.8, 3603.4, and 5072.7 
cal. mole! for He, T. and HT, and it has the 
value 348.2 cal. mole. The possibility of isotopic 
electronic shifts? ®® has not been considered in 
this and the following section. Table VI gives 
values of K at several temperatures. 


THE DISSOCIATION OF ISOTOPIC 
HYDROGEN MOLECULES 


The dissociation. constants of T, and HT have 
been calculated from a relationship of the form 
of Eq. (2). The value of AE»°® is found from 
D,(Hz) and the difference in the zero point 
energies of Hz and the isotopic molecule. The 
work of Beutler! provides the accurate value of 
36,116+6 cm or 102,800+20 cal. mole for 
D,(H2). The values for the dissociation constants 
of He, HD and D,"? are included for comparison. 
They need not be corrected since they were based 
on a value for Do(H2) fortuitously identical to 
the: present value. Dissociation constants are 
presented in Table VII. The corresponding per- 


(1938). C. Urey and D. Rittenberg, J. Chem. Phys. 1, 137 
® Beutler, Zeits. Physik. Chemie B29, 315 (1935); 
Beutler and Jiinger, Zeits. Physik 101, 304 (1936). 
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centages of dissociation for a total pressure of 
one atmosphere are given in Table VIII. 

The writer wishes to express his thanks for the 
computing assistance rendered by the following 
members of Dr. Bengt Carlson’s computing 
group: Josephine Elliott, Max Goldstein, Mar- 
garet Johnson, Nell Lane, Marian Peterson, and 
Alice Snowden. 


This paper is based on work performed at the 
Los Alamos Laboratory of the University of 
California under Government Contract W-7405- 
eng-36 and the information contained therein 
will appear in Division V of the National Nuclear 
Energy Series (Manhattan Project Technical 
Section) as part of the contribution of the Los 
Alamos Laboratory. 
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Fugacity Determinations of the Products of Detonation 


ME vin A. Cook! 
Eastern Laboratory, Explosives Department, E. I. du Pont de Nemours and Company, Gibbstown, New Jersey 


(Received June 2, 1948) 


By employing the equation of state derived from the hydrodynamic theory and observed 
detonation velocities, the thermodynamic equilibrium constant K,(T) is derived in terms of 
the concentrations of the various gases comprising the products of detonation and the ‘co- 
volume a(v), which may be evaluated for any known or assumed density of the products of 
detonation. A reiteration method is developed for solving simultaneously as many equilibria 
as is necessary to define completely the composition of the products of detonation. The applica- 
tion of the theory is illustrated by calculations for TNT which are compared with experimental 
data as well as results obtained from the approximate partial pressure method. Calculations 


are presented for several additional explosives. 


HE application of the hydrodynamic theory 

of detonation in studies of explosives re- 
quires, in general, a knowledge of the products of 
detonation from which to calculate the essential 
thermodynamic properties, e.g., the number of 
moles (”) of gas produced per kilogram of 
explosive, the heat of explosion at constant 
volume (Q), the constant volume heat capacity 
(C,), etc. The calculation of products of detona- 
tion in most explosives is not a simple matter, 
however, because the mixture may frequently 
consist of quite a large number of different gases 
interacting with each other, presumably accord- 
ing to the laws of thermodynamic equilibrium. 
It is thus necessary in an accurate treatment to 
solve simultaneously a. series of equations in- 
volving the equilibrium constants of all of the 
reactions involved in determining the equilibrium 
mixture. The problem is further complicated by 
the fact that in a thermodynamically exact solu- 
tion one must make use of the activity of the 


1 Now at the University of Utah, Salt Lake City, Utah. 


various components rather than the concentra- 
tion or partial pressure in the thermodynamic 
equilibrium constants. In the particular case of 
gaseous explosives the activity or partial fugacity 
is practically equal to the partial pressure be- 
cause the temperature and pressure conditions 
in these cases are such that the ideal equation 
of state, 


po=nRT, (1) 


applies with the required accuracy for practical 
applications. Lewis and Friauf? carried out 
calculations of products of detonations for some 
gaseous explosives taking into account the neces- 
sary equilibria, and obtained quite accurate 
correlations between calculated detonation ve- 
locities obtained by application of these data and 
experimental velocities. This work helped to 
demonstrate not only the validity of the hydro- 
dynamic theory and of the ideal gas law as 
applied to gaseous explosives, but also that 


2 B. Lewisand J. B. Friauf, J. Am. Chem. Soc. 52, 3905 
(1930). 
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chemical equilibrium exists between the products 
_ of detonation in the detonation wave of gaseous 
explosives.’ 

Regarding condensed explosives, it is not valid 
to employ partial pressures in carrying out 
calculations of the equilibrium composition of 
the product gases because detonation pressures 
are thousands of times higher than in gaseous 
explosives. As a matter of fact, in some cases the 
total fugacity may be of the order of 10‘ times 
higher than the total pressure under the condi- 
tions occurring in the detonation wave of some 
high density explosives. Despite this fact, results 
obtained by the use of partial pressures rather 
than activities have been surprisingly successful 
in many applications of theoretical calcula- 
tions.4~* In many cases, however, serious errors 
have resulted from the use of the approximation 
methods. 

Brown‘ has developed a systematic method for 
taking into account as many as eleven equi- 
librium equations although his method employs 
the partial pressure approximation. Since Brown’s 
work, the hydrodynamic theory has been worked 
out in considerable detail, making possible the 


TABLE I. Values of y and e4”” for various densities of 
products of detonation.* 


Density a y ev eu eu/2 

0 0 a 1 1 
0.5 1.089 0.95 2.59 6.69 1.62 
0.60 0.829 1.22 3.39 11.47 1.84 
0.70 0.633 1.52 4.57 20.9 2.14 
0.80 0.5 1.85 6.36 40.5 2.53 
0.90 0.396 2.22 9.21 84.8 3.03 
1.00 0.32 2.67 14.44 209 3.82 
1.10 0.261 3.13 22.9 524 4.81 
1.20 0.216 3.65 38.5 1.48 x 10 6.23 
1.30 0.179 4.20 66.7 4.45x 108 8.17 
1.40 0.153 4.79 120.3 1.45 x 104 11.02 
1.50 0.130 5.39 219 4.80 X 104 14.88 
1.60 0.112 6.05 425 1.81105 20.7 
1.70 0.095 6.78 881 7.76X105 = 29.7 
1.80 0.085 7.44 1.71 108 2.91108 41.3 
1.90 0.077 8.01 3.01 103 9.06 108 55.2 
2.00 0.070 8.60 5.43 X 108 2.95 X 107 73.7 
2.10 0.064 9.15 9.43 X 108 8.88107 97.5 
2.20 0.059 9.69 1.62 X 104 2.61108 127.7 


* These data were calculated by Mr. R. E. Lunn, Eastern Laboratory. 


3B. Lewis and G. von Elbe, Combustion, Flames, and 
— of Gases (The Macmillan Company, New York, 
“ A Schmidt, Zeits. ges. Schiess-Springstoffw. 29, 259 
6 J. Roth, Zeits. ges. Schiess-Springstoffw. 35, 193 (1940). 
Ff. Brown, Bureau of Mines Tech. Pub. No. 632. 


evaluation of the equation of state applying to 
the product gases in detonation.7—” 
Considerations discussed in reference (7), 
particularly a comparison of experimental and 
calculated detonation temperatures, support the 
assumption that an equation of state of the form 


pv=nRT+a(v)p (2) 


applies in detonation. Here a(v) is a co-volume 
factor depending only on the specific volume, v. 
The present article describes the development of 
a method of calculation of the products of 
detonation, similar to Brown’s method, but 
employing through the equation of state (2) 
partial fugacities instead of partial pressures in 
the thermodynamic equilibrium constants. The 
method is based on the process of reiteration in 
which any desired number of chemical equilibria 
may be taken into consideration. Fugacity calcu- 
lations of the products of detonation have also 
been carried out by Brinkely and Wilson® and 
others, based, however, on different equations of 
state than (2). 


THERMODYNAMIC RELATIONS 


The various chemical equilibria may be ex- 
pressed by the equacion 


(3) 


where the 7,4; refer to reactants, and the m,B, 
to products. The thermodynamic equilibrium 
constant may then be written in the form 


fBo™ i 


where the Lewis fugacity f is used for the 
activity, f; being the partial fugacity of the zth 
component of the gaseous mixture, and »; the 
coefficient of this component in (3). The total 
fugacity is related to the pressure and specific 


7M. A. Cook, J. Chem. Phys. 15, 518 (1947). 

8G. B. Kistiakowsky and E. B. Wilson, Jr., O.S.R.D. 
Report No. 69. 

*S. R. Brinkley and E. B. Wilson, Jr., O.S.R.D. Reports 
Nos. 905, 1231, 1510, 1707. 

10S. Paterson, Research 1, 221 (1948). 

" G. Morris and H. Thomas, Research 1, 132 (1947). 

” P. Caldirola, J. Chem. Phys. 14, 738 (1946). 

1G. N. Lewis and M. Randall, Thermodynamics and 
Free Energy of Chemical Substances (McGraw-Hill Book 
Company, Inc., New York, 1923). 
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Ki Ki Ke K; Ks Kes Kr Ks 
Reaction 
 4$H.+CH= CO+H,0= 2NO= 2CO= 
H,O H=3H, O=}02 CO.+H, N=—4Ne CO2+C 
2000 7.2 X10-* 9.5 4.4 X10-5 1.7 xX10-5 4.55 4.0X 26 
2400 8.3 «1075 1.4 X10-4 3.9 2.1 5.8 0.0025 112 
2800 4.7 9.0 1.83 10-3 1.2 6.6 0.0092 370 
3200 1.66 X 1073 3.64 X 1073 5.80 X 107% 4.42 10-3 7.3 7.0X 10-5 0.0244 815 
3600 4.52 10-3 1.09 x 107 0.0143 0.0129 8.1 3.0 10-4 0.0502 1410 
4000 9.9 K1073 2.54 107? 0.0292 0.0270 8.2 1.0 10-3 0.100 2360 
4400 0.028 0.049 0.053 0.052 8.3 0.166 3750 
4800 0.082 0.075 0.083 0.087 8.55 0.247 5550 
5200 0.143 0.120 0.131 8.85 
5600 0.512 0.165 9.15 
Ki Ks Kw Ku Ki Kis Kus Ku Ki 
Reaction 
0+2CO= CH= CH+H0= Ku HCN=C 
Temp. — °K 2CO; C+2H: CO+3H: Ks KK; =NH; +4N2+$He2 
2000 0.680 0.120 0.0254 4.5X 10-5 32 0.0515 
2400 0.324 0.018 0.003 8.3X10-* 10-5 0.170 
2800 0.192 0.005 7.3X10-* 6.3X10-° 10" 0.440 
- 3200 1.66104 0.0350 2.0X 10-3 0.900 
3600 0.00141 0.0905 8.6X 10-3 1.560 
4000 0.0077 0.225 0.028 97 2.350 
4400 0.104 0.406 0.132 
4800 0.470 0.641 0.327 
5200 
5600 22.00 
volume by the equation Equation (4) may now be written 
d Inf Av 
nRT: 0 
which from Eq. (2) becomes where 
d Inf nRT (11) 
=—-+a. (6) 
dp/r p In order to make use of Brown’s developments 


Integration of (6) between pressure limits of 
zero and p gives 


—=exp—— 
p PARTY, ” 
The partial fugacity of the ith component is 
related to the total fugacity by the equation 


fi=fQNi, (8) 


where JN; is the mole fraction of the ith component 
in the mixture. Multiplying and dividing the 
left side of (7) by N; and taking the »;th power 
of both sides gives 


fi vi P 
= adp. 9 


a new form of the equilibrium constant may be 
defined as follows: 


1.2181 4” 
K.-K,(——) (12 
Now 
NmRT A; T 
a 1.2181 


where a=v—a, and A; may be any of the A,’s or 
B;’s in Eq. (3) expressed in moles per 100 g of 
explosive while ‘‘a’’ has the dimensions of liters 
per kilogram as ition by Brown ( expressed 
here in atmospheres). Substituting Eqs. (12) 
and (13) in (10) gives 


expAvy, (14) 
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TABLE III. Definitions of equilibrium constants of Table IT. 


Ki =((H2)(O2)4/(H20)) 
Ke=((H2)#(OH)/H:20) -e4/?/at 
Ks3=((H)/(H2)!) -ev/*/at 
Ks=(CO)(H20) /(CO2) (He) 
Ke=((N)/(N2)4) -e4/?/at 
K7=(NO)?/(O2)(N2) 
Ks= (CO) -e%x/a 
(Ox) -e%x/a 
Kio= ((CH4)/(H2)*) -a/ev 
Kui = ((CH4)(H20) /(CO) (H2)8) -a?/e™ 


6 
Ku=(KoK 


Kis= ev/a 


‘where 

Equation (14) differs from the approximate 
partial pressure equation by the factor expAvy 
which actually is unity for the most important 
equilibrium involved in a large number of explo- 
sives, namely, the water gas equilibrium. In 
many explosives, in fact, the partial pressure and 
fugacity methods give practically identical re- 
sults for this reason. 

In view of the observation discussed in refer- 
ence (7), that @ is the same function of v for all 
explosives, a conclusion which has been verified 
in a large number of cases involving widely 
different type of explosives, y (Eq. (15)) will also 
be the same function of v for all explosives and 


may be evaluated from the a vs: v data of Fig. 1 
of reference ‘7. The values of y corresponding to 


various values of v have been evaluated by R. E. . 


Lunn from a plot of a against 1/(v—a) whose 
results are given in Table I. 

Table II gives K;(7) data taken largely from 
Brown‘ and Schmidt,‘ although some of the data 
were extrapolated from curves drawn from lower 
temperature data and should not be regarded as 
entirely accurate. 


METHOD OF CALCULATION 


The following method of calculating the 
products of detonation is modeled somewhat 
after the method of Brown® which is especially 
convenient for use with the equation of state (2). 
The equilibrium constants of Table II are in a 


form such that one need only introduce the 
factor (a) as shown by Eq. (4) to convert from 
partial pressures to concentrations. We will 
consider in the present study only explosives 
containing carbon, oxygen, hydrogen, and nitro- 
gen. Other types of explosives may be worked 
into the theory readily but, in general, it is 
unnecessary to do so, since when other atomic 
constituents are present one may generally 
decide what products of detonation are formed 
from them without taking their thermodynamic 
equilibria into consideration. 

The atom balances of the explosive composi- 
tion may be expressed as follows: 


Co=CO2+CO+CHit+HCN 


+2C2N2+:--, (16a) 


+HCN+3NH3;+::-, (16b) 


+NHs+---, (16c) 


Oo = 
+H,0+OH+:--. (16d) 


Four key equations may be derived from (16) 
by solving for CO, He, Ne, and a new oxygen 
balance expression, giving 


CH,+C+HCN +2C2N2.+ eee 
Co 
(17a) 
H+0OH+4CH.+HCN+3NH; 
1 
Ho 
(17b) 
N =| 
Ny 
(17c) 


H 


OH 
(17d) 


| 

H)=2H.0+2H.+H+OH+4CH, 

I 

fr 

tic 

de 

fi 

cal 

A= 

| 

apy 

wit 


1) 
6) 


b) 


/c) 


id) 
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Here A is the oxygen balance, 


Ho 
Cox 2 
B,= ’ Bz= ’ 
1+x 


and x =CO/COs». The derivatives of these func- 
tions are given by Brown. The dimensions of (a) 
used by Brown was liters/kg, while the atom 
balances were expressed in g atoms/100 g of 
explosive. This convention is retained here. In 
terms of the K; of Table III, Eqs. (17a), (b), (c) 
may be written — 


K12(CO) 
a’Co 


co=s,(1 
(18a) 


Kis 
x 
He. =B, 1 


4K (CO) (H 2) 
Hoa? 


Ho 


N 
Ns ~ 
2 Noe’? xNoew/? 


(18c) 


(18b) 


The values of (a) and e4” may be determined 
from the conditions under which the computa- 
tions are to be made, i.e., the known or assumed 
density. One may then proceed as follows: 


Step I. A trial temperature is taken. It is convenient 
first to estimate this temperature from some approximate 
calculation, e.g., by writing approximate chemical equations 
for the explosive reaction, or one may make a rough 
calculation assuming only the water gas equilibrium (Ks) 
to be important. 

Step II. Find the value of x (by a process of successive 
approximations) by means of Eqs. (18a), (b), (c) together 
with whatever auxiliary equilibrium equations (see Table 
III) are needed to make a complete set, such that Eq. 


(17d) gives the same value of A as the equation 


A= (19) 


Step III. When the value of x is found such that Eqs. 
(17d) and (19) give identical results for the oxygen balance, 
one may compute from the composition obtained the heat 
of explosion Q and the average heat capacity C, between 
To and T which will then allow one to derive T for adiabatic 
explosion, by the equation 


T= Q/C.+To (20) 


or, when dealing with the detonation wave conditions, by 
the Eq. (16), reference 7, namely, 


(21) 


Step IV. Repeat steps II and III until the temperature 
calculated from (20) or (21) agrees with that used in the 
computations in step II. 


In case of large oxygen deficiency, the method 
of calculation may conveniently be altered to 
take into account the relative importance of the 
coke oven equilibrium, (Ks). The selection of a 
particular x then immediately fixes the values of 
both CO and COs. The remaining Cy may be 
then distributed: (a) as CH, which may be 
computed from Eq. (7), Table III; (b) the carbon 
remaining may be considered to go into HCN 
until Eq. (16), Table III, is satisfied; (c) any Co 


TABLE IV. Products of detonation of TNT from theo- 
retical fugacity and partial pressure methods, and Haid 
and Schmidt’s experimental method. 


B. “Detonation” 


A. “Explosion” Exp. wave state 
state Haid partial 
I. =1.59 g/cc fugacity partial and fugacity pressure 
products method pressure Schmidt method method 
CO (moles/kg) 9.4 15.9 9.39 1.7 15.5 
COz 6.8 3.2 5.47 9.7 3.4 
He 0.6 2.9 1.63 0.07 2.7 
H:0 3.4 4.3 6.09 3.1 4.2 
Ne 5.0 $2 4.4 4.0 
CHa 2.2 1.25 0.42 0.7 
NH; 1.1 0.13 1.5 3.4 0.1 
HCN 2.0 2.65 0.32 1.0 5.3 
Carbon 10.4 7.38 14.6 16.3 5.9 
nm (moles gas/kg) 30.5 35.6 24.5 35.8 
@) 1120 860 1105 770 
Co 0.342 0.317 0.338 0.326 
Ts 3420 2975 4150 (T2) 3300 (T2) 
II. =1.0 g/cc 
co 16.8 18.2 17.0 16.3 23.1 
COz 3.5 2.0 2.84 2.8 0.9 
He 2.5 4.3 2.49 1.5 44 
H:0 2.6 3.2 3.77 4.5 1.3 
Ne 5.5 Sa 4.75 4.6 5.1 
CHa 1.9 12 0.1 2.3 1.8 
NH; 0.8 0.1 2.85 0.7 0.2 
HCN 1.4 1.72 0.47 3.2 2.9 
Carbon a 7.3 10.0 5.9 a4 
n (moles gas/kg) 35.1 36.6 35.9 39.7 
QO 865 7380 870 660 
Co 0.317 0.304 0.332 0.303 
T3 3000 2830 3340 (T2) 2280 


» 
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TABLE V. Explosion properties of various types by fugacity and partial pressure methods. 
(a=fugacity method; b=partial pressure method.) 


PETN (density =1.6) Cyclonite (1.6) LNG (1.6) Tetryl (1.6) 60% St. Dyn. (1.32) 
Products (a) (b) (a) (b) (a) (b) (a) (b) (a) (b) 
CO (moles/kg) 4.7 6.5 8.0 10.1 1.1 2.6 10.6 18.1 6.7 7.0 
CO, 11.1 9.3 5.5 3.4 12.1 10.6 6.6 2.9 8.3 8.0 
Hz 0.5 0.9 1.3 35 0.1 0.3 0.5 2.9 0.9 a4 
H,0 10.9 11.2 8.0 9.8 10.8 10.2 a2 4.0 10.0 10.1 
Ne 5.9 6.1 12.1 13.5 6.4 6.0 7.0 7.0 4.9 4.9 
OH 0.1 0.6 — 0.3 0.2 0.9 —_— — —_— 0.1 
NO 0.1 0.5 0.5 0.01 0.02 
n 34.1 35.4 37.7 40.7 32.6 Se 32.3 38.3 31.1 32.3 
g 1420 1310 1315 1250 1420 1240 1013 910 985 975 
3 4770 4420 4530 4330 4780 4350 3650 3470 3690 3590 


remaining may be considered to go to free carbon. 
This calculation is necessarily approximate in 
view of the absence of heat data for other 
possible products of detonation involving carbon. 


RESULTS 


Table IV presents calculated results for TNT 
at densities of 1.59 g/cc (cast) and 1.0 g/cc for 
(A) “adiabatic” or ‘‘explosion”’ conditions, and 
(B) for conditions occurring at the Chapman- 
Jouguet plane. Also included in Table IV are 
experimental results obtained by Haid and 
Schmidt" in a method designed to “freeze” the 
equilibrium corresponding to the “explosion” 
state. A careful examination of the comparative 
results between experimental and calculated data 
reveals some definite consistencies if one inter- 
prets the experimental data on the basis that the 
experimental method does not entirely ‘‘freeze”’ 
the equilibria at those corresponding to “explo- 
sion”’ state. The excellent agreement in CO is 
probably not fortuitous; this type of agreement 
is quite general as has been found in numerous 
other comparisons by the author of calculated 
and experimental data. In fact, one has theo- 
retical justification for the consideration that a 
moderate shift in the water gas equilibrium from 


4 A. Haid and A. Schmidt, Zeits. ges. Schiess-Spring- 
stoffw. 26, 253, 293 (1931). 


- the ‘‘explosion”’ state (by decreasing 7) should 


not affect the CO concentration appreciably 
when CH, and HCN are present in the gas 
mixture. The failure to freeze the water gas 
equilibrium completely is reflected, however, in 
low CO, and high H,0O in the experimental data. 
CH, will exist in relatively large amounts only 
because the pressure is very high. One would 
expect the pressure to drop more rapidly than 
the temperature, -after detonation, because of 
the co-volume effect, so that CH, would tend 
to decompose into C and He. HCN would tend 
to decompose as the temperature decreases. Both 
effects would result in high C and He: concentra- 
tions as observed. Since NH; is very sensitive 
to Hz concentrations, one accounts for the excess 
NH; found experimentally on the basis that it 
formed during the decomposition of CH, and 
HCN. 

Table V lists results calculated by the fugacity 
and partial pressure methods. When the temper- 
ature is high and/or the oxygen balance is quite 
negative the two calculations lead to quite 
different results. However, as the oxygen balance 
approaches zero the differences become less, and, 
if the “explosion” temperature is also low the 
difference between the partial pressure and the 
fugacity methods is negligible at or near zero 


O.B. 


J 
1 
( 


F 
cl 
al 
at 
of 
| Ve 
cu 
na 
of 
co 
me 
ar 
na 
mi 
| pre 
ate 
co! 
inc 
vay 
: 
( 
car 
pre 
as | 
Bas 
tru 
at 1 
Ste: 
Bas 
con 
rela 


THE JOURNAL OF CHEMICAL PHYSICS VOLUME 16, NUMBER 11 NOVEMBER, 1948 


Letters to the Editor 


HIS section will accept reports of new work, provided 

these are terse and contain few figures, and especially 
few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
Contributions to this section should not exceed 600 words in 
length and must reach the office of the Managing Editor not 
later than the first of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
authors. The usual publication charge ($4.00 per page) will 
not be made and no reprints will be furnished free. - 


The Vapor Pressure and Heat of 
Sublimation of Graphite 


L. H. LonGc 
Washington Singer Laboratories, University College, Exeter, England 
August 31, 1948 


N arecent communication under the above title, Brewer, 
Gilles, and Jenkins! describe an equilibrium effusion 
method for measuring the vapor pressure of graphite. 
From the results of experiments below 3000°K they con- 
clude that the vapor reaches a pressure of one atmosphere 
at 4630°K, and that the heat of sublimation is 170 kcal./g- 
atom. 

However, apart from the fact that several indirect lines 
of evidence discussed elsewhere? will not admit of the latter 
value, these figures encounter a number of serious diffi- 
culties when compared with much evidence of a more direct 
nature, as pointed out below: 

(i) Setting aside experiments involving the temperature 
of the positive crater in the carbon arc, which do not give 
consistent results, there are only three other sets of experi- 
ments in which equilibrium temperatures and pressures 
are measured for graphite. The results of two of these, 
namely, those of Thiel and Ritter® (first two of the three 
methods described) and Basset,‘ provide temperature- 
pressure data closely supporting the figure 125 kcal./g- 
atom for the heat of sublimation (the lowest alternative 
considered by Brewer et al.). The third, due to Ribaud and 
Bégué,® is in poor agreement with the other two and 
indicates a lower vapor pressure; but in all three cases the 
vapor-pressure data are higher than those of Brewer et al. 
by several powers of ten. 

(ii) Basset* found the coordinates of the triple point of 
carbon to be 4000°K and 100 atmospheres, below which’ 
pressure the carbon sublimed without melting. Even if, 
as Brewer et al. suppose, the temperature measurements of 
Basset are seriously in error, the same is not likely to be 
true of his pressure measurements. Further, the pressure 
at the triple point has been independently determined by 
Steinle* who also finds that it is 100 atmospheres. Both 
Basset and Steinle provide photographic evidence for their 
conclusions. But, according to the temperature-pressure 
relation of Brewer et al., the vapor pressure reaches 100 


atmospheres above 6000°K. This would place the melting 
point of graphite at an impossible figure well above those 
of the refractory metal carbides, whereas, for example, 
graphite can be fused in a crucible of tantalum carbide 
which has a reputed melting point of 4150°K. Other 
experimental figures for the melting point of graphite are 
variously given as 3800,7 3900,® 3760,° and 3845°K.” It thus 
does not appear that Basset’s figure is an underestimate, 
nor can it be presumed that all of the various investigators 
cited here have erred in their temperature measurements by 
such an enormous and approximately constant amount. 

(iii) If the vapor pressure of graphite reaches 1 atmos- 
phere at 4630°K, it is surprising that temperatures reported 
for the positive crater of the carbon arc, which are generally 
regarded to be limited by the vapor pressure, are never 
sensibly above 4000°K, even under elevated pressures. 

(iv) Apart from absolute temperature measurements, 
Mott" has determined the relative volatilities of a number 
of substances by a variety of methods. His placing of 
carbon offers support to the measurements of Basset and 
of Thiel and Ritter, but is decidedly out of harmony with 
those of Brewer et al. 

The consensus of opinion is thus that graphite possesses 
a very much higher vapor pressure and consequently 
lower heat of sublimation than Brewer et al. affirm. It is 
relevant to note in conclusion that these authors in their 
experiments with replaceable graphite tubes found that 
rapid volatilization occurred which limited the life of the 
tubes. This observation is scarcely compatible with their 
own temperature-pressure relation, which in all cases 
necessitates pressures well below 0.1 mm at the tempera- 
tures of their experiments. 

1L. Brewer, P. W. Gilles, and F. A. Jenkins, J. Chem. Phys. 16, 
797 (1948)... 

2 L. H. Long and R. G. W. Norrish, Proc. Roy. Soc. A187, 337 (1946); 
L. H. Long, forthcoming paper. 

3 A, Thiel and F. Ritter, Zeits. f. anorg. Chemie 132, 125 (1923). 

4]. Basset, J. Phys. Radium (vii) 10, 217 (1939); Chimie et Industrie 
46, No. 3b, 7 (1941); Brennstoff-Chemie 23, 127 (1942). 

5G. Ribaud and J. Bégué, Comptes Rendus 221, 73 (1945). 

6 H. Steinle, Zeits. f. angew. Mineralogie 2, 344 (1940). 

7K. Fajans and E. Ryschkewitsch, Naturwiss. 12, 304 (1924); 
E. Ryschkewitsch, Zeits. f. Elektrochemie 31, 54 (1925). 

8 A. Hagenbach and W. P. Liithy, Naturwiss. 12, 1183 (1924). 

9H. Alterthum, W. Fehse, and M. Pirani, Zeits. f. Elektrochemie 31, 


313 (1925). 
10 E, Ryschkewitsch and F. Merck, Zeits. f. Elektrochemie 32, 42 


(1926). 


1 W. R. Mott, Trans. Am. Electrochem. Soc. 34, 255 (1918). 


Elution Equations for Adsorption and Ion 
Exchange in Flow Systems 


NeEvIN K. HIESTER AND THEODORE VERMEULEN 
University of California, Berkeley, California 
September 7, 1948 


LUTION, the process of separating from a granular 
E adsorbent bed one or more components previously 
adsorbed upon it, is an essential step in any regenerative 
cycle of ion exchange, chromatography, or gas adsorption. 
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The typical ‘‘elution bell’ i.e., the concentration-time curve 
which displays a maximum, has been fitted approximately 
by the stepwise-discontinuous method of equilibrium 
stages.1-$ For the special case where the bed is initially 
fully saturated, Thomas‘ has provided exact kinetic 


Fic. Concentration curves for a saturating component in an 
eluting uid stream passing a given cross section of granular 
calculated from Eq. (4). 


equations which include the result of Boyd, Myers, and 
Adamson.* By an extension of Thomas’s method, we have 
now derived general kinetic equations for the elution bell 
for any fractional extent of previous saturation. 

Our study has shown, incidentally, that Thomas’s equa- 
tions for saturation of an initially empty bed® contain the 
equations proposed by several other investigators. For the 
concentrations of saturating component in the fluid (c) 
and on the bed (p), Thomas obtained 


£____ (1) 
Co rs) rt)’ 


THE EDITOR 


rs) (2) 
pa ot, 1s, 71)’ 


where we utilize his notation except for r, s, t, and p. which 
have the values given in Table I. Only when r=1 can c/co 
and p/p. be evaluated exactly from the heat-transfer 
charts of Furnas.” The functions F(c, p) in Table I contain 
in each case the rate constants for a second-order forward- 
second-order reverse reaction (typified by exchanging ions 
of equal valence). For kinetically simpler cases, the con- 
stant factors (ka) or (R2co) must be replaced in the tabu- 
lated terms by whatever rate constants are appropriate. 

Consider now that a fluid containing a single saturating 
component at concentration ¢» has passed through the bed 
during a time interval, tsat, which corresponds to a par- 
ticular value of y and thus of #;;(=7;;). At this instant, an 
eluting fluid (the same fluid, but without saturating com- 
ponent) begins to pass through the bed at an unchanged 
rate of flow. (For ion exchange, the fluid is assumed to 
contain an eluting component identical with the species on 
the bed prior to saturation; its concentration will tempo- 
rarily be taken equal to that of the saturating component, 
Co, in the preceding fluid.) 

The equations for elution will depend, first, upon the bed 
concentration at each point at passage of the (assumedly 
sharp) new fluid front. Thus at y’ = (r’—1/R)/m=0, where 
7’ is time reckoned from the start of regeneration, we have 


; ¢(T, rs) (3) 
Po 15) +I0(2(rsT)4)+ rT)" 


In addition, since the component of interest is no longer 
entering the bed, c=0 at x =0. Thomas’s conservation and 


TABLE I. Rate, distance and time parameters for Thomas equations. 


Previous 
investi- 
Reaction order Rate, F(c, p) peo Dimensionless parameters gators 
Forward Reverse 
(i) Gj) rij Si bij 
Surface reaction 
controlled 
kic(a—p) a ke/ki mkiax mkicoy a 
—k2p(co—c) 
hie(a—p) 
2 0 kic(a—p) a 0 mkax i mkicoy c-f 
(Rkia)e kia 
1 1 —(kace)p 1 m(kia)x m(k2co)y b,e 
1 0 (Ria)c 1 m(kia)x 
Diffusion 
controlled 
1a 1a 1 
1 1 ko 1 mk ax mk py h-j 
* See reference 6. 
b See reference 4. 
e G. S. Bohart and E. Q. Adams, J. Am. Chem. Soc. 42, 523 (1920). 
4L. G. Sillen, Arkiv. f. Kemi. Min. Geol. A22, No. 15 (1946). 
76 E. Walter, J. Chem. Phys. 13, 332 (1945). 
f Danby, Davoud, Everett, Hinshelwood, and Lodge, J. Chem. Soc. 68, 918 (1946). 
® Drew, Spooner, and Douglas: D. M. Klotz, Chem Rev. 39, mal (1946), reference 10. 
b See reference 5. 
iO. A. Hougen and W. R. Marshall, Chem. Eng. Progress 43, 4 (1947). 
i R. H. Beaton and C. C. Furnas, Ind. Eng. Chem. 33, 1500 (1941). 
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rate conditions are then applied, with replacement of y by 
y’, and lead to the equations: 


eTLo(T+?’, rs) +Io(2(rsT+rst’)*)] 
<= —L[elt’, rs)-+ (4) 
eT rs) +-1o(2(rsT +rst’)4) 
+¢(s, rT +rt')]+[e(rt’, s)— et’, rs)] 


2. e~*[o(T+’, rs) ]—Lelt’, rs)] (5) 
bo rs) 1o(2(rsT +rst’)#) 
+¢(s, s)— elt’, rs)] 


where r and s remain as before, and ¢’;; corresponds to the 
values of ¢;; in Table I with y replaced by y’. Figure 1 
presents typical elution curves as calculated from these 
equations, with r;; and s;; constant. At large T the equa- 
tions become identical with those of Thomas,‘ as is indi- 
cated in the figure; for instance: 


g(rs, t’) (6) 


As the new fluid front travels down the bed, the discon- 
tinuity in c/co decreases and rapidly becomes negligible; 
in actual practice, it may be obscured by diffusional effects 
even at small distances. For ion exchange, the above equa- 
tions may be extended to regenerant concentrations other 
than Co (i.e., Co’) by replacing co by ¢o’ (except in the calcu- 
lation of T which is a “‘memory”’ term); in this case a dis- 
continuity is always encountered at the front. Figure 1 
shows also how regeneration time is reduced by doubling 
the concentration level of regenerant. 

The curve for p is continuous, but otherwise of the same 
general shape as for c; its maximum represents a true 
equilibrium condition, at the time that net reaction at that 
particular cross section finally shifts from saturation to 
elution. 

Where r=1, Eqs. (4) and (5) take on an especially 
simplified form: 


T(t’, rs) +Io(2(rst)t)], (7) 


+r) o( T+’, rs) [(t’, (8) 


The two terms in each of the above equations represent 
pairs of values in the Furnas charts that are displaced 
along the ¢’ scale by a value of 7. 

Saturation and elution equations have also been derived 
for cases involving more than one “‘memory” term. We 
are preparing tables of ¢ to facilitate exact evaluation of 
the above equations, and also are charting solutions of 
Eq. (1) for values of r other than unity. Complete deriva- 


tions will be reported and‘ applications discussed in a 
forthcoming paper. 


1A. V. P. Martin and R. L. M. Synge, Biochem. rer aes (1941). 
2S. Stene, Arkiv. f. Kemi. Min. Geol. 18, No. 18 (1945). 

asa). . Mayer and E. R. Tompkins, J. Am. Chem. Soc. 69, 2866 
4H. . Thomas, Annals N. Y. Acad. Sci. 49, 161 (1948). 


5G. E. Boyd, .& Myers, Jr., and A. W. Adamson, J. Am. Chem. 
Soc. 69, 2849 (19479. 
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Absorption Spectrum of Hexavalent Plutonium 


R. H. Betts AND B. G. HARVEY 


Atomic Energy Project, National Research Council, 
Chalk River, Ontario, Canada 


September 7, 1948 


HE chemical similarity between uranium and plu- 

tonium now appears to be firmly established.'? 
Recently, Kasha* has used this analogy to explain certain 
features of the absorption spectrum of PuO,t* in 0.1N 
perchloric acid. Now the absorption spectrum of UO,*+ 
is well known;‘ it contains eleven regularly spaced bands 
in the region 3400-5000A. These bands are believed to 
arise from symmetrical vibrations in the U—O bonds of 
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Fic. 1. Absorption spectrum of Pu(VI). 
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Fic. 2. Absorption spectrum of Pu(VI) in HNOs solution. 
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the UO,*t* ion. Kasha found a series of four regularly 
spaced bands for PuO.** in the region 3900-4300A, and 
ascribed them to symmetrical vibrations in the Pu—O 
bonds. In addition, an exceedingly sharp band at 8330A 
was found, which was attributed to electronic transitions 
within the 5f electron shell. The UO.** spectrum has no 
corresponding sharp lines; this is to be expected, since 
uranium contins no 5f electrons in this valence state.? 

We have measured the absorption spectrum of PuO.+* 
in 0.9N nitric acid solution, using a Beckman quartz spec- 
trophotometer. The spectrum over the region 3600—-10,000A 
is shown in Fig. 1. This result is very similar to that de- 
scribed for PuO.** in 0.1N perchloric acid. That portion 
of the spectrum ascribed by Kasha to vibration is shown on 
an enlarged scale in Fig. 2. Four regularly spaced bands 
may be distinguished in the region 3900-4300A (indicated 
by arrows). Starting from the ultraviolet end, the suc- 
cessive frequency intervals for these bands are 630, 670, 
and 670 cm™, giving an average of 657 cm™. This agrees 
within experimental error with the value of 708 cm™ 
found by Kasha. 

A qualitative comparison of the uranyl and plutonyl 
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spectra at higher pH values tends to confirm this distinc- 
tion between the vibrational and electronic portions of the 
Pu(VI) spectrum. For U(VI), Sutton‘ has shown that as 
the pH is raised above 2.0, the banded structure gradually 
disappears, and is replaced by a region of more or less 
continuous absorption below 5000A. We have determined 
the absorption spectrum of Pu(VI) in dilute ammonium 
hydroxide solution. (Ammonium plutonate is insoluble, in 
general, in dilute ammonium hydroxide solution, but we 
have found on occasion that a solution of Pu(VI) deepens 
to an orange color on addition of excess ammonia without 
the appearance of a precipitate.) The result is given in 
Fig. 1. The peak at 8330A is broadened slightly, and shifted 
to longer wave-lengths, without a marked change in the 
extinction coefficient. However, the band structure has 
been replaced by a region of general absorption with much 
higher values of the extinction coefficients. 

1 Harvey, Heal, Maddock and Rowley, J. Chem. Soc. 1010 (1947). 

2B. G. Harvey, Nucleonics 2, 30 (1948). 

3M. Kasha, MDDC-591, U. S. Atomic Energy Commission, Oak 
Ridge, Tennessee (and references therein). 


4J. Sutton, Report No. 1612, National Research Council of Canada, 
Ottawa, Canada. 
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